Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,487,163	0	C	T	42.6%	58.5 / 34.0	54	intergenic (+200/‑72)	hrpA/ydcF	putative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base C (5/26);  new base T (23/0);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 9.29e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGG  >  NC_000913/1487021‑1487313                                                                                                                                                 |                                                                                                                                                            attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTAAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAAAAATAA‑G‑AATttctcc                                                                                                                                                         >  2:8471/1‑147 (MQ=255)   taATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCa                                                                                                                                                    <  2:21897/150‑1 (MQ=255)     aTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCagag                                                                                                                                                 <  1:128841/151‑1 (MQ=255)       tCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAC‑G‑AATTTCTCCtttttgag                                                                                                                                                 >  2:222659/1‑144 (MQ=255)          aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTg                                                                                                                                               <  2:688816/148‑1 (MQ=255)           aaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTcc                                                                                                                                           <  1:802520/151‑1 (MQ=255)              cGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAA‑G‑ATTTTCTCCTTTCTTAGGGCTcc                                                                                                                                           >  2:939115/1‑148 (MQ=255)               ggTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGACATAC‑A‑ATTTCTTTCTTTTCTAGTGCt                                                                                                                                             >  1:417692/1‑145 (MQ=255)                       tAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGAAAAAC‑G‑ATTTTCTCCTTTTAGAGTGCTCCGCCTCTCACtat                                                                                                                               >  2:209686/1‑151 (MQ=255)                         aaaaTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGATAAAGACAAAA‑G‑ATTTTCTCCCTTTAGAGGGGTCCGCTTCTCTCtttt                                                                                                                              >  1:953595/1‑150 (MQ=255)                           aaTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAA‑A‑AATTTCTCCTTTTTTAGTGCTCCCCCTCTCCCtatt‑at                                                                                                                           >  1:567736/1‑150 (MQ=255)                                  gAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCa                                                                                                                   <  2:41624/151‑1 (MQ=255)                                    aTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGACAAAA‑G‑AATTTCTCCCTTTCGAGGGCTCCCCCTCTCACtattttttccccc                                                                                                                     >  2:46654/1‑139 (MQ=255)                                      tCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTa                                                                                                                <  1:350035/150‑1 (MQ=255)                                       cAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAATACAA‑ATTTCTTCCTTT‑TTAGGGCTTCGCTCTTCTCTATT‑ATTTCACTCAGtt‑tt                                                                                                             >  1:166602/1‑151 (MQ=255)                                           aaaaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTa                                                                                                        <  1:498749/151‑1 (MQ=255)                                                 cAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGaa                                                                                                  <  1:704284/151‑1 (MQ=255)                                                    ttCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAA‑G‑AATTTCTCCTTTTTTAGAGCTCCCCCTCTCCCTATT‑TTTTCACGCCGTA‑TTT‑TTAAagg                                                                                                    >  2:571184/1‑146 (MQ=255)                                                     tctcTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAA‑C‑AATTTTTTCTTTTTGAGGGCGTCCCTTTTCACTATT‑TTCTCACCCAGTA‑TTT‑TTAAaggga                                                                                                  >  2:606232/1‑145 (MQ=255)                                                        ctgctgCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGAt                                                                                              <  1:545165/148‑1 (MQ=255)                                                          gctgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAGGACATAC‑G‑ATTTTCTCCTTTTCGAGGGGTCCCCCTCTCCCTCTT‑TTTTTCCCCCGGA‑TTC‑TTAAAGGAACGATaa                                                                                            >  2:715272/1‑148 (MQ=255)                                                           ctgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAC‑A‑ATTTTCTTCTTTTAGAGAGGTCCCCCTTTCCCTATT‑TTCTCCCGCAGGATTTT‑TTAAGGGAA‑GATaa                                                                                            >  2:132063/1‑147 (MQ=255)                                                             gcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAaggag                                                                                        <  1:503863/149‑1 (MQ=255)                                                              cAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGaa                                                                                     <  1:835052/151‑1 (MQ=255)                                                              cAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGAAATAA‑A‑AATTTTTTCTTTTAGAGGGGTCCGCTTCTCACTATT‑TTTTTACGCAGGA‑TTT‑TTTAGGGGAAGAAAaggagg                                                                                       >  1:226682/1‑149 (MQ=255)                                                               aaGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAAc                                                                                    <  1:604815/151‑1 (MQ=255)                                                                   tGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAGAAAA‑G‑ATTTTTTCCCTTCAGAGGGCTTCCCCCCTCCCTATT‑ATCTCACCCCGTA‑TTCTTTAAGGGAAGGTAAAggggga                                                                                      >  2:597808/1‑146 (MQ=255)                                                                    gAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑Tg                                                                               <  2:380700/150‑1 (MQ=255)                                                                     aataatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGaa                                                                             <  2:713019/151‑1 (MQ=255)                                                                       taatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGCTAAAGGAATAA‑A‑ATTTTTTCCTTTTTGAGTGGTCCGCGTCCTCCTCTT‑TTTTTTCGCCGAA‑TTC‑TTTAGGGGAACATAAAGAGGAAACA‑TGaaa                                                                            >  1:733739/1‑149 (MQ=255)                                                                        aatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACa                                                                           <  1:740323/150‑1 (MQ=255)                                                                         atGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGACAAAA‑A‑AATTTCTCCTTTCTTAGTGGTCCGCCTCTCACTATT‑TTTTCACCCCGGA‑TTC‑TTTAAGGGACGATAAAGAGGGAACC‑TGAACa                                                                           >  2:730707/1‑149 (MQ=255)                                                                           gTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑c                                                                      <  2:81036/151‑1 (MQ=255)                                                                            tttATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAA‑G‑ATTTTCTCCTTTTAGAGTGCTCCCCCTCTCCCTATT‑TTCTCCCCCAGTA‑TTCTTTAAGGGAA‑GAGAAGGAGGaaaaacgaaaa                                                                            >  2:900862/1‑137 (MQ=255)                                                                             ttATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑cc                                                                     <  1:358154/150‑1 (MQ=255)                                                                                  tACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGCTTAAGAAAAAA‑G‑ATTTTCTCCTTTTTGAGTGGTTCCCTTTTCTCTATT‑ATTTCTCTCAGTATTTT‑TTAAGGGAA‑GATTAAGAGGAAGCA‑TTAAACTTA‑cccccctt                                                               >  2:42526/1‑147 (MQ=255)                                                                                     aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAGGAAATAA‑A‑AATTTTTTCTTTCTGAGTGGTCCCCTTTTTCATATT‑TTTTCAACCCGGA‑TTT‑TTTTGGGGAAGATAAAGGGGGAACA‑TGAAAAATA‑CCCCGctt                                                               >  1:244879/1‑148 (MQ=255)                                                                                     aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGCTTAAGACAAAA‑GCTTTTCCTCCTTCTGGAG‑GTTCCGCTTTTCCCCATTCTTCTC‑CGCCGGT‑TACTTTAAGGGAA‑GGTAAGGGAGGACCA‑TGAAAAATA‑CCCCCTTTcc                                                             >  2:150062/1‑150 (MQ=255)                                                                                       cATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAA‑G‑ATTTTCTCCTTTTTGAGTGGTCCCCCTCTCACTATT‑TTTTCACCCCGGATTTT‑TTAAGGGAA‑GAAAAAGAGGAAAAATTAAAAATAA‑CCCCTTTTCCgc                                                           >  2:932350/1‑150 (MQ=255)                                                                                          ttCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGACAAAA‑G‑ATTTTCCCCCTTCTTAGGGGCCCCCTTCTCCCCATT‑TTTTTCCCCCGTATTTT‑TTAAGGGAA‑GATAAGGAGGAAACA‑TTAACATTACCCCCGTTTc                                                              >  2:58319/1‑145 (MQ=255)                                                                                               tAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTc                                                    <  1:252305/149‑1 (MQ=255)                                                                                                    gATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAATAA‑C‑ATTTTCTCCCTTCAGAGTGCTCCCCTTCTCACTTTT‑ATTTCCCCCAGAA‑TTC‑TTAAAGGAACGATAAAGAGGAAACA‑AGAAAAATA‑ACCCCCTTCCGACGCgttttcc                                                 >  2:583571/1‑143 (MQ=255)                                                                                                    gATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAA‑G‑GTTTTCTCCCTTTAAAATGGTCCGCTTTTCACTATT‑TTTTCACCCCGGATTTT‑TTAAGGGAA‑GATAAAGAGGGGACA‑ACAAAAATA‑ACCCGTTTTCGACCCGTTTGCCgg                                               >  2:722400/1‑149 (MQ=255)                                                                                                     atatGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCaa                                            <  2:315300/151‑1 (MQ=255)                                                                                                        tGTCAGGCTTGCGGAGTGGCGGTTAAGGAAAAAC‑G‑ATTTTCTCCTTTTTGAGGGCTCCCCCTCTCACTATT‑TTCTCACCCCGGA‑TTC‑TTAAAGGAAAGATAAAGAGGGAACA‑TGAACATTA‑CCCCGCTTCCGACGCTTTCGCCCGCAAc                                           >  2:241752/1‑149 (MQ=255)                                                                                                            aGGCTTGCGGAGTGGCGGTTTAGGAAATAA‑A‑ATTTTCTCCTTTCTGAGTGCCCCCCTTCTCTCTCTT‑TTCTCACCCAGTATTTT‑TTAAGGGAA‑GAAAAAGAGGGAACA‑TGAACATTA‑CCCCGCTTCCTACGCtct                                                     >  1:123585/1‑133 (MQ=255)                                                                                                             ggCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGAt                                    <  2:835599/151‑1 (MQ=255)                                                                                                                     gAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATaaa                            <  1:830421/151‑1 (MQ=255)                                                                                                                         ggCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGtt                        <  2:16362/151‑1 (MQ=255)                                                                                                                         ggCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGt                         <  2:211120/150‑1 (MQ=255)                                                                                                                            ggTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATc                     <  2:400617/151‑1 (MQ=255)                                                                                                                              ttAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTAGCCGGCAACTATAGATGCCATAAATGTTATCg                    <  2:328856/150‑1 (MQ=255)                                                                                                                                aaGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGAc                 >  2:153461/1‑151 (MQ=255)                                                                                                                                   gACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGt              <  1:950746/151‑1 (MQ=255)                                                                                                                                    aCATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtg             >  1:379919/1‑151 (MQ=255)                                                                                                                                                tcctccTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGgcgc     >  1:880483/1‑149 (MQ=255)                                                                                                                                                 cctcctTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAgg  <  2:569148/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                            ATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAC‑G‑ATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATT‑ATCTCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCA‑TGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGG  >  NC_000913/1487021‑1487313

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