Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,502,7790AT30.2% 72.2 / 13.6 43L218* (TTA→TAA) insGIS4 family putative transposase
Reads supporting (aligned to +/- strand):  ref base A (8/22);  new base T (12/1);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.13e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTGTATCCGGT  >  NC_000913/4502570‑4503015
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
gATGGCAGACTTTTCCTTTGCGCGTCACGGGCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTGTCGTGCCGGCGGGCTGTTTTTCAGCTTCACCAGATGATCTCCGGTACCCAGTTTGCTGATCTCTTCATATTGCGCCCCCTTTCGGAGAGGTATCATCCAGTGGCGGTGGTCTCCCGCCAGTCC‑CCAGGCATTTAGCAG‑GCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCaga                                                                                                                                                                                                       <  1:296171/251‑1 (MQ=255)
       gACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGTGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTAt                                                                                                                                                                                                <  2:268893/251‑1 (MQ=255)
            ttCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTTTTATCCCGCCCGGCTTCCAG‑CAATTAAAAG‑TCCCAGAGAGTAATAA‑CCTTTATTCATTAAAGTCAAAGTGGTATCGCC‑g                                                                                                                                                                                          >  1:26435/1‑251 (MQ=255)
              ccTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCCGGCCTTTTACCG‑GCCCCGTGGGGAATAA‑ACCTTATTCATTAACGTCCGAGTGGTATCGCC‑ggg                                                                                                                                                                                        >  1:142143/1‑250 (MQ=255)
                  tGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTg                                                                                                                                                                                     <  1:430784/250‑1 (MQ=255)
                      cgTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCta                                                                                                                                                                                <  1:389405/251‑1 (MQ=255)
                      cgTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCGTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGGGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCta                                                                                                                                                                                <  1:141296/251‑1 (MQ=255)
                       gTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTATTTTCCCCCCAGGCT‑CCCGGCATTTTACAG‑TCCCACAGAGGAATAA‑CCCTTATCCCTTTACGTCAGAGTGTTATCGCC‑CGGTTGTTCtat                                                                                                                                                                               >  1:418269/1‑251 (MQ=255)
                            ggTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCCCCCCGCT‑CCCGGCCTTTAAAAG‑TCCCCGGGAGTAATAA‑ACTTTTTCCCTTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATaaatt                                                                                                                                                                          >  1:236604/1‑248 (MQ=255)
                                  cagGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCAAGTGGCGGTGGTCTCCCCACACGCT‑CCCCGCATTTTACAG‑TCCCCGAGAGTAAAAA‑CCTTTTTTCCTTAACATCAGGGTGTTATCGCC‑GGTTTGttgttttatttgctcc                                                                                                                                                                     >  2:416804/1‑236 (MQ=255)
                                   gggCGGGCAGTCCCTTCATTTCCCAGTCCCTGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCGCCCTTTCCGAGAGGTATCACCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATGTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGc                                                                                                                                                                   <  2:15242/250‑1 (MQ=255)
                                   agGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCCCCCGGCC‑CCAGGCATTTTACAG‑TCCCCGTGGGTAATAA‑CCTTTTTCCCTTAACGTCCGAGGGTTATCGCC‑GGGTTTTTCTATTAGTTGCTCAGc                                                                                                                                                                   >  2:127948/1‑251 (MQ=255)
                                      cGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCCTCCAGTGGCGGTGTTTTCCCGCCCGGGT‑TCAGGCCTTTTACAG‑TCCCCATGAGTAATAA‑ACCTTATCCCCTAACGTCAGAGGGGTTTCGCC‑GGGTTTTTCTATCAGTTGCTccgcca                                                                                                                                                                 >  2:81535/1‑245 (MQ=255)
                                       gggCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCAGGAT‑CCAGGCATTTAACAGACCCCAG‑GAGGAAATA‑CCTTTATCCATTAACATTCGAGTGGTTTCGCC‑GGGTTGTTCTCAAAGATGCTTAGCAAGc                                                                                                                                                               >  1:282364/1‑251 (MQ=255)
                                       gggCAGTCACTTCATTTCCCAGTCCCGGCCACGTTTTTCGTGCCGGCGGGCTGGTGTTCAGCTTCACCAGATGATCGCCGTTACCCAGTGTTCTGATCTCTTCATATTGCGCTCCCGTTCGGAGAGGTATCATCCAGGGGCGGTGGTCTCCCGCCAGGCG‑CCAGGCATGTAACAG‑TCCCAGTGAGGAATAA‑CCTTTACCCAGTAACGTCAGAGTGTTATCGcc                                                                                                                                                                                            <  2:437472/223‑1 (MQ=255)
                                            gTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCCGTGGGGAATAA‑CCTTTATCCCTTAAAGTCCGAGGGTTATCGCCCGGTTT‑TTCTATAAGTTGCTCAGCCAGCTCA‑tt                                                                                                                                                         >  1:226009/1‑251 (MQ=255)
                                                      ttCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGtt                                                                                                                                                <  2:373244/250‑1 (MQ=255)
                                                                 ggCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATc                                                                                                                                                                                               <  1:78907/194‑1 (MQ=255)
                                                                   ccACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCCGGCA‑CCAGGCCTTTTACAG‑ACCCCGGGAGGAATAAACCTTT‑TACATTAACGTCAGAGTGTTATCGCC‑GGGTTGTTCTTTAAATTGGTCAGCAAACCCA‑TTTTTTCTGTTCTTCATCGTGcccc                                                                                                                                  >  1:17925/1‑249 (MQ=255)
                                                                            tCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTAGAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCa                                                                                                                          <  2:267770/250‑1 (MQ=255)
                                                                                 ccTGCGGACTGGGTTTCAGCTTCACCAGCTGATCGCCTTTACCCCGTTTTCTGATCTCTTCATATTGCGCTCCATTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCCCCCGGCA‑CCCGGGCTTTTACAA‑TCCCCATGAGGGCTAA‑CCTTTTTCCATTTAAGTCAACGTGTTTTCACA‑CGTTTTTTCTCTAACTTTCTCAGTAACCTCC‑TTTTCGCTGTTCTTCATCGTGCCGAAggttggcaccggc                                                                                                                    >  2:266802/1‑240 (MQ=255)
                                                                                      gggCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCCCCCTGCC‑CCCGGGCTTTTACAA‑TCCCCGGGAGTAATAA‑ACTTTTTCCCTTAACGTCCGGGGGGTATCCCC‑CGTTTTTTCCATAAGTTGCTTCGCCAGCTCAATTTTC‑CTTTTCTTCATTCTGCCGacagcgtccacccgccgcac                                                                                                               >  2:221128/1‑232 (MQ=255)
                                                                                       ggCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGa                                                                                                              <  1:416804/251‑1 (MQ=255)
                                                                                               ttCAGCTTCACCAGATGATCGCCTTTACCCCGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTTCCCCCAGGCT‑CCAGGCATTTTACAG‑TCCCAGGGCGTAATTA‑ACTTTTTTCATTAACGTCAGAGGGTTTTTGCC‑GGTTTTTTTTTTAAGtt                                                                                                                                                                          >  2:192240/1‑184 (MQ=255)
                                                                                               ttCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTGCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGtt                                                                                                                                                                          <  1:192240/184‑1 (MQ=255)
                                                                                                 cAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCa                                                                                                    <  2:226009/251‑1 (MQ=255)
                                                                                                   ccTTCCCCAGATGATCGCCTTTACCCTGTTTTCTGATCTCTTCATATTTCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCATGCT‑CCAGGCATTTAACAG‑TCCCAGTGATTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGt                                                                                                  <  2:322613/250‑1 (MQ=255)
                                                                                                            gatgatCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCCGGCT‑CCAGGCATTTAACAG‑TCCCCGTGAGTAATAA‑ACTTTTTCCATTAACGTCAGAGTGTTATCGCC‑GGTTTTTTTTATTAATTGCTCAGCAAAGTCC‑TTTTTGCTGTTCTTCATCCTGCCGAAAGCTGCCGCCGTCAGCAGATGGCTGGTCAGTTCCATCttg                                                                                         >  1:122165/1‑249 (MQ=255)
                                                                                                                  cGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCCCCAGGCG‑CCAGGCATTTAAAAG‑TCCCAGTGAGTTATAA‑ACTTTTTCCATTCACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTTCTCAGCAAGATCC‑TTTTTGCTGTTTTTCCTTGTGCCCAAGGCTGCAGCCCTCAGCCCGTGGCTGGTCAGTTCCCTCTGGccgcac                                                                                   >  1:357474/1‑246 (MQ=255)
                                                                                                                        tACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCGCA‑TTTTCGCTGGTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCAtttt                                                                              <  2:236604/250‑1 (MQ=255)
                                                                                                                                    tGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCCCCCGGCC‑CCAGGGATTTTACAA‑TCCCAGTGGGGAATAA‑CCTTTTTCCATTAACGTCAGAGGGTTATCGCC‑GGTTTTTTTTATAAGTTTCTCAGCAAGCTCA‑TTTTTGCTTTTTTTCTTTGTGCCGCAAGCTGCAGCCGTCAGCAGATGGGTGGGCCGTTTCATCTGGCAGACCCtttttaactccgggtgg                                                                 >  2:462309/1‑238 (MQ=255)
                                                                                                                                      aTCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCCGTGAGTAATAA‑CCTTTATCCCTTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCggggag                                                                 >  1:8719/1‑246 (MQ=255)
                                                                                                                                       tctcTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCCCCCAGCT‑CCAGGCATTTTACAG‑TCCCAGTGAGGAATAA‑CCTTTTTCCATTAACCTCCGAGGGGTATCGCC‑GGTTTTTTCTATCAATTTCTCCGCCAGCTCA‑TTTTCGCTTTTTTTTATCTTGCCGCAGGCCGCGGCCGTCAGCAGATGGCTGTTCAGTTCCTTCTGGCGGACCATTTTGTCCTGCGGGTagggg                                                              >  2:178068/1‑248 (MQ=255)
                                                                                                                                                      cTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGGTCTCCCCCCCGGCC‑CCAGGCATTTTACAGATCCC‑GGGAGTAATAA‑CCTTTTTCCATTAACGTCAGCGTGTTATCCCC‑GGGTTGTTCTTTAAGTTTCTCAGCAAGCTCC‑TTTTTGCTGTTTTTCTTCGTGCCGAAGGCTTCCGCCGTCATCAGATGGCTGGTCCGTTTCCTCTCGCCGACCATTTTGTCCTCCGGGGAGAGCGCCCGGTTccccgcc                                               >  2:328408/1‑247 (MQ=255)
                                                                                                                                                            ttCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATgtgt                                          <  1:219275/250‑1 (MQ=255)
                                                                                                                                                            ttCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCAGCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCTGGTTCCCGGCATgtgt                                          <  2:476642/250‑1 (MQ=255)
                                                                                                                                                                   gggtgtcagccAGTGGCGGTGTTCTCACGCCAGGCG‑CCAGGCATTTAACAG‑GCCCAGTGAGTTATAG‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGTAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCCTGTGTCTGGCgtg                                  <  2:351573/242‑3 (MQ=255)
                                                                                                                                                                           cccAGTGGCGGTGTGCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTTATTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTAGCGCC‑GGTTTGTTCTATATGTTGCGCAGCAAGCTCA‑TTGTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGGCAGCTGATGGCTGGGTAGTTCCATCTGGCAGATCATTTTGTCCTGCGGGTAGAGCACCGGGTTCCCGGCATGGGTCTGGCGGGGGAAGGCt                          <  2:293019/250‑1 (MQ=255)
                                                                                                                                                                                 ggCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATc                     <  1:462309/250‑1 (MQ=255)
                                                                                                                                                                                 gccGGTTTTCCCCCGCCCGGCT‑CCAGTCATTTTACAG‑TCCCGGTGAGTAATAA‑CCTTTATCCATTACCGTCAGCGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTTATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGGCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATc                     <  2:366887/248‑1 (MQ=255)
                                                                                                                                                                                              cGCCAGGCT‑CCAGGCATTTAACAG‑TCCCATTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTGTAt       <  1:170449/251‑1 (MQ=255)
                                                                                                                                                                                                ccAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGAGTATcc     <  1:127948/251‑1 (MQ=255)
                                                                                                                                                                                                   ggCT‑CCAGGAATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGCAGGCTGCCTCGTTCTCTGGGGTCTCCggg  >  2:396426/1‑250 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
GATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCT‑CCAGGCATTTAACAG‑TCCCAGTGAGTAATAA‑CCTTTATCCATTAACGTCAGAGTGTTATCGCC‑GGTTTGTTCTATAAGTTGCTCAGCAAGCTCA‑TTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGAAGGCTGCATCGTTCTCTGGTGTATCCGGT  >  NC_000913/4502570‑4503015

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