Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,034 C→G 6.2% T184T (ACG→ACC)  yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,901,0340CG 6.2% 300.9 / 9.2 97T184T (ACG→ACC) yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base C (54/36);  new base G (2/4);  total (57/40)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

CGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGG  >  NC_000913/1900888‑1901176
                                                                                                                                                  |                                                                                                                                               
cGACGCGAAATCTGGTTAAACACTTCGATGATGATCTAGAAGCAAATCGCGGGATACAGGTAAGCTTTAGGAATGTGGAAACCGAAACCTTCTGGCCCCAGACTCAGACCAATCAATAACAGGTAGCTCAGACAGAGGACCACCACCGTg                                                                                                                                              <  1:84699/150‑1 (MQ=255)
      gcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGc                                                                                                                                       >  1:28375/1‑150 (MQ=255)
        aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTg                                                                                                                                     <  1:1337499/150‑1 (MQ=255)
        aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtc                                                                                                                                     <  2:537161/150‑2 (MQ=255)
         aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTgg                                                                                                                                             >  2:1502546/1‑141 (MQ=255)
         aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTgg                                                                                                                                             <  1:1502546/141‑1 (MQ=255)
         aTCAGGTTAAACACTTCGATGATGATCGAGAAGGCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTgg                                                                                                                                    <  2:961133/150‑1 (MQ=255)
          tCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGTCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtctg                                                                                                                                   <  2:1706449/150‑4 (MQ=255)
           cTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGtt                                                                                                                                  >  1:1619346/1‑150 (MQ=255)
            tgtttAAACACTTCGATGAAGATCGAGAAGCCAATCGCGGCATACAGGTAACCATTCGGAATGTGGAAACCGAAACCTTCTGCACCCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCCGGG‑GGTCCAGGTTa                                                                                                                                 >  2:736591/4‑150 (MQ=255)
            tGGTTAAACACTTCGATGAGGAGCGAGAAGCAAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACGGAAAACTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtctgtc                                                                                                                                 <  2:1171732/150‑6 (MQ=255)
              gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAAc                                                                                                                               <  2:1307106/150‑1 (MQ=255)
              gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAAc                                                                                                                               <  1:706937/150‑1 (MQ=255)
              gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAAc                                                                                                                               >  1:754693/1‑150 (MQ=255)
                 aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGaa                                                                                                                            >  2:621165/1‑150 (MQ=255)
                 aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGaa                                                                                                                            <  2:1613266/150‑1 (MQ=255)
                          gatgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGATTGTGGAAACCGAAACCTTCTGCCACCAGACTCCGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCCCCCGTGG‑GGTGCTGGTTAACGAATCGCGTctg                                                                                                                   >  2:1552191/1‑148 (MQ=255)
                           atgatgatCGAGAAGCCAATCGCTGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGc                                                                                                                  <  1:859365/150‑1 (MQ=255)
                               tgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGtt                                                                                                              >  1:742791/1‑150 (MQ=255)
                               tgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAAAAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTTACGAATCGCGTCAGCGGtt                                                                                                              >  2:550568/1‑150 (MQ=255)
                                tatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                             <  2:408779/149‑1 (MQ=255)
                                gatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                             >  2:1268894/1‑150 (MQ=255)
                                   cGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTggt                                                                                                          >  2:1430212/1‑149 (MQ=255)
                                     agaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATg                                                                                                        <  2:28375/150‑1 (MQ=255)
                                     agaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATg                                                                                                        >  2:21376/1‑150 (MQ=255)
                                     agaAGCCAATCGCGGACTACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCAGCAGACTCAGACGAATCATTAACTGGAAGCTCAGACAGAGCACCACCACCGTGT‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATg                                                                                                        <  1:163410/150‑1 (MQ=255)
                                       aaGCCAATCGCGGCCTACAGGGAACCTAGCGGAATGTGGTAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGcc                                                                                                      <  2:250564/150‑1 (MQ=39)
                                        aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     <  1:1268894/150‑1 (MQ=255)
                                        aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     >  1:1238917/1‑150 (MQ=255)
                                        aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     >  2:1102874/1‑150 (MQ=255)
                                        aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     <  1:363547/150‑1 (MQ=255)
                                             aTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAAc                                                                                                <  1:279484/150‑1 (MQ=255)
                                                gcgGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATa                                                                                             <  1:91440/150‑1 (MQ=255)
                                                gcgGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATa                                                                                             <  1:343322/150‑1 (MQ=255)
                                                  ggCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCAACACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCCGTATGGATGCCAGCAACATAAc                                                                                           >  2:1626270/1‑150 (MQ=255)
                                                   gCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                          <  2:754693/150‑1 (MQ=255)
                                                   gCAAACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGTCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                          <  1:858308/150‑1 (MQ=255)
                                                    cATACAGCTAACCTATCGGCATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCACCGGATTGGATGCCAGCAACATAACCg                                                                                         <  1:1571841/150‑1 (MQ=255)
                                                          ggTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCg                                                                                   >  2:1346314/1‑150 (MQ=255)
                                                              aCCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAAt                                                                               <  1:1348155/150‑1 (MQ=255)
                                                               ccTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GATGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                              <  1:1528323/150‑1 (MQ=255)
                                                               ccTTTCGGAATGAGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                              <  2:857290/150‑1 (MQ=255)
                                                                 tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCAc                                                                            <  1:1226017/150‑1 (MQ=255)
                                                                 tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCAc                                                                            >  1:1270121/1‑150 (MQ=255)
                                                                 tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATAATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCAAGCAACATAACCTCCATCGCAATCAc                                                                            >  2:1125108/1‑150 (MQ=255)
                                                                  ttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCAca                                                                           >  2:561552/1‑149 (MQ=255)
                                                                   tCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATAACTa                                                                          >  1:1510877/1‑150 (MQ=255)
                                                                    cGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAc                                                                         >  1:1481460/1‑150 (MQ=255)
                                                                     ggAATGTGGAAACCGAAACCTTCTGCCACCAGCCTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGAAGCCAGCAACATAACCGCCATCGCAATCACTAcc                                                                        >  2:1163135/1‑150 (MQ=255)
                                                                     ggAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAcc                                                                        >  2:1626825/1‑150 (MQ=255)
                                                                     ggAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAcc                                                                        >  1:1676606/1‑150 (MQ=255)
                                                                         tgtgGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGCAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAccgccg                                                                    <  1:1520193/150‑1 (MQ=255)
                                                                            ggAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACGACCGCCGCca                                                                 >  2:1045194/1‑150 (MQ=255)
                                                                                 ccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCatcatc                                                            >  2:995498/1‑150 (MQ=255)
                                                                                   gaaaccTTCTGCCACCAGCCTCAGACAAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTCCTGGTTAACGAATCTCGTCACCGGTTTGGATGCCAGCAACATCACCGCCATCGCAATCACTACCGCCCCcatccacac                                                          >  1:58911/1‑145 (MQ=255)
                                                                                   gaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCAc                                                          >  1:1341473/1‑150 (MQ=255)
                                                                                           cTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGAt                                                  <  1:1620165/150‑1 (MQ=255)
                                                                                            tGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATg                                                 <  2:878706/150‑1 (MQ=255)
                                                                                             gCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCCGCAACATAACCGCCATCGCAATCACTACCGCCGCCAGCCTCCCCGGCAGATgg                                                >  2:1245532/1‑150 (MQ=255)
                                                                                              ccaccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACAGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGt                                               >  1:1239295/1‑150 (MQ=255)
                                                                                               caccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                              >  2:576400/1‑150 (MQ=255)
                                                                                               caccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                              <  2:161175/150‑1 (MQ=255)
                                                                                                      cTCAGACCAATCATTAACAGGAAGCTCAGACTGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCCCCGGCAGATGGTTAACCCTc                                       >  1:1266471/1‑150 (MQ=255)
                                                                                                      cTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGACTCGCGTCCGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCCCTACCGCCGCCCTCATCCCCGGCCGATGGTTGACCATc                                       >  2:101048/1‑150 (MQ=255)
                                                                                                      cTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTACCCATc                                       >  2:1086884/1‑150 (MQ=255)
                                                                                                      cTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATc                                       >  2:1367825/1‑150 (MQ=255)
                                                                                                      cTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTAAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATACCCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATc                                       >  1:29674/1‑150 (MQ=255)
                                                                                                       tCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATcc                                      <  1:621165/150‑1 (MQ=255)
                                                                                                        cAGAGCAATCATTAACAGGAAGCTCAGACAGAGCACCACCCCCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGGTAACCATccc                                     >  2:1241666/1‑150 (MQ=255)
                                                                                                        cAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATccc                                     >  1:509166/1‑150 (MQ=255)
                                                                                                          taCCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTAGGATCACAGCAACATAACCGCCATCGCAATCACTACAGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTa                                   >  1:1116170/2‑150 (MQ=255)
                                                                                                          gACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCCCTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTa                                   >  2:945064/1‑150 (MQ=255)
                                                                                                          gACCAATCATTAACAGGAAGCTCAGACAGAGCACAACCAACGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCACCCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATACCCGCCAGATGGTTAACCATCCCTa                                   >  1:1116174/1‑150 (MQ=255)
                                                                                                            ccAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACt                                 <  2:1238917/150‑1 (MQ=255)
                                                                                                              aaTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGATTGGATGCCAGCAACAACACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGAAGGGTCACCATCCCTTCTGc                               >  2:1726321/1‑150 (MQ=255)
                                                                                                               aTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCa                              <  2:30739/150‑1 (MQ=255)
                                                                                                                 cATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGt                            >  2:600786/1‑150 (MQ=255)
                                                                                                                  aTTAACAGGAAGCTCAGACATAGCACCAACACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTa                           <  1:1296957/150‑1 (MQ=255)
                                                                                                                  aTTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTa                           >  2:167418/1‑150 (MQ=255)
                                                                                                                  aTTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTa                           >  2:574125/1‑150 (MQ=255)
                                                                                                                     aaCATGAAGCTCAGAAAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATAGCTTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAACGCCGCCATCATCACCGGCAGATTGTTAACCATCCCTACTGCAGTAAtt                        <  1:87312/150‑1 (MQ=255)
                                                                                                                     aaCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  1:1370167/1‑150 (MQ=255)
                                                                                                                     aaCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  2:1163380/1‑150 (MQ=255)
                                                                                                                     aaCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  2:1369751/1‑150 (MQ=255)
                                                                                                                      aCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTa                       >  1:930941/1‑150 (MQ=255)
                                                                                                                         ggAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACAATCCCTACTGCAGTAAATACCg                    >  2:108327/1‑150 (MQ=255)
                                                                                                                                cagacagaGCACCACCCCCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCaa             >  2:1029760/1‑150 (MQ=255)
                                                                                                                                cagacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTTCCGCCGCCATCATCACCGGCCTCTGGTGAACCCTCCCGACGGCAGTCATTACCGCATCCaa             >  2:249480/1‑150 (MQ=255)
                                                                                                                                cagacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCaa             >  2:508984/1‑150 (MQ=255)
                                                                                                                                 agacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc            <  1:1015887/150‑1 (MQ=255)
                                                                                                                                 agacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc            <  1:1169731/150‑1 (MQ=255)
                                                                                                                                  gacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACg           >  2:1637389/1‑150 (MQ=255)
                                                                                                                                    cagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTACCCATCCCTACTGCAGTAATTACCGCATCCAACgag         >  1:1049893/1‑150 (MQ=255)
                                                                                                                                    cagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCATCgag         <  2:522649/150‑1 (MQ=255)
                                                                                                                                     ccagCACC‑CCACGGTGGTGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACgaga        >  1:644097/3‑150 (MQ=255)
                                                                                                                                       agCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGagaag      >  1:1470120/1‑150 (MQ=255)
                                                                                                                                       agCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGAGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGagaag      <  1:235870/150‑1 (MQ=255)
                                                                                                                                           ccaccaccGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACgg  >  2:1662602/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGG  >  NC_000913/1900888‑1901176

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: ‑