Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,935,301 C→A 5.5% intergenic (‑54/+1092) USA300HOU_1798 ← / ← hsdS2 hypothetical protein/type I site‑specific deoxyribonuclease specificity subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,935,3010CA 5.5% 215.0 / 2.6 75intergenic (‑54/+1092)USA300HOU_1798/hsdS2hypothetical protein/type I site‑specific deoxyribonuclease specificity subunit
Reads supporting (aligned to +/- strand):  ref base C (29/39);  new base A (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

GGATTTGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAATATTCGGAATA  >  CP000730/1935168‑1935441
                                                                                                                                      |                                                                                                                                             
ggATTTGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCgg                                                                                                                                        <  2:851688/141‑1 (MQ=255)
ggATTTGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAAGTCTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCgg                                                                                                                                        <  2:836960/141‑1 (MQ=255)
   tttGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTaa                                                                                                                                     >  2:317315/1‑141 (MQ=255)
           aTTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAATCTTTAATCGATAATCTTCTTCCTTTTTATTTGTTAAATACTCTTTTTTATAGTGTCCTCTATTTTTTTCTTTTTTCACTTTGCGGTAACTTGaata                                                                                                                             >  1:683364/1‑141 (MQ=255)
           aTTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACAAGTTGCGGTAACTTGaata                                                                                                                             >  1:147174/1‑141 (MQ=255)
           aTTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACAAGTTGCGGTAACTTGaata                                                                                                                             >  1:664715/1‑141 (MQ=255)
               ttATTTTTAATATTATTATTTTGATTACCGGTAAATATAAACTTTAATCGATAACCTTCTTCCAGTTTATGTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTTCTTTTTACACGTTGCGGTAACTTGaatactaa                                                                                                                         >  1:158667/1‑141 (MQ=255)
                     ttAATATTATTATTTTGATTTCCTATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGaataatgatgtcga                                                                                                                   <  2:367855/141‑6 (MQ=255)
                     ttAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTAt                                                                                                                   <  2:855297/141‑1 (MQ=255)
                     ttAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTAt                                                                                                                   <  1:569636/141‑1 (MQ=255)
                         tattattatTTTGTTTTCCAATAAATATAAACTTTATTCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGt                                                                                                               <  1:597832/141‑1 (MQ=255)
                            tattattTTGATTTCCAATCAATATAAACTTTACTCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatt                                                                                                            <  2:443928/141‑1 (MQ=255)
                            tattattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGTTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtctt                                                                                                            >  2:790031/1‑141 (MQ=255)
                             attattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatta                                                                                                           >  2:203784/1‑141 (MQ=255)
                             attattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatta                                                                                                           <  1:844171/141‑1 (MQ=255)
                             attattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatta                                                                                                           <  2:426559/141‑1 (MQ=255)
                             attattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatta                                                                                                           >  1:970350/1‑141 (MQ=255)
                                  tttGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCtat                                                                                                      <  1:953838/141‑1 (MQ=255)
                                                  tataAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGt                                                                                     >  1:845279/1‑141 (MQ=255)
                                                  tataAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGt                                                                                     <  1:1044804/141‑1 (MQ=255)
                                                   ataAACTTTACTCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTTCACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTg                                                                                    <  2:1010581/141‑1 (MQ=255)
                                                   ataAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTg                                                                                    <  2:1058316/141‑1 (MQ=255)
                                                      aaCTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCa                                                                                  >  1:597518/1‑141 (MQ=255)
                                                      aaCTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCAt                                                                                 >  1:104641/1‑141 (MQ=255)
                                                      aaCTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCAt                                                                                 <  1:778650/141‑1 (MQ=255)
                                                       aCTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCAtt                                                                                <  2:891400/141‑1 (MQ=255)
                                                        cTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCAttt                                                                               <  2:941157/141‑1 (MQ=255)
                                                          ttAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCAttttt                                                                             <  1:500001/141‑1 (MQ=255)
                                                           tAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCAtttttt                                                                            <  2:272946/141‑1 (MQ=255)
                                                             aTCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTTTTTTTTTTATTTTTTTGTTGTTTTAATAATGATTTTATTCTTATTTTTTATTATATGTAATGCTGGTGTCtttttttt                                                                          >  2:912886/1‑140 (MQ=255)
                                                              tCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTTTTTTTTTTTTTTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGt                                                                         >  2:11372/1‑141 (MQ=255)
                                                                gATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATtttttttttt                                                                       >  2:590348/1‑141 (MQ=255)
                                                                gATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGttt                                                                       <  2:406526/141‑1 (MQ=255)
                                                                  tAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGttttt                                                                     >  2:1059899/1‑141 (MQ=255)
                                                                      cttcttCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTttatt                                                                 <  1:985136/141‑1 (MQ=255)
                                                                         cttcTATTTTATTTGTTAAATACTCTTTTTTTTTTTGTTTTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTttttgttttttgttttt                                                              >  2:951927/1‑138 (MQ=255)
                                                                          ttcCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGa                                                             <  2:1079245/141‑1 (MQ=255)
                                                                            ccATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGAtt                                                           <  1:479823/141‑1 (MQ=255)
                                                                            ccATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGAtt                                                           <  1:251018/141‑1 (MQ=255)
                                                                              aTTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATtct                                                         >  1:122138/1‑141 (MQ=255)
                                                                                tttATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCa                                                       <  2:970350/141‑1 (MQ=255)
                                                                                   aTTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGATAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTTCATTTTTTTGTTTTTTATTATGATTCTCaatt                                                    >  2:486158/1‑139 (MQ=255)
                                                                                      tGTTAAATACTCTTTTTTTTTTTTTCTTCTAATTTTTGCGGGTTATACGTTGTGGTAACTTGAGTAATGATTTTATTGGGATTATCTATAAGATCAAATTCTGTTGCATTTTTTTGTTTTTtattatttttttcttgtgat                                                 >  2:52957/1‑127 (MQ=255)
                                                                                      tGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTTTTTTTTTTTTTTTTTTGTAACTTGAATAATAATTTTATTCGTATTATCTATATGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATTATtcttattttat                                                 >  2:404856/1‑133 (MQ=255)
                                                                                      tGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTTTTTTTTTTTATGATTCTCAAGTGaa                                                 >  2:902013/1‑141 (MQ=255)
                                                                                       gTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTTTTTTTTTTTCTTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCaatttatt                                                >  2:400070/1‑135 (MQ=255)
                                                                                            aTACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCa                                           >  2:257787/1‑141 (MQ=255)
                                                                                              aCTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCata                                         <  2:95962/141‑1 (MQ=255)
                                                                                               ctctTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCatat                                        <  2:868838/141‑1 (MQ=255)
                                                                                                 ctTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATAc                                      <  2:276394/141‑1 (MQ=255)
                                                                                                 ctTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATAc                                      <  1:163884/141‑1 (MQ=255)
                                                                                                    tttttAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACaa                                    <  2:997564/141‑1 (MQ=255)
                                                                                                    tttttAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACaa                                    <  1:138688/141‑1 (MQ=255)
                                                                                                       ttAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAatat                                 >  1:780349/1‑141 (MQ=255)
                                                                                                       ttAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAatat                                 <  2:504282/141‑1 (MQ=255)
                                                                                                       ttAGAGTGTCCTCTATTTTTTGCTTTTTACAAGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAatat                                 <  2:147174/141‑1 (MQ=255)
                                                                                                       ttAGAGTGTCCTCTATTTTTTGCTTTTTACAAGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAatat                                 <  2:664715/141‑1 (MQ=255)
                                                                                                        tAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTTTTATTCTTAAGTGAATACCCATATACattatt                                >  2:379712/1‑141 (MQ=255)
                                                                                                         agagTGTCCTCTATTTTTTGCTTTTTACACGTTGTGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATAttt                               <  2:646367/141‑1 (MQ=255)
                                                                                                         agagTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATAttt                               <  1:19785/141‑1 (MQ=255)
                                                                                                         agagTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATAttt                               <  1:826141/141‑1 (MQ=255)
                                                                                                           agTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTaa                             <  2:2760/141‑1 (MQ=255)
                                                                                                           agTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTaa                             >  1:719606/1‑141 (MQ=255)
                                                                                                            gtgtCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAAc                            >  1:530269/1‑141 (MQ=255)
                                                                                                              gtCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAg                          <  2:1028286/141‑1 (MQ=255)
                                                                                                                      ttttttGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATTATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaa                  >  2:427676/1‑141 (MQ=255)
                                                                                                                      ttttttGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaa                  >  1:820663/1‑141 (MQ=255)
                                                                                                                       tttttGCTTTTTACTCGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATTAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATTATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaaa                 <  2:154290/141‑1 (MQ=255)
                                                                                                                        ttttGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaaaa                >  1:255304/1‑141 (MQ=255)
                                                                                                                        ttttGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaaaa                >  1:262495/1‑141 (MQ=255)
                                                                                                                           tGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaaaaaaa             <  1:971137/141‑1 (MQ=255)
                                                                                                                            gCTTTTTACACGTTGAGGTAGCTTGAATAATAATTTTATTCGTATTATCTAAAAGATCATATGCTAGTGCATTTTATTATTTTTTATTAGGATTCGCAAGTAAAAACCCATCAAAAATAATTAACAGGTATCTAAAAAAat            <  2:923656/141‑1 (MQ=255)
                                                                                                                             cTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAata           >  1:453790/1‑141 (MQ=255)
                                                                                                                              tttttACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAatat          >  1:931675/1‑141 (MQ=255)
                                                                                                                                 ttACTCGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTATAAAAAAATATTCg       >  2:636757/1‑141 (MQ=255)
                                                                                                                                    cacGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAATATTCGGaa    <  1:714218/141‑1 (MQ=255)
                                                                                                                                     acGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAATATTCGGAat   >  1:506507/1‑141 (MQ=255)
                                                                                                                                      cGTTGTGGTAACTTGAATAATAATTTTATTCGTATTATCTACATATTCATTTCTTTGGTGTTTTTGGTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAATAGGTCTCTAAAAAAATATTCGGAata  <  1:431954/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                             
GGATTTGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAATGCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAATATTCGGAATA  >  CP000730/1935168‑1935441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: