Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,487,153	0	C	A	38.2%	119.3 / 24.7	55	intergenic (+190/‑82)	hrpA/ydcF	putative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base C (5/29);  new base A (21/0);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 7.81e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGAC  >  NC_000913/1487006‑1487298                                                                                                                                                      |                                                                                                                                                           aGCGCCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATa                                                                                                                                                        <  2:939128/151‑1 (MQ=255)               tattaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCtt                                                                                                                                         <  2:764721/150‑1 (MQ=255)                   aaTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCag                                                                                                                                     <  2:1179182/150‑1 (MQ=255)                    aTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCagag                                                                                                                                   <  1:1245451/151‑1 (MQ=255)                      tCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGt                                                                                                                                  <  2:1199690/150‑1 (MQ=255)                        tAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTTAAGAA‑ATAAGAA‑TTTTTTCCttttc                                                                                                                                      >  2:436433/1‑143 (MQ=255)                        tAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTATAGAA‑AAACG‑A‑GTTTCTCCTTTCata                                                                                                                                    >  1:406398/1‑144 (MQ=255)                         aaaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCt                                                                                                                               <  1:131136/150‑1 (MQ=255)                                  aGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑GTTAGGGA‑ATACG‑ATTTTCTTCCttttttattcggtcc                                                                                                                            >  2:339528/1‑133 (MQ=255)                                  aGGGTAAAAATTCAGGAATCCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTTAAGAA‑AAAAG‑A‑ATTTCTCCTTTTCGAGTGGTcc                                                                                                                             >  1:287473/1‑143 (MQ=255)                                     gTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACta                                                                                                                  <  2:1100929/151‑1 (MQ=255)                                         aaaTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGGTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TAATGGAA‑ATAGG‑A‑TTTTCTTCTTTcttagaggg                                                                                                                               >  1:1478869/1‑126 (MQ=255)                                             tCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑GTAAAGAA‑ATAAA‑A‑TTTTTTTCTTTTTTAGTGCTCCCCTTCTCCCtctttt                                                                                                              >  2:829726/1‑145 (MQ=255)                                               atgAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCAc                                                                                                        <  1:1214454/148‑1 (MQ=255)                                                  aaTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAg                                                                                                    <  2:173055/151‑1 (MQ=255)                                                   aTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGt                                                                                                   <  1:1121521/151‑1 (MQ=255)                                                   aTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAAGAA‑ATACA‑A‑TTTTTTCCTTTTTTAGAGGTCCGCCTCTCCCtcttttttcccaca                                                                                                      >  2:1089105/1‑139 (MQ=255)                                                      cAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑t                                                                                                <  1:228364/150‑1 (MQ=255)                                                      cAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑t                                                                                                <  1:1278972/150‑1 (MQ=255)                                                       aGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGG‑TTAAGAAA‑AAAAA‑A‑TTT‑TTCCTTTTAGAGAGGTCCCCTTCTCCCTATTATT‑TCACCCACTA‑t                                                                                                >  2:1104398/1‑149 (MQ=255)                                                        gAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑t                                                                                            <  1:734486/151‑1 (MQ=255)                                                            aaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGATTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAg                                                                                        <  1:233325/151‑1 (MQ=255)                                                                   ttCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTTAGGAA‑ATAAA‑A‑ATTTCTTCTTTCTTAGTGGGTCGCCTCTCCCTATTATT‑TTACGCCGCA‑TTC‑tt                                                                                           >  2:1098729/1‑141 (MQ=255)                                                                    tctcTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTAGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGAt                                                                                >  2:606696/1‑151 (MQ=255)                                                                     ctctGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GGTTTAAGAAA‑AAACC‑A‑‑TTCCTCCCTTTTAAGGGGTTCCCCTCTCCCTATTATT‑TCTCGCCGGA‑TTC‑tttatggga                                                                                    >  1:340971/1‑139 (MQ=255)                                                                       ctgctgCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATa                                                                               <  1:1134047/149‑1 (MQ=255)                                                                           tgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTTAAGAA‑AAAAG‑A‑TTTCCTCCCTTTTTGGGGCTTCCCCTCTCCCTATTTTT‑TCACGCCGTC‑TTC‑TTTAAGGGACGATAagg                                                                            >  1:1241718/1‑148 (MQ=255)                                                                             cAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGaa                                                                       <  2:710617/151‑1 (MQ=255)                                                                                     ataatGTTTATCTTCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑GTTAAGAA‑AAAAA‑A‑TTTTCTCTTTTTATCAGGGTCCGCCTCTCCCTATTATT‑TTTCTCCAGA‑TTT‑TTTAGGGGACAATAaggagggaacaaaaa                                                                >  2:413083/1‑141 (MQ=255)                                                                                     ataatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAc                                                              <  1:1202252/151‑1 (MQ=255)                                                                                        atGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑Att                                                          <  1:839081/151‑1 (MQ=255)                                                                                          gTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑c                                                       <  2:725968/151‑1 (MQ=255)                                                                                              aTCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑cccc                                                    <  1:1069277/150‑1 (MQ=255)                                                                                                  aCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGttt                                                <  2:1363389/150‑1 (MQ=255)                                                                                                  aCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTc                                               <  1:505011/151‑1 (MQ=255)                                                                                                     gCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑GTAAAGACAAAACA‑A‑TTTCTTCC‑TTCTGAGTGGCGCGCCTCTTCCTATTATT‑TCTCACCGGC‑TTCTTTAGGGGAAGGAAAAGGGGGAAAAATA‑ACA‑ATTA‑CCCCCtt                                                 >  1:199101/1‑147 (MQ=255)                                                                                                      cATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAAGAA‑ATACA‑A‑TTTTCTCCTTTTTTAGTGGTCCGCCTCTCCCTTTTATT‑TCACCCAGAATTTT‑TTAAGGGAA‑GATAAGGAGGGAACAAG‑AAA‑AATA‑cccccctttc                                              >  2:874804/1‑142 (MQ=255)                                                                                                       aTTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAAGAA‑ATAAG‑A‑TTTTCTCCCTTCTTAGTGCTCCGCTTCTTACTATTTTT‑TCACGCCGGA‑TTC‑TTAAGGGGAAGAAAAAGAGGAAACACGAAAA‑ATAA‑CCCCTTTTccac                                            >  1:1026445/1‑148 (MQ=255)                                                                                                              tAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTc                                     <  2:878749/149‑1 (MQ=255)                                                                                                              tAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTc                                     <  2:30811/149‑1 (MQ=255)                                                                                                               aaaaGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGcc                                  <  2:1052229/151‑1 (MQ=255)                                                                                                                aaaGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTTAAGAA‑AAAAA‑A‑ATTTCTCCTTTTTCAGGGCTCCCCCTCTCCCTATTTTT‑TCACGCCGGATTTT‑TTAAGGGAA‑GATAAAGAGGGAACATG‑AAA‑AATA‑ACCCCCTTTCGACGCctt                                      >  1:444237/1‑146 (MQ=255)                                                                                                                aaaGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAAGAA‑AAAAC‑A‑ATTTCTCCTTTTTGAGGGTTCCCCCTCTCACTCTTTTCTTCCCCCGGAT‑TTT‑TTAAGGGAA‑GATAAAGAGGGAACAAG‑AAA‑AATATCCCCCTTTCC‑ACGCtttttcc                                  >  1:728485/1‑146 (MQ=255)                                                                                                                aaaGATATGTCAGGCTTGCGGAGTGGC‑GG‑GTTAGGAC‑AAACG‑A‑ATTTTTCCTTTTTTAGGGGTCCCCCTCTCCCTTTTATT‑TTCCCCCGTA‑TTT‑TTTATGGGACGATAAAGAGGGACCAAG‑AAC‑ATTATCCCCCTTTCC‑ACGCttt                                      >  2:779095/1‑146 (MQ=255)                                                                                                                  aGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTTAAGAA‑ATACG‑A‑ATTTTTCCCTTTCGAGTGTTCCGCCTCTCACTTTTTTC‑TTACCCAGGA‑TTC‑TTTAGGGGACGATAAGGAGGAAACATG‑AAC‑ATTA‑CCCCCCTTTCCACGCTTTCGCCgg                                >  2:1097550/1‑150 (MQ=255)                                                                                                                    atatGTCAGGCTTGCGGAGTGGC‑GG‑TTTAAGAA‑AAAAC‑C‑TTTTCTCCCTTTATAGTGCTCCGCCTTTCCCTATTTTT‑TCACGCAGTA‑TTT‑TTTAGGGGAAGATAAAGAGGGAACATG‑AACAATAA‑CCCCCTTTccccgcttt                                       >  1:34093/1‑135 (MQ=255)                                                                                                                    atatGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAA‑ATAAA‑A‑ATTTCTCCTTTTCGAGGGCTCCGCTTCTCCCTATTTTT‑TTCCCCCGTATTTT‑TTAAGGGAA‑GATAAGGAGGAACCATG‑AAA‑AATA‑CCCCCTTTTCCCCGCTTTTGcc                                  >  1:1094450/1‑146 (MQ=255)                                                                                                                              cTTGCGGAGTGGC‑GG‑TTATGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGTGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGAt                     <  1:380572/149‑1 (MQ=255)                                                                                                                               ttGCGGAGTGGC‑GG‑TTAAAGAA‑AAAAG‑A‑TTTTCTCCTTTCAGAGGGGTCCGCTTCTCCCTATTATC‑TTCCGCAGTG‑TTT‑TTATGGGAACGATAAGGAGGAACCATT‑AAC‑AATA‑ACCCGTTTTCGACGCTTTTCCCGGCAACTATAGATTcc                  >  2:408551/1‑151 (MQ=255)                                                                                                                                tGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGc                   <  1:240184/149‑1 (MQ=255)                                                                                                                                       tGGC‑GG‑TTAAGGAA‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGAGGCCATaa              >  2:135072/1‑147 (MQ=255)                                                                                                                                         gC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTa        <  2:971169/151‑1 (MQ=255)                                                                                                                                          c‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTAt       >  1:805053/1‑151 (MQ=255)                                                                                                                                               ttAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTCTTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATc      >  1:1151283/1‑149 (MQ=255)                                                                                                                                                 aaGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGAc  <  2:1465494/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                           AGCGCCGGGTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGC‑GG‑TTAAGGAC‑ATACG‑A‑TTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATC‑TCACGCAGTA‑TTC‑TTAAGGGAACGATAAGGAGGAACCATG‑AAC‑ATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGAC  >  NC_000913/1487006‑1487298

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