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

AAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGG  >  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)
             gcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGc                                                                                                                                       >  1:28375/1‑150 (MQ=255)
               aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTg                                                                                                                                     <  1:1337499/150‑1 (MQ=255)
               aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtc                                                                                                                                     <  2:537161/150‑2 (MQ=255)
                aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTgg                                                                                                                                             >  2:1502546/1‑141 (MQ=255)
                aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTgg                                                                                                                                             <  1:1502546/141‑1 (MQ=255)
                aTCAGGTTAAACACTTCGATGATGATCGAGAAGGCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTgg                                                                                                                                    <  2:961133/150‑1 (MQ=255)
                 tCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGTCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtctg                                                                                                                                   <  2:1706449/150‑4 (MQ=255)
                  cTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGtt                                                                                                                                  >  1:1619346/1‑150 (MQ=255)
                   tgtttAAACACTTCGATGAAGATCGAGAAGCCAATCGCGGCATACAGGTAACCATTCGGAATGTGGAAACCGAAACCTTCTGCACCCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCCGGGGGTCCAGGTTa                                                                                                                                 >  2:736591/4‑150 (MQ=255)
                   tGGTTAAACACTTCGATGAGGAGCGAGAAGCAAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACGGAAAACTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtctgtc                                                                                                                                 <  2:1171732/150‑6 (MQ=255)
                     gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                               <  2:1307106/150‑1 (MQ=255)
                     gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                               >  1:754693/1‑150 (MQ=255)
                     gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                               <  1:706937/150‑1 (MQ=255)
                        aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGaa                                                                                                                            <  2:1613266/150‑1 (MQ=255)
                        aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGaa                                                                                                                            >  2:621165/1‑150 (MQ=255)
                                 gatgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGATTGTGGAAACCGAAACCTTCTGCCACCAGACTCCGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCCCCCGTGGGGTGCTGGTTAACGAATCGCGTctg                                                                                                                   >  2:1552191/1‑148 (MQ=255)
                                  atgatgatCGAGAAGCCAATCGCTGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGc                                                                                                                  <  1:859365/150‑1 (MQ=255)
                                      tgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGtt                                                                                                              >  1:742791/1‑150 (MQ=255)
                                      tgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAAAAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTTACGAATCGCGTCAGCGGtt                                                                                                              >  2:550568/1‑150 (MQ=255)
                                       tatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                             <  2:408779/149‑1 (MQ=255)
                                       gatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                             >  2:1268894/1‑150 (MQ=255)
                                          cGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTggt                                                                                                          >  2:1430212/1‑149 (MQ=255)
                                            agaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATg                                                                                                        <  2:28375/150‑1 (MQ=255)
                                            agaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATg                                                                                                        >  2:21376/1‑150 (MQ=255)
                                            agaAGCCAATCGCGGACTACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCAGCAGACTCAGACGAATCATTAACTGGAAGCTCAGACAGAGCACCACCACCGTGTGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATg                                                                                                        <  1:163410/150‑1 (MQ=255)
                                              aaGCCAATCGCGGCCTACAGGGAACCTAGCGGAATGTGGTAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGcc                                                                                                      <  2:250564/150‑1 (MQ=39)
                                               aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     <  1:363547/150‑1 (MQ=255)
                                               aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     <  1:1268894/150‑1 (MQ=255)
                                               aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     >  2:1102874/1‑150 (MQ=255)
                                               aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     >  1:1238917/1‑150 (MQ=255)
                                                    aTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAAc                                                                                                <  1:279484/150‑1 (MQ=255)
                                                       gcgGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATa                                                                                             <  1:343322/150‑1 (MQ=255)
                                                       gcgGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATa                                                                                             <  1:91440/150‑1 (MQ=255)
                                                         ggCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCAACACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCCGTATGGATGCCAGCAACATAAc                                                                                           >  2:1626270/1‑150 (MQ=255)
                                                          gCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                          <  2:754693/150‑1 (MQ=255)
                                                          gCAAACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGTCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                          <  1:858308/150‑1 (MQ=255)
                                                           cATACAGCTAACCTATCGGCATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCACCGGATTGGATGCCAGCAACATAACCg                                                                                         <  1:1571841/150‑1 (MQ=255)
                                                                 ggTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCg                                                                                   >  2:1346314/1‑150 (MQ=255)
                                                                     aCCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAAt                                                                               <  1:1348155/150‑1 (MQ=255)
                                                                      ccTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGATGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                              <  1:1528323/150‑1 (MQ=255)
                                                                      ccTTTCGGAATGAGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                              <  2:857290/150‑1 (MQ=255)
                                                                        tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCAc                                                                            >  1:1270121/1‑150 (MQ=255)
                                                                        tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCAc                                                                            <  1:1226017/150‑1 (MQ=255)
                                                                        tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATAATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCAAGCAACATAACCTCCATCGCAATCAc                                                                            >  2:1125108/1‑150 (MQ=255)
                                                                         ttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCAca                                                                           >  2:561552/1‑149 (MQ=255)
                                                                          tCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATAACTa                                                                          >  1:1510877/1‑150 (MQ=255)
                                                                           cGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAc                                                                         >  1:1481460/1‑150 (MQ=255)
                                                                            ggAATGTGGAAACCGAAACCTTCTGCCACCAGCCTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGAAGCCAGCAACATAACCGCCATCGCAATCACTAcc                                                                        >  2:1163135/1‑150 (MQ=255)
                                                                            ggAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAcc                                                                        >  2:1626825/1‑150 (MQ=255)
                                                                            ggAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAcc                                                                        >  1:1676606/1‑150 (MQ=255)
                                                                                tgtgGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGCAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAccgccg                                                                    <  1:1520193/150‑1 (MQ=255)
                                                                                   ggAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACGACCGCCGCca                                                                 >  2:1045194/1‑150 (MQ=255)
                                                                                        ccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCatcatc                                                            >  2:995498/1‑150 (MQ=255)
                                                                                          gaaaccTTCTGCCACCAGCCTCAGACAAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTCCTGGTTAACGAATCTCGTCACCGGTTTGGATGCCAGCAACATCACCGCCATCGCAATCACTACCGCCCCcatccacac                                                          >  1:58911/1‑145 (MQ=255)
                                                                                          gaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCAc                                                          >  1:1341473/1‑150 (MQ=255)
                                                                                                  cTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGAt                                                  <  1:1620165/150‑1 (MQ=255)
                                                                                                   tGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATg                                                 <  2:878706/150‑1 (MQ=255)
                                                                                                    gCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCCGCAACATAACCGCCATCGCAATCACTACCGCCGCCAGCCTCCCCGGCAGATgg                                                >  2:1245532/1‑150 (MQ=255)
                                                                                                     ccaccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACAGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGt                                               >  1:1239295/1‑150 (MQ=255)
                                                                                                      caccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                              >  2:576400/1‑150 (MQ=255)
                                                                                                      caccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                              <  2:161175/150‑1 (MQ=255)
                                                                                                             cTCAGACCAATCATTAACAGGAAGCTCAGACTGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCCCCGGCAGATGGTTAACCCTc                                       >  1:1266471/1‑150 (MQ=255)
                                                                                                             cTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGACTCGCGTCCGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCCCTACCGCCGCCCTCATCCCCGGCCGATGGTTGACCATc                                       >  2:101048/1‑150 (MQ=255)
                                                                                                             cTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTACCCATc                                       >  2:1086884/1‑150 (MQ=255)
                                                                                                             cTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATc                                       >  2:1367825/1‑150 (MQ=255)
                                                                                                             cTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTAAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATACCCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATc                                       >  1:29674/1‑150 (MQ=255)
                                                                                                              tCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATcc                                      <  1:621165/150‑1 (MQ=255)
                                                                                                               cAGAGCAATCATTAACAGGAAGCTCAGACAGAGCACCACCCCCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGGTAACCATccc                                     >  2:1241666/1‑150 (MQ=255)
                                                                                                               cAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATccc                                     >  1:509166/1‑150 (MQ=255)
                                                                                                                 taCCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTAGGATCACAGCAACATAACCGCCATCGCAATCACTACAGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTa                                   >  1:1116170/2‑150 (MQ=255)
                                                                                                                 gACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCCCTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTa                                   >  2:945064/1‑150 (MQ=255)
                                                                                                                 gACCAATCATTAACAGGAAGCTCAGACAGAGCACAACCAACGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCACCCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATACCCGCCAGATGGTTAACCATCCCTa                                   >  1:1116174/1‑150 (MQ=255)
                                                                                                                   ccAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACt                                 <  2:1238917/150‑1 (MQ=255)
                                                                                                                     aaTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGATTGGATGCCAGCAACAACACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGAAGGGTCACCATCCCTTCTGc                               >  2:1726321/1‑150 (MQ=255)
                                                                                                                      aTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCa                              <  2:30739/150‑1 (MQ=255)
                                                                                                                        cATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGt                            >  2:600786/1‑150 (MQ=255)
                                                                                                                         aTTAACAGGAAGCTCAGACATAGCACCAACACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTa                           <  1:1296957/150‑1 (MQ=255)
                                                                                                                         aTTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTa                           >  2:167418/1‑150 (MQ=255)
                                                                                                                         aTTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTa                           >  2:574125/1‑150 (MQ=255)
                                                                                                                            aaCATGAAGCTCAGAAAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATAGCTTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAACGCCGCCATCATCACCGGCAGATTGTTAACCATCCCTACTGCAGTAAtt                        <  1:87312/150‑1 (MQ=255)
                                                                                                                            aaCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  2:1369751/1‑150 (MQ=255)
                                                                                                                            aaCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  2:1163380/1‑150 (MQ=255)
                                                                                                                            aaCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  1:1370167/1‑150 (MQ=255)
                                                                                                                             aCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTa                       >  1:930941/1‑150 (MQ=255)
                                                                                                                                ggAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACAATCCCTACTGCAGTAAATACCg                    >  2:108327/1‑150 (MQ=255)
                                                                                                                                       cagacagaGCACCACCCCCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCaa             >  2:1029760/1‑150 (MQ=255)
                                                                                                                                       cagacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTTCCGCCGCCATCATCACCGGCCTCTGGTGAACCCTCCCGACGGCAGTCATTACCGCATCCaa             >  2:249480/1‑150 (MQ=255)
                                                                                                                                       cagacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCaa             >  2:508984/1‑150 (MQ=255)
                                                                                                                                        agacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc            <  1:1015887/150‑1 (MQ=255)
                                                                                                                                        agacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc            <  1:1169731/150‑1 (MQ=255)
                                                                                                                                         gacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACg           >  2:1637389/1‑150 (MQ=255)
                                                                                                                                           cagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTACCCATCCCTACTGCAGTAATTACCGCATCCAACgag         >  1:1049893/1‑150 (MQ=255)
                                                                                                                                           cagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCATCgag         <  2:522649/150‑1 (MQ=255)
                                                                                                                                            ccagCACCCCACGGTGGTGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACgaga        >  1:644097/3‑150 (MQ=255)
                                                                                                                                              agCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGagaag      >  1:1470120/1‑150 (MQ=255)
                                                                                                                                              agCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGAGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGagaag      <  1:235870/150‑1 (MQ=255)
                                                                                                                                                  ccaccaccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACgg  >  2:1662602/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGG  >  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: