Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,092,703 A→G F41L (TTC→CTC)  frvA ← putative PTS enzyme IIA component FrvA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,092,7030AG96.2% 152.2 / ‑2.1 52F41L (TTC→CTC) frvAputative PTS enzyme IIA component FrvA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (20/30);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

TGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATT  >  NC_000913/4092466‑4092937
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
tGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGTGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTc                                                                                                                                                                                                                                <  1:628988/250‑1 (MQ=255)
  cGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCCGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCggg                                                                                                                                                                                                                             >  2:564372/1‑250 (MQ=255)
     cacaGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACAAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGAt                                                                                                                                                                                                                          <  1:374686/251‑1 (MQ=255)
      cacgtttgccgggtactttGTCCTGATCCTCGTCGCCGCTTTGCGGCACGCCGAGGCAGATCAAGCTATTGACGGCTGCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAATTCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGGTGGTGTGAGTCGGCGATa                                                                                                                                                                                                                         <  2:574398/235‑1 (MQ=255)
       caGTGTGCCGATGATTTTGACCTGATCCTCTTCGCAGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGACCCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGGGTGGGTCGGTgtta                                                                                                                                                                                                                         >  2:644443/1‑247 (MQ=255)
          tgtgCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAAc                                                                                                                                                                                                                     >  2:134676/1‑251 (MQ=255)
             gCCGATGATGTTGCCCTGAACCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCATCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGt                                                                                                                                                                                                                  <  2:335649/251‑1 (MQ=255)
               cGATGATTTTGACCTGATCCTCGTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGttt                                                                                                                                                                                                                <  1:660280/251‑1 (MQ=255)
                gatgatTTTGCCCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCCGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAAACGtttt                                                                                                                                                                                                               >  2:628620/1‑251 (MQ=255)
                gatgatTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGtttt                                                                                                                                                                                                               >  1:23997/1‑251 (MQ=255)
                             tGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGt                                                                                                                                                                                                  >  1:655949/1‑251 (MQ=255)
                                 ccTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGcc                                                                                                                                                                                              <  1:161146/251‑1 (MQ=255)
                                   tctTCGCCGCTTTGCGGCACGCCGCGGCAGATCAAGCAATTTACGTCTTCGCCCTCGCTGGCTTTCCAGTCAATCGCCTGCGCTTTGTGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGGCACCGGAACCAAATCAAGTAGGGTGCATTTTTTCTCGCAGCAGCAGGGTTTGCAGGAGATGGTGTGAGTCGGTGATAAAACCGTTTTGTTACGATATTGTCGCCa                                                                                                                                                                                             >  2:456768/1‑250 (MQ=255)
                                      tCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACt                                                                                                                                                                                          <  1:70542/250‑1 (MQ=255)
                                        gccgcTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCt                                                                                                                                                                                       <  1:86141/251‑1 (MQ=255)
                                                 cGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGATTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGATTCAGAACgg                                                                                                                                                                              >  2:149589/1‑251 (MQ=255)
                                                      cGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATa                                                                                                                                                                          >  1:618711/1‑250 (MQ=255)
                                                              cAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCAtt                                                                                                                                                                 >  1:271727/1‑251 (MQ=255)
                                                               agagacAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGGTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACAGCGGAACCACATCCAGTAGAGTGCAGTGTTTCGCGCAGCAGCAGGGTTTGCAGGAGCGGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTg                                                                                                                                                                <  1:421438/246‑1 (MQ=255)
                                                                   ccAGCAATTAACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGCGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGCTTGCAGGAGCTGGGGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCt                                                                                                                                                            >  2:504204/1‑251 (MQ=255)
                                                                    cAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTg                                                                                                                                                           <  1:134676/251‑1 (MQ=255)
                                                                    cAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGAGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAAAGGAATAAGGGCCATTGCCCTg                                                                                                                                                           >  2:110008/1‑251 (MQ=255)
                                                                                     gCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATAc                                                                                                                                          >  1:321416/1‑251 (MQ=255)
                                                                                          cGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGCACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCtt                                                                                                                                     >  1:562812/1‑251 (MQ=255)
                                                                                                    ccAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAagag                                                                                                                           <  1:560102/251‑1 (MQ=255)
                                                                                                        tCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGc                                                                                                                       >  1:400662/1‑251 (MQ=255)
                                                                                                           aTAGCCTGCGCTTTGCGGGCGAAGAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCAt                                                                                                                    <  2:240155/251‑1 (MQ=255)
                                                                                                                cTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTact                                                                                                                >  1:26399/1‑248 (MQ=255)
                                                                                                                  gcgcTTTGCGGGGGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTAAGAACGGAATAAGCGCCATTGCCCTGAATATGAAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATc                                                                                                               <  2:535460/249‑1 (MQ=255)
                                                                                                                  gcgcTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCgg                                                                                                             <  1:535083/251‑1 (MQ=255)
                                                                                                                  gcgcTTTGCGGGAGAATAATACGAACTGGTGTTTAACGCAGGCGCTTTTTCCGTGCGGTACGGCGACACCGGAACCAAATCCAGTAGAGTTCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCTTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATGCAGGTCAATACAGCTTGCAGTAAGATCTGCCATTAATc                                                                                                               <  2:563767/249‑1 (MQ=255)
                                                                                                                       ttGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTc                                                                                                         <  1:266909/250‑1 (MQ=255)
                                                                                                                            ggCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATaaa                                                                                                   <  1:708813/251‑1 (MQ=255)
                                                                                                                                aataatACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatat                                                                                               <  1:58587/251‑1 (MQ=255)
                                                                                                                                 attattcGAACGGGTGTTTAACGCAGGCGCTTTTACCGGGCTGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTTTGAGTCGGAGATAAAACTGTTTTTTAACGCTATTGTCGACAACTGCTGCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATAAAGCTGGCAGTAAGAGCTGCAATTAATCGGACTCCATAAAatata                                                                                              <  2:512183/245‑1 (MQ=255)
                                                                                                                                        gAACGTTTGTTTAACGCGAGCGCTTCTACCGGTCGGCTCAGCGACACCGGAACCAAATCCAGGAGAGTGCATTTTTTCGCGCAGCAGCAGGTTTTGCAGGAGCTGGTGTGAGGCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTACTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATAtttgact                                                                                       <  2:157694/251‑5 (MQ=255)
                                                                                                                                          aCGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCTACACCGGAACCAAATCCAGTAGAGTGCTTTTTGTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGAGAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTaagca                                                                                     <  2:497265/251‑4 (MQ=255)
                                                                                                                                           cGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCa                                                                                     >  1:107538/1‑250 (MQ=255)
                                                                                                                                            ggTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAAt                                                                                   <  1:255885/251‑1 (MQ=255)
                                                                                                                                                    aaCGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCCGTAAGCGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtgg                                                                           >  1:81495/1‑251 (MQ=255)
                                                                                                                                                                tttACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTa                                                                <  1:325450/250‑1 (MQ=255)
                                                                                                                                                                tttACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTa                                                                <  1:628620/250‑1 (MQ=255)
                                                                                                                                                                 ttACCGTGCGGCACGGCAACACCGGAACCAAATCCAGTAGAGGGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTAAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTaa                                                               <  2:140778/250‑1 (MQ=255)
                                                                                                                                                                    ccGTGCGGATCGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGATTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTGGCAGTAAGAGCGGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtga                                                           <  2:291092/251‑1 (MQ=255)
                                                                                                                                                                     cGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtgat                                                          >  1:491458/1‑251 (MQ=255)
                                                                                                                                                                                             cAAATCCAGTAGAGTGCATTTTGTCGCGCAGCAGCATGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACaa                                  <  2:4022/251‑1 (MQ=255)
                                                                                                                                                                                                       agagTGCCTTTTTTCGCGCAGCAGCAGGGTTTTCAGGAGCCTGTGTGAGTCGGTGATAAAACCGTTTTGGAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTg                        <  1:644443/251‑1 (MQ=255)
                                                                                                                                                                                                       agagTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAttt                         <  1:21878/250‑1 (MQ=255)
                                                                                                                                                                                                                tttttCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACagag                <  1:677179/250‑1 (MQ=255)
                                                                                                                                                                                                                tttttCGCGAAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGAGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAACATCTATTTTACTCAATTGTAATGTTGGAGAGATTAAGTGATGAGAGGCGATATGTCTGTAAACAATGAAAATTTACGACagagg               >  2:493494/1‑250 (MQ=255)
                                                                                                                                                                                                                           gcagcagGGTGTGAAGGAGCGGGGGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCTCTGAATATTCAGGGCAATACAGCTTGAAGTAAGAGCTGGCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGa                                                                       <  1:404647/184‑1 (MQ=255)
                                                                                                                                                                                                                             agcagGGTTTGCAGAAACTGGTGTGAGTCGGTGATAAAACAGTTTTGTAACGCTATTGTCGACAACTGCTTCAGAAAGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGGGAAAGGACAGAtt  >  2:236604/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
TGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATT  >  NC_000913/4092466‑4092937

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑