Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 662,685 A→C 5.8% intergenic (‑43/+67) ybeD ← / ← dacA UPF0250 family protein/D‑alanyl‑D‑alanine carboxypeptidase (penicillin‑binding protein 5)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913662,6850AC 5.8% 387.6 / 9.9 132intergenic (‑43/+67)ybeD/dacAUPF0250 family protein/D‑alanyl‑D‑alanine carboxypeptidase (penicillin‑binding protein 5)
Reads supporting (aligned to +/- strand):  ref base A (67/57);  new base C (6/2);  total (73/59)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 135 total aligned reads displayed.

ACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCAACA  >  NC_000913/662539‑662834
                                                                                                                                                  |                                                                                                                                                     
aCCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAgg                                                                                                                                                    >  2:258725/1‑150 (MQ=255)
aCCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAgg                                                                                                                                                    >  2:1282103/1‑150 (MQ=255)
 ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGc                                                                                                                                                   <  1:1710580/150‑1 (MQ=255)
  cACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTGCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCg                                                                                                                                                  >  1:1681346/1‑150 (MQ=255)
   aCTTCAACCACCTGATCAACCAGCTCAGGTAAAGCCTGCCCCATAACTTTGAAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAAATAAATAGCAACTCCCGCCAGAAGGCgg                                                                                                                                                 >  1:1187615/1‑150 (MQ=255)
   aCTTCAACCACCTGATCAACCAGCTCAGGCAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCgg                                                                                                                                                 <  1:904268/150‑1 (MQ=255)
     ttCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTACCGCCAGAAGGCGGGa                                                                                                                                               >  2:132122/1‑150 (MQ=255)
      tCAACCACCTGATCAACCAGCTCAGGTAACTCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCTTCTGGCGGGAg                                                                                                                                              <  2:342217/150‑1 (MQ=255)
        aaCCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCTTCTGGCGGGAGtt                                                                                                                                            <  2:1012949/150‑1 (MQ=255)
         accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttt                                                                                                                                           <  2:232848/150‑1 (MQ=255)
            accTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttttt                                                                                                                                        >  1:1610505/1‑150 (MQ=255)
             ccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTTAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATTTCAGCTCAGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                                                                                                                                       >  1:1647109/1‑150 (MQ=255)
             ccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                                                                                                                                       <  2:1561081/150‑1 (MQ=255)
               tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACt                                                                                                                                     <  2:461768/150‑1 (MQ=255)
                 aTCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGa                                                                                                                                   <  1:1818506/150‑1 (MQ=255)
                 aTCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGa                                                                                                                                   <  1:1455459/150‑1 (MQ=255)
                   cAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                                                                                                                                 >  2:561725/1‑150 (MQ=255)
                      ccAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCtt                                                                                                                              <  1:149213/150‑1 (MQ=255)
                      ccAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCtt                                                                                                                              >  1:966364/1‑150 (MQ=255)
                       cAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTa                                                                                                                             <  2:549525/150‑1 (MQ=255)
                        aGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAg                                                                                                                            <  1:16093/150‑1 (MQ=255)
                        aGCTCAGGTAACGCCTGCCCCATAACTATGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAATTAGAAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAg                                                                                                                            <  2:1179050/150‑1 (MQ=255)
                         gCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTGCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                                                                                                                           <  1:744182/150‑1 (MQ=255)
                          cTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                                                                                                                          >  2:1395955/1‑150 (MQ=255)
                            cAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGtata                                                                                                                        >  2:1540127/1‑150 (MQ=255)
                               gTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTATTTTAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTTTATgg                                                                                                                     <  2:617294/150‑1 (MQ=255)
                                 aaCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgggg                                                                                                                   >  2:322684/1‑150 (MQ=255)
                                   cGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc                                                                                                                 >  1:1031590/1‑150 (MQ=255)
                                   cGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc                                                                                                                 >  1:1388279/1‑150 (MQ=255)
                                    gCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACg                                                                                                                >  1:635444/1‑150 (MQ=255)
                                     ccTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATTGGGACgg                                                                                                               <  1:1563687/150‑1 (MQ=255)
                                        gCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTACCGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGaaa                                                                                                            <  1:1214634/150‑1 (MQ=255)
                                            cATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTac                                                                                                        >  1:243190/1‑150 (MQ=255)
                                            cATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTac                                                                                                        <  2:556163/150‑1 (MQ=255)
                                                aCTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACtt                                                                                                    >  2:1358390/1‑150 (MQ=255)
                                                  tttGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACGCttta                                                                                                  >  1:1149944/1‑149 (MQ=255)
                                                     gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAg                                                                                               <  2:78903/150‑1 (MQ=255)
                                                      taagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgt                                                                                              >  1:954969/1‑150 (MQ=255)
                                                       aagtaaAAGGAGTATGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtg                                                                                             <  2:1682162/150‑1 (MQ=255)
                                                        agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTTAATAGCATCTCCCGCCAGCAGGCGTGAGTTTTTTACTGATGCGTAGTAGATGGGGACGGAAATTACACTTTCAAgtgg                                                                                            >  2:169398/1‑149 (MQ=255)
                                                        agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtgt                                                                                            >  1:821036/1‑150 (MQ=255)
                                                        agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtgt                                                                                            >  2:1030978/1‑150 (MQ=255)
                                                         gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGtt                                                                                           >  2:500706/1‑150 (MQ=255)
                                                         gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGtt                                                                                           >  1:571764/1‑150 (MQ=255)
                                                           aaaaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTa                                                                                         >  2:1801344/1‑150 (MQ=255)
                                                              aGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAtt                                                                                      <  1:1707426/150‑1 (MQ=255)
                                                               ggtgTAGGGATTTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGTAAATTACACTTTCATGTGTTTAAttt                                                                                     >  1:898243/4‑150 (MQ=255)
                                                                    aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaac                                                                                >  2:1553498/1‑150 (MQ=255)
                                                                    aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAAGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTATTaaa                                                                                >  2:1683912/1‑149 (MQ=255)
                                                                        aaTTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaa                                                                            <  2:1294255/150‑1 (MQ=255)
                                                                          ttCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGTGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaacc                                                                          <  1:335922/150‑1 (MQ=255)
                                                                            cAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATTCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTATTaaccaaaccat                                                                        >  2:1138400/1‑149 (MQ=255)
                                                                            cAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTGTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAg                                                                        >  2:1306610/1‑150 (MQ=255)
                                                                              aGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTg                                                                      >  2:1510382/1‑150 (MQ=255)
                                                                              aGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTAGATGGGGACGGAAATTACCCTTTCAAGTGTTTAATTTTTAACCAAACCAGTg                                                                      >  2:1302423/1‑150 (MQ=255)
                                                                                 aGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTAAAGTGTTTAATTTTTAAACAAACCAGTGATg                                                                   >  2:1070375/1‑150 (MQ=255)
                                                                                  gTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATgg                                                                  <  2:1643482/150‑1 (MQ=255)
                                                                                  gTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATgg                                                                  <  2:1648958/150‑1 (MQ=255)
                                                                                   ttCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGa                                                                 <  1:910190/150‑1 (MQ=255)
                                                                                    tCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGaa                                                                >  1:426906/1‑150 (MQ=255)
                                                                                    tCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGATTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGaa                                                                >  1:1411002/1‑150 (MQ=255)
                                                                                        tAAGTTTGGTTTTCATGTCAGCTCCTGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCATAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAtt                                                            <  1:1455422/150‑1 (MQ=255)
                                                                                         aaGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAtta                                                           >  2:516717/1‑150 (MQ=255)
                                                                                         aaGTTGGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAtta                                                           <  2:355314/150‑1 (MQ=255)
                                                                                                gTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaat                                                    <  2:1635723/150‑1 (MQ=255)
                                                                                                  tttCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGt                                                  <  2:230146/150‑1 (MQ=255)
                                                                                                   ttCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTa                                                 >  2:441480/1‑150 (MQ=255)
                                                                                                         tCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAAt                                           <  1:1297527/150‑1 (MQ=255)
                                                                                                         tCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAAt                                           >  2:126540/1‑150 (MQ=255)
                                                                                                          cAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGGGATGGAACATTAATTTAATGTAATCAATg                                          <  1:1687221/150‑1 (MQ=255)
                                                                                                                  gCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGcc                                  >  1:1738960/1‑150 (MQ=255)
                                                                                                                    gtaacgtaaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCga                                >  2:908689/1‑150 (MQ=255)
                                                                                                                    gtaacgtaaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCga                                <  1:34877/150‑1 (MQ=255)
                                                                                                                      aacgtaaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCgaag                              <  2:1324635/150‑1 (MQ=255)
                                                                                                                          taaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTGTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGATCATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGt                          >  1:319527/1‑150 (MQ=255)
                                                                                                                          taaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGGATATGGGGACGTAAATTTCACTTTCTAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGt                          <  2:966364/150‑1 (MQ=255)
                                                                                                                           aaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATTTAATCAATGATTTTGCCGAAGAAGtt                         >  2:1149710/1‑150 (MQ=255)
                                                                                                                              tAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTAcc                      >  1:497328/1‑150 (MQ=255)
                                                                                                                                 aTAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTc                   >  1:1545722/1‑150 (MQ=255)
                                                                                                                                 aTAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTc                   >  1:667962/1‑150 (MQ=255)
                                                                                                                                   aaaaaCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCg                 >  2:19396/4‑150 (MQ=255)
                                                                                                                                   aaaaaCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCg                 >  2:739806/4‑150 (MQ=255)
                                                                                                                                     cAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCTAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCggg               <  2:243190/150‑1 (MQ=255)
                                                                                                                                       aCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAt             >  2:1784015/1‑150 (MQ=255)
                                                                                                                                       aCTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAt             <  2:1388279/150‑1 (MQ=255)
                                                                                                                                        cTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAtt            >  1:1200483/1‑150 (MQ=255)
                                                                                                                                        cTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAtt            >  1:431586/1‑150 (MQ=255)
                                                                                                                                         tCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACAAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAttt           >  1:1085408/1‑150 (MQ=255)
                                                                                                                                          cccGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCTGGATTTc          >  1:1111387/1‑150 (MQ=255)
                                                                                                                                            cGCCAGAAGTCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCtt        <  1:1282229/150‑1 (MQ=255)
                                                                                                                                            cGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCtt        >  2:399569/1‑150 (MQ=255)
                                                                                                                                             gCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTg       >  2:1044348/1‑150 (MQ=255)
                                                                                                                                              cgctaaGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTTCACTTTCAAGTGTTTAATTTTATAACAAACCAGTGATGGAACATTAATTTTATGTAATCAATGATTTTGACGAAGAAGTTACATTGCGGGATTTCTTGc      <  2:1487834/146‑1 (MQ=255)
                                                                                                                                              ccAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTTATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGTTTTCTTGc      <  2:1489762/150‑1 (MQ=255)
                                                                                                                                              ccAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGc      <  1:1540127/150‑1 (MQ=255)
                                                                                                                                              ccAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGc      <  2:1031590/150‑1 (MQ=255)
                                                                                                                                               cAGAAGTCGGTAGTTTTTTACTGATGCTTAGTATTTGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCa     <  2:524572/150‑1 (MQ=255)
                                                                                                                                               cAGAAGGCGGGAGTTTTTTACTGATGCTTATTATATCTGGACTGAAATTACACTTTCAAGTGTTTAATTTTAAACCAAACCAGTGATGGATCATTAATTTAAGGTAATCAATGATTTTGCCGAAGAAGGTACCTTCCGGGATTTCTTGCa     <  2:277839/150‑1 (MQ=255)
                                                                                                                                                  aaGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCAaca  >  1:1728088/1‑150 (MQ=255)
                                                                                                                                                  aaGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCAaca  <  1:1395955/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCAACA  >  NC_000913/662539‑662834

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: