Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,323,492 C→A 5.2% V193V (GTG→GTT)  glmM ← phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,323,4920CA 5.2% 378.8 / 7.3 123V193V (GTG→GTT) glmMphosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base C (61/55);  new base A (3/4);  total (64/59)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

TCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTC  >  NC_000913/3323343‑3323641
                                                                                                                                                     |                                                                                                                                                     
 cgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCgcagca                                                                                                                                                      >  2:1636444/1‑150 (MQ=255)
  gccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCAc                                                                                                                                                     >  2:621310/1‑150 (MQ=255)
        gAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTGCGCGCAGCACGTTCgg                                                                                                                                               <  1:849663/150‑1 (MQ=255)
           ggCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGgcgc                                                                                                                                            >  2:822400/1‑150 (MQ=255)
           ggCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCCCAGCACGATCTgcgc                                                                                                                                            >  2:352437/1‑150 (MQ=255)
             cAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGa                                                                                                                                          >  1:1953074/1‑150 (MQ=255)
              aaTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGCTGGTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGAACGTGCTGCGCGaa                                                                                                                                         <  1:47910/150‑2 (MQ=255)
              aaTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGATTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACTTTCGGCGCGAt                                                                                                                                         >  1:1493289/1‑150 (MQ=255)
              aaTACCGAGATCCACTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAt                                                                                                                                         >  1:65954/1‑150 (MQ=255)
               aTACCGCGATCCGCTTATTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGAACGTGCTGCGCGaac                                                                                                                                        <  2:500921/150‑3 (MQ=255)
               aTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGAACGTGCTGCGCGaac                                                                                                                                        <  1:1686202/150‑3 (MQ=255)
               aTACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtg                                                                                                                                        >  1:1061666/1‑150 (MQ=255)
                tACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtgt                                                                                                                                       <  1:1320829/150‑1 (MQ=255)
                tACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGAtgt                                                                                                                                       <  1:193934/150‑1 (MQ=255)
                   cGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGAt                                                                                                                                    >  2:516796/1‑150 (MQ=255)
                    gagaTCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATa                                                                                                                                   <  2:1600757/150‑1 (MQ=255)
                    gagaTCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGATCGCCCCTAGTTCGCGCCGAACGTGCTGCGCGaactgggg                                                                                                                                   <  1:421815/150‑8 (MQ=255)
                      gaTCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAg                                                                                                                                 <  2:657619/150‑1 (MQ=255)
                      gaTCCGCTGTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCCCTTCGGCATTGATGTTTACACCGTTTGGCTCTCAACCGATAGAGATAACGTTCGCCCCCAGTTCGCGCAGACCGTTCGGCGCGATGTGATAAg                                                                                                                                 >  2:1075119/1‑150 (MQ=255)
                         ccGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCATTGCGGGCAGCCCGTTCGGCGCGATGTGATAAGTCg                                                                                                                              >  2:988833/1‑150 (MQ=255)
                           gCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCa                                                                                                                            <  2:41993/150‑1 (MQ=255)
                            cTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCAc                                                                                                                           <  2:586261/150‑1 (MQ=255)
                             tttttCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCAcc                                                                                                                          >  1:1804419/1‑150 (MQ=255)
                               tttCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGt                                                                                                                        >  1:142829/1‑150 (MQ=255)
                                   agccagcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGc                                                                                                                    <  1:1416643/150‑1 (MQ=255)
                                    gccagcTCACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATGGATGTTTACACCGTATGGCTCACAACCGATAGCGATAACGTGCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATATGTGGCACCGTTTGca                                                                                                                   <  2:1432437/150‑1 (MQ=255)
                                    gccagcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGca                                                                                                                   <  2:1256205/150‑1 (MQ=255)
                                    gccagcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAACACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGca                                                                                                                   >  1:978629/1‑150 (MQ=255)
                                      cagcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGcaca                                                                                                                 >  2:1960438/1‑150 (MQ=255)
                                        gcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAAt                                                                                                               <  1:249674/150‑1 (MQ=255)
                                            acGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATccac                                                                                                           <  1:937937/150‑1 (MQ=255)
                                             cGAGCCTGGAGCGCGCGAACGTCGGTAGACCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCGGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATccacc                                                                                                          >  1:103978/1‑150 (MQ=255)
                                                gCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGCTAAGTCGCACCGTTTGCACAATCCACcaca                                                                                                       >  2:43442/1‑150 (MQ=255)
                                                gCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCAATGCTGTTTACACCGTTTGGGTCACAACCGATAGCGATAACGTACGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTAGCAGAATCCACcaca                                                                                                       >  2:287998/1‑150 (MQ=255)
                                                  cTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAAt                                                                                                     >  1:487421/1‑150 (MQ=255)
                                                    ggAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCt                                                                                                   <  1:710392/150‑1 (MQ=255)
                                                      aGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTc                                                                                                 <  2:93909/150‑1 (MQ=255)
                                                         gcgcgAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGt                                                                                              <  2:1419231/150‑1 (MQ=255)
                                                         gcgcgAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGt                                                                                              >  1:965309/1‑150 (MQ=255)
                                                          cgcgAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGtt                                                                                             >  2:2014637/1‑150 (MQ=255)
                                                          cgcgAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGtt                                                                                             >  1:911949/1‑150 (MQ=255)
                                                              aaCGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACt                                                                                         <  1:822401/150‑1 (MQ=255)
                                                                cGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGa                                                                                       <  2:1611074/150‑1 (MQ=255)
                                                                 gTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCG‑CCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACCATCCACCACAATCTTCAGTTCACTGAgg                                                                                     >  2:88357/1‑150 (MQ=255)
                                                                   cGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGc                                                                                    <  2:1661990/150‑1 (MQ=255)
                                                                        gCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgt                                                                               >  1:1448089/1‑150 (MQ=255)
                                                                        gCCCCCACTTCGGCATTGATGTTTACACAGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCCCAATCCACCACAATCTTCAGTTCACTGAGGCTAAgt                                                                               >  1:845466/1‑150 (MQ=255)
                                                                         cccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCCCCACAATCTTCAGTTCACTGAGGCTATgtt                                                                              >  2:1510802/1‑150 (MQ=255)
                                                                         cccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgtt                                                                              >  1:1539444/1‑150 (MQ=255)
                                                                         cccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgtt                                                                              >  1:575010/1‑150 (MQ=255)
                                                                          ccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgttc                                                                             >  2:1084284/1‑150 (MQ=255)
                                                                            ccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGAACGTTCGGCGAGATGTGATAAGTCGCACCGATTGCACAATCCACCACAAACTTCAGTTCACTGAGCCTAAgttcgt                                                                           >  2:1060450/1‑150 (MQ=255)
                                                                               cTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGttcgttcg                                                                        <  1:1315504/150‑1 (MQ=255)
                                                                               cTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGAACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGttcgttcg                                                                        >  2:720065/1‑150 (MQ=255)
                                                                                    gCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAAc                                                                   >  1:118010/1‑150 (MQ=255)
                                                                                     cATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACg                                                                  <  1:1132763/150‑1 (MQ=255)
                                                                                       ttGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTg                                                                >  1:1769617/1‑150 (MQ=255)
                                                                                        tGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCCGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTgg                                                               >  2:115066/1‑150 (MQ=255)
                                                                                         gATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGCACGTGGc                                                              >  2:269904/1‑150 (MQ=255)
                                                                                            gTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCTTTCGGGAACGTGTCttt                                                           >  2:1261975/1‑150 (MQ=255)
                                                                                            gTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCttt                                                           >  2:1384028/1‑150 (MQ=255)
                                                                                             tttACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTg                                                          <  1:555468/150‑1 (MQ=255)
                                                                                              ttACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGc                                                         <  1:814254/150‑1 (MQ=255)
                                                                                              ttACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGc                                                         <  1:717572/150‑1 (MQ=255)
                                                                                               taaacCGGTTGGCACACAACCGAAAGCGTAAACCTTCCCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGGCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCa                                                        <  1:988834/147‑1 (MQ=255)
                                                                                                     gTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTc                                                  <  2:1702237/150‑1 (MQ=255)
                                                                                                      tttGGCTCACACCCGCTAGAGATAACGTTCGCCCCCAGTTCGCGCAGCACGTGCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGGTCGGGAACGTGGCTTTGCAAAACTCg                                                 <  2:212916/150‑1 (MQ=255)
                                                                                                       ttGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGa                                                <  2:1804418/150‑1 (MQ=255)
                                                                                                          gCTCACACCCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCCCGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCCCTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAACCTCGatat                                             >  1:921039/1‑150 (MQ=255)
                                                                                                          gCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGatat                                             >  2:808479/1‑150 (MQ=255)
                                                                                                             cacaACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGc                                          <  1:1835930/150‑1 (MQ=255)
                                                                                                             cacaACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGc                                          >  2:759876/1‑150 (MQ=255)
                                                                                                             cacaACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGc                                          <  2:1861411/150‑1 (MQ=255)
                                                                                                                      tAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGGGATAAGTCGCACCGGTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCgg                                 <  1:64361/150‑1 (MQ=255)
                                                                                                                        gCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCa                               >  2:1192431/1‑150 (MQ=255)
                                                                                                                        gCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCa                               >  1:1917412/1‑150 (MQ=255)
                                                                                                                         cGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGCTAAGTCGCACCGTTTGCCCAATCCACCACACTCTTCCGTTCCCTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCAt                              >  1:1581733/1‑150 (MQ=255)
                                                                                                                         cGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGCCCCGCGGCAt                              >  1:1125478/1‑150 (MQ=255)
                                                                                                                         cGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCAt                              >  1:13089/1‑150 (MQ=255)
                                                                                                                          gATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATc                             <  2:1651636/150‑1 (MQ=255)
                                                                                                                          gATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATc                             <  2:1275966/150‑1 (MQ=255)
                                                                                                                               cGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGa                        <  2:58594/150‑1 (MQ=255)
                                                                                                                                gTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGAt                       >  2:1905468/1‑150 (MQ=255)
                                                                                                                                gTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGAt                       >  1:1067523/1‑150 (MQ=255)
                                                                                                                                gTTCGCCCCCAGTGCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGAt                       <  2:1275151/150‑1 (MQ=255)
                                                                                                                                 ttCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATa                      <  2:1159927/150‑1 (MQ=255)
                                                                                                                                   cGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACg                    >  2:852926/1‑150 (MQ=255)
                                                                                                                                   cGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACg                    >  1:1160618/1‑150 (MQ=255)
                                                                                                                                     cccccAGTTCGCGCAGCACGTTCGTGGCGATGTGATAAGTCGCACGGTTTGCACAATGCAGCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAAATCGATATAGCGACCCGCGGCATCAACGATACggc                  <  2:1463794/150‑1 (MQ=255)
                                                                                                                                      ccccAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggct                 <  1:1636445/150‑1 (MQ=255)
                                                                                                                                       cccAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggctg                <  2:487421/150‑1 (MQ=255)
                                                                                                                                       atcacatCGCGCCGAACGTGCTGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggctg                >  1:1243298/7‑150 (MQ=255)
                                                                                                                                         cAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggctggc              <  2:1959996/150‑1 (MQ=255)
                                                                                                                                         cAGTTCGCGCAGCACGTTCGGCGCGAAGTGATAAGTCGCACCGGTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggctggc              <  1:172160/150‑1 (MQ=255)
                                                                                                                                            ttCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAGCGATACGGCTGGCttt           >  1:1577925/1‑150 (MQ=255)
                                                                                                                                             tCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTa          >  2:506473/1‑150 (MQ=255)
                                                                                                                                                  cagcaCGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGGTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCACCGATACGGCTGGCTTTACCCAg     >  2:1331289/1‑150 (MQ=255)
                                                                                                                                                   agcaCGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGt    <  1:1084285/150‑1 (MQ=255)
                                                                                                                                                     caCGTTCGGCGCGATGTGATAAGTCGCTCCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTc  <  2:1453200/150‑1 (MQ=255)
                                                                                                                                                     caCGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACCATCCACCACAATCTTAAGTTCAATGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATAGGGCTGGCTTTACCCAGTTc  <  1:1075120/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTC  >  NC_000913/3323343‑3323641

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