Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,029 Δ1 bp 6.2% coding (557/1557 nt) 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,0290A. 6.2% 281.2 / 12.4 97coding (557/1557 nt)yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base A (54/36);  new base . (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 = 9.15e-01

AAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGG  >  NC_000913/1900881‑1901176
                                                                                                                                                    |                                                                                                                                                    
aaaGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGcaccac                                                                                                                                                     >  2:1226017/1‑150 (MQ=255)
   gTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGcaccaccac                                                                                                                                                  >  2:1620165/1‑150 (MQ=255)
       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:754693/1‑150 (MQ=255)
                     gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAAc                                                                                                                               <  1:706937/150‑1 (MQ=255)
                        aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGaa                                                                                                                            <  2:1613266/150‑1 (MQ=255)
                        aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGaa                                                                                                                            >  2:621165/1‑150 (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:363547/150‑1 (MQ=255)
                                               aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     <  1:1268894/150‑1 (MQ=255)
                                               aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     >  2:1102874/1‑150 (MQ=255)
                                               aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     >  1:1238917/1‑150 (MQ=255)
                                                    aTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAAc                                                                                                <  1:279484/150‑1 (MQ=255)
                                                       gcgGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATa                                                                                             <  1:343322/150‑1 (MQ=255)
                                                       gcgGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATa                                                                                             <  1:91440/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:1270121/1‑150 (MQ=255)
                                                                        tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCAc                                                                            <  1:1226017/150‑1 (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                        >  2:1369751/1‑150 (MQ=255)
                                                                                                                            aaCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  2:1163380/1‑150 (MQ=255)
                                                                                                                            aaCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  1:1370167/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)
                                                                                                                                                    |                                                                                                                                                    
AAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGG  >  NC_000913/1900881‑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: ‑