Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,323,490 A→C 5.5% L194R (CTG→CGG)  glmM ← phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,323,4900AC 5.5% 279.0 / 4.7 87L194R (CTG→CGG) glmMphosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (43/39);  new base C (3/2);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCC  >  NC_000913/3323343‑3323636
                                                                                                                                                   |                                                                                                                                                  
tcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCgcagc                                                                                                                                                  <  1:900574/150‑1 (MQ=255)
tcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCgcagc                                                                                                                                                  <  2:1170973/150‑1 (MQ=255)
 cgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCgcagca                                                                                                                                                 <  2:812249/150‑1 (MQ=255)
 cgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCgcagca                                                                                                                                                 <  2:133220/150‑1 (MQ=255)
  gccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCAc                                                                                                                                                >  1:1047298/1‑150 (MQ=255)
   ccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACg                                                                                                                                               >  2:395368/1‑150 (MQ=255)
      tCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTc                                                                                                                                            <  1:496088/150‑1 (MQ=255)
          aGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGgcg                                                                                                                                        >  1:1512016/1‑150 (MQ=255)
           ggCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGgcgc                                                                                                                                       <  1:656493/150‑1 (MQ=255)
           ggCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGgcgc                                                                                                                                       <  1:1340242/150‑1 (MQ=255)
            gCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGAACGTGCTGcgcg                                                                                                                                      <  1:1560491/150‑1 (MQ=255)
             cAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGa                                                                                                                                     >  2:488903/1‑150 (MQ=255)
              aaTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAt                                                                                                                                    <  2:790683/150‑1 (MQ=255)
               aTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGAACGTGCTGCGCGaac                                                                                                                                   <  1:1316700/150‑3 (MQ=255)
               aTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtg                                                                                                                                   <  1:48570/150‑1 (MQ=255)
               aTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCG‑‑AACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtgtg                                                                                                                                 <  1:931598/150‑1 (MQ=255)
                 aCCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGAACGTGCTGCGCGaactg                                                                                                                                 >  2:1405769/1‑146 (MQ=255)
                 aCCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtgtg                                                                                                                                 <  2:1351602/150‑1 (MQ=255)
                 aCCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtgtg                                                                                                                                 <  1:414386/150‑1 (MQ=255)
                   cGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGAt                                                                                                                               >  2:155307/1‑150 (MQ=255)
                    gagaTCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATa                                                                                                                              <  2:1031903/150‑1 (MQ=255)
                        tCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTc                                                                                                                          >  1:902563/1‑150 (MQ=255)
                          cGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGc                                                                                                                        <  2:1249220/150‑1 (MQ=255)
                               tttCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGt                                                                                                                   >  2:559297/1‑150 (MQ=255)
                                  cagccagcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTg                                                                                                                <  1:128242/150‑1 (MQ=255)
                                       agcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACaa                                                                                                           <  1:1095670/150‑1 (MQ=255)
                                       agaaaacGAGCCTCGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACAACGTTTGGCTCACAAACGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGATTGCACaa                                                                                                           >  2:221673/6‑150 (MQ=255)
                                               aGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCcaccac                                                                                                   <  2:100505/150‑1 (MQ=255)
                                               aGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCcaccac                                                                                                   <  2:785663/150‑1 (MQ=255)
                                                  cTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAAt                                                                                                >  2:872334/1‑150 (MQ=255)
                                                   tGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACCCCGTTTGGCTCACACCCGATAGCGATACCGTTCGCCCCCAGTTCGCGCCGCACGTTCGGCGCGATGTGATAAGTCGCCCCGTTTGCCCAATCCACCCCAAtg                                                                                               >  1:784241/1‑149 (MQ=255)
                                                    ggAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCt                                                                                              <  2:180463/150‑1 (MQ=255)
                                                      aGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTc                                                                                            <  1:1489532/150‑1 (MQ=255)
                                                      aGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTc                                                                                            >  2:322098/1‑150 (MQ=255)
                                                            cgAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCa                                                                                      >  1:81630/1‑150 (MQ=255)
                                                            cgAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCa                                                                                      <  2:246174/150‑1 (MQ=255)
                                                            cgAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCa                                                                                      >  1:1453746/1‑150 (MQ=255)
                                                            cgAACGTCGGTAGCCACCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATACCGATAACGTTCGCACCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCa                                                                                      >  2:1183935/1‑150 (MQ=255)
                                                                    ggTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCt                                                                              >  1:1458061/1‑150 (MQ=255)
                                                                     gTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTa                                                                             <  1:687665/150‑1 (MQ=255)
                                                                     gTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTa                                                                             <  2:568975/150‑1 (MQ=255)
                                                                     gTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTa                                                                             <  2:932277/150‑1 (MQ=255)
                                                                     gTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTa                                                                             >  1:1544589/1‑150 (MQ=255)
                                                                       atcccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAg                                                                           <  2:784241/148‑1 (MQ=255)
                                                                         cccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgtt                                                                         >  2:632615/1‑150 (MQ=255)
                                                                         cccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgtt                                                                         >  1:320090/1‑150 (MQ=255)
                                                                          ccccACTTCGGCATTGATGTTTACACCGTTTGGCTAACAACCGATAGAGATAACGTTCGCCCACAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATACACCACAATCTTCAGTTCACTGAGGCTAAgtta                                                                        >  2:1110193/1‑149 (MQ=255)
                                                                           cccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgttcg                                                                       >  2:1508280/1‑150 (MQ=255)
                                                                              aCTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgttcgttc                                                                    <  1:472861/150‑1 (MQ=255)
                                                                               cTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGttcgttcg                                                                   <  1:155307/150‑1 (MQ=255)
                                                                               cTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGttcgttcg                                                                   <  2:902563/150‑1 (MQ=255)
                                                                               cTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGttcgttcg                                                                   <  1:872334/150‑1 (MQ=255)
                                                                                   ggCGTTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGaa                                                               >  1:419682/1‑150 (MQ=255)
                                                                                     cATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACg                                                             >  1:1128238/1‑150 (MQ=255)
                                                                                       ttGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTg                                                           <  2:1544589/150‑1 (MQ=255)
                                                                                           tGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCtt                                                       >  1:1207544/1‑150 (MQ=255)
                                                                                             tttACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGTGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTg                                                     <  2:81630/150‑1 (MQ=255)
                                                                                              ttACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGc                                                    >  1:767771/1‑150 (MQ=255)
                                                                                              ttACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCAACACAATCTTCAGTTAAATGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGc                                                    >  2:636833/1‑150 (MQ=255)
                                                                                               tACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCa                                                   >  2:463173/1‑150 (MQ=255)
                                                                                                 cacCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGACCGTGGCTTTGCaaa                                                 >  1:1049721/1‑150 (MQ=255)
                                                                                                 cacCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGAAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCaaaa                                                >  2:1445204/1‑150 (MQ=255)
                                                                                                   ccGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAAc                                               >  2:1381047/1‑150 (MQ=255)
                                                                                                   ccGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAAc                                               <  1:559297/150‑1 (MQ=255)
                                                                                                     gTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTc                                             >  2:722481/1‑150 (MQ=255)
                                                                                                      tttGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCg                                            <  2:320090/150‑1 (MQ=255)
                                                                                                      tttGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCg                                            <  2:324006/150‑1 (MQ=255)
                                                                                                        tGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGat                                          >  2:246285/1‑150 (MQ=255)
                                                                                                        tGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGat                                          >  2:261396/1‑150 (MQ=255)
                                                                                                            tCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAg                                      <  2:1363893/150‑1 (MQ=255)
                                                                                                            tCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAg                                      >  1:193092/1‑150 (MQ=255)
                                                                                                             cacaACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGc                                     >  2:262116/1‑150 (MQ=255)
                                                                                                                   cGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCcg                               <  1:261396/150‑1 (MQ=255)
                                                                                                                    gATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCcgc                              >  2:934693/1‑150 (MQ=255)
                                                                                                                       aGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGc                           <  1:463173/150‑1 (MQ=255)
                                                                                                                         cGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCAt                         >  1:816737/1‑150 (MQ=255)
                                                                                                                         cGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCAt                         >  2:790293/1‑150 (MQ=255)
                                                                                                                          gATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATc                        >  2:372617/1‑150 (MQ=255)
                                                                                                                          gATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATc                        <  1:1508280/150‑1 (MQ=255)
                                                                                                                               cGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGa                   >  1:903215/1‑150 (MQ=255)
                                                                                                                               cGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGa                   >  1:1462813/1‑150 (MQ=255)
                                                                                                                               cGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGa                   >  1:490089/1‑150 (MQ=255)
                                                                                                                                gTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGAt                  <  1:722481/150‑1 (MQ=255)
                                                                                                                                 ttCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGa                                           <  2:837496/124‑1 (MQ=255)
                                                                                                                                 ttCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGa                                           >  1:837496/1‑124 (MQ=255)
                                                                                                                                     cccccAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggc             <  1:632615/150‑1 (MQ=255)
                                                                                                                                      ccccAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCTCCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggct            >  1:55421/1‑150 (MQ=255)
                                                                                                                                       cccAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggctg           >  1:186967/1‑150 (MQ=255)
                                                                                                                                           catCGCGCCGAACGTGCTGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCtt       >  1:613181/3‑150 (MQ=255)
                                                                                                                                               gcgcAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTAcc   <  1:1381047/150‑1 (MQ=255)
                                                                                                                                               gcgcAGCACGTTCGGCGCGATGTGAAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTAccc  <  1:1445204/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCC  >  NC_000913/3323343‑3323636

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