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→ACCyoaEputative 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

CGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGG  >  NC_000913/1900888‑1901176
                                                                                                                                                  |                                                                                                                                               
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:706937/150‑1 (MQ=255)
              gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                               >  1:754693/1‑150 (MQ=255)
                 aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGaa                                                                                                                            >  2:621165/1‑150 (MQ=255)
                 aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGaa                                                                                                                            <  2:1613266/150‑1 (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:1268894/150‑1 (MQ=255)
                                        aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     >  1:1238917/1‑150 (MQ=255)
                                        aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     >  2:1102874/1‑150 (MQ=255)
                                        aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                     <  1:363547/150‑1 (MQ=255)
                                             aTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAAc                                                                                                <  1:279484/150‑1 (MQ=255)
                                                gcgGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATa                                                                                             <  1:91440/150‑1 (MQ=255)
                                                gcgGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATa                                                                                             <  1:343322/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:1226017/150‑1 (MQ=255)
                                                                 tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCAc                                                                            >  1:1270121/1‑150 (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                        >  1:1370167/1‑150 (MQ=255)
                                                                                                                     aaCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  2:1163380/1‑150 (MQ=255)
                                                                                                                     aaCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                        >  2:1369751/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)
                                                                                                                                                  |                                                                                                                                               
CGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGG  >  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: