Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,192	0	G	C	28.8%	99.1 / 12.9	63	intergenic (‑27/+8)	leuV/leuP	tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (9/33);  new base C (0/17);  total (10/53)
Fisher's exact test for biased strand distribution p-value = 4.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAA  >  NC_000913/4606042‑4606337                                                                                                                                                        |                                                                                                                                                  aaaTCTAAGAAGAAACATCTTATGAAATAAAGTTAATTTGTGCCAGCCGGGGC‑ATTTAAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGAGTCTACCAATTCCGCCAACTTAGCACAGTCATCTTACTTTTTTTTATATCg                                                                                                                                                   <  7:1579409/151‑1 (MQ=12) aaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGC‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCg                                                                                                                                                   <  8:503253/151‑1 (MQ=255)  aaTCAAGGAAGAAACAACAACGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                                                                                                                                                  <  1:709562/151‑1 (MQ=255)  aaTCAAGCAACAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                                                                                                                                                  <  6:1124678/151‑1 (MQ=255)   aTCAAGGAAGAAACAAGAAAGGACGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGcc                                                                                                                                                 <  3:321035/151‑1 (MQ=255)     cAAGGAATAAACAAGAAAGGAAGTAAAGATAATTGGGGCGTGGGTGGGC‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTATCGCGTCTACCAATTCCGCCACCTTCTCACAGTCATCTTACTTTTTTTGATATCGCCTc                                                                                                                                               <  5:1211290/151‑1 (MQ=25)     cAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc                                                                                                                                               <  8:1350186/151‑1 (MQ=255)      aacgaagaaACAAGAAAGGAACTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                              <  3:1033596/148‑1 (MQ=255)       aGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGt                                                                                                                                             <  2:871232/151‑1 (MQ=255)       aGGAACAACCAACAACGGACGTAAAGATATTTCGTGCGCGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGt                                                                                                                                             <  2:1081248/151‑1 (MQ=18)         taacaaccAAGAATGGATGTAAACATATTTGGTGCGCGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGtt                                                                                                                                            <  6:582080/143‑1 (MQ=18)         gaagaaACAAGAAAGGAAGTAAAGATAATTGGGGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                           <  8:452351/151‑1 (MQ=255)         gaacaaACAAGAACGGAATTAAAGATAATGGGGGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                           <  2:307612/151‑1 (MQ=255)          aagaaACAACAAAGCAAGTAAACATAATTGGGGCGAGGGGGGGGCACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                           <  3:1093955/151‑1 (MQ=21)            aaaaCAAGAAAGGAAGTAAAGTTATTTGGTGGGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                                                                                                        <  5:720277/150‑1 (MQ=255)             aaccAAGAAAGGAAGTAAAGATAATTGGGGCGCGGGGGGCG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                       <  4:902713/148‑1 (MQ=255)             aaaaaaGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                       <  4:455076/151‑1 (MQ=255)             aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                       <  3:301490/151‑1 (MQ=255)             aaaCAAGAAAGGAAGTAAAGATAATTGGGGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                       <  7:814348/151‑1 (MQ=255)             aaaCAAGAAAGCAAGTAAACATAATTGGTGCGCGGGGGGGC‑ACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                       <  2:926808/151‑1 (MQ=255)               aaaaGAAAGGATTTATATTTAATTGGTGGCAGGGGGGGC‑ACTTCACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg                                                                                                                                     <  8:1447482/149‑1 (MQ=18)                 aaGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                                                                                                   <  1:309172/151‑1 (MQ=255)                 aaGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                                                                                                   <  7:855525/151‑1 (MQ=255)                   caaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                                                                                                 <  4:733687/150‑1 (MQ=255)                      aGGAATTAAACATAATTGGGGCCTGGGGGGGG‑ACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                                                                                                              <  7:668050/151‑1 (MQ=21)                      aGGAAGTAAAGATAATTGGGGCGCGGGGGGGG‑CCTTGAACCCCCACGGCCGTAAGGACACTACCACCTGAAGCTAGCGCGGCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                                                                                                              <  6:116431/151‑1 (MQ=25)                      aGGAAGAAAAGATAGTTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGGCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                                                                                                              <  8:994870/151‑1 (MQ=255)                      aGGAAGAAAAGATAATTGGGGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                                                                                                              <  8:741759/151‑1 (MQ=255)                                      gggggcgcGGGGGGGC‑CCTACACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc                                                                                                              <  7:371797/147‑1 (MQ=25)                                          ggcgtttgtgtggtcTTGAATCCCCACGTGCTTTATGAGACTAACACCTGCAGCTTTCTCGTCTATCAATTCCGGCACCATCGCACAGTGATCTTACTTTTTTTGATTTCTCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTc                                                                                                           <  1:832950/137‑1 (MQ=19)                                            gAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACATGAAGCTAGCGAGTCTAACAATTCAGCCAAATTCGAACAGTAATCTTACTTTTTTTGATATAGAATAGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGt                                                                                                        >  4:794210/1‑151 (MQ=21)                                              ggggtggg‑CCTTCAACCCCCACGGCCGTTATTACACTAACACCTGAAGATCCCCCGTCTACCAATTCCGCCACCTTCCCACAGTCATCTAACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTaa                                                                                                      <  3:1014161/151‑1 (MQ=25)                                                   gcgcacgcccccccccACGGCCGCAACGACACAAACCCCTGAAGCTAGCGCGGCTACCAATTCCCCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                                                                                                  <  3:791494/140‑1 (MQ=37)                                                    gggactaccaCCCCCACGGCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCCCACAGGCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                                                                                                 <  3:779284/142‑1 (MQ=37)                                                     g‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTATTTTTGATATCTCCTCGTTTGGTGCGAGGGGGGGGACTTTAACCCACACGTACGTAATAACACt                                                                                               >  6:85313/1‑151 (MQ=25)                                                      ccctccccccccccCGGCCCTACCGACCCTACCACCTGACGCTCCCGCTTCAACCAATGCCGCCACCTACCAACAGGCATCTAACTTTTTTTGATCTCCCCTCGTTTGGGGCGCGGGGGGGGACTTCAACCCCCACGTCCGTAAGAACACt                                                                                               <  8:356471/143‑1 (MQ=255)                                                       cctgccaCCCGCACGGCCATAAGGACCCGAAACCGGGACGCTCGGCCGTATAACAAATCCGCCACCAACTCACATTCAACTTACTTTTTTGGACCTACCCTACTTGGGTGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                                                                                              <  8:1382398/145‑1 (MQ=9)                                                          tcaaCCCCCACGTCCGTAAGTACAGTAAAACCTGAAGCTAGCGCTGCTAGCAATTCCGCCACCTTACCACAGTCATCAAACTTTATTTGATATCGCCTACTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                             <  5:1389648/147‑1 (MQ=21)                                                           caaCCCCCACGGCCTTTTTTACACTAACACATGATGCTTTCTATTATACCAATTCCGCCACCTAACAACAGTCATCTTACTTTTTTTGATATACCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAacac                                                                                          <  4:475865/150‑1 (MQ=16)                                                            accccccAAGTCGGTTATGACAATAACACCTGTCGCAAGCTCCGATACCAATTCCGCCAACAAAAAACATTCAACTAACATTTTATGATATAACCTACTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                                                         <  7:1365555/149‑1 (MQ=9)                                                                   ccggccttatttacacTAACACCTGAAGATACCGCTTATACAAATTCCGCCACCTACCAACAGTCATATAACTTTTTTTTATATCCCCTCTTTTGGTGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                                                                  <  5:1347163/139‑1 (MQ=18)                                                                      tCCGTATTGACACTAACTCATGAAGCTACCGAGTATACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGAAATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                 <  5:22990/149‑1 (MQ=25)                                                                             gttcacGAAAACATGAAGCTAGCACTTAAACCAATTCCGCCAGCTAACAAAACGCAACAAACTTTATTTGATTTACCCTTTTTTGGGGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                                                         <  5:744419/147‑1 (MQ=2)                                                                               acagTAAAACATGATGCTCGCGAGTAAACCAATTCAGGCACCATCCCACATTCATCTAACTTTTTTTGAAATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTAc                                                                       <  3:224070/150‑1 (MQ=17)                                                                               acacTAACACATGCTGATAGCCCGTAAACCAATTCCGCCACCAAAAAACATGCATATAACTTTATTTGATTTCCCCTTTTTTGGTGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                       <  8:604601/150‑1 (MQ=9)                                                                                                        aatctaaTTCCGCCACATAACAACATCCATCTAACTTTTTATTATATACCCTTTTTTGGGGCGTGGGGGGGGACTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                              <  7:854064/145‑1 (MQ=14)                                                                                                              aatCCGGCACCAACACACAGGCATCTAACATTAATTTATTTCACCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                       <  7:295731/149‑1 (MQ=14)                                                                                                              aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGGTTGGTGCGAGGGGGGGGACTTGCACCCCCACGTCCGTAAGAAAACTAACACCTGAAGCTAGCTCGTCTACCAATTCCGCCACCTTCGcaaactcggaagact                                       >  7:301357/1‑137 (MQ=21)                                                                                                              aTTCCGCCACAAACCAACAGTCAACAAACTTTATTTGATCTCCCCACTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                       <  8:637024/151‑1 (MQ=17)                                                                                                               ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAATAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaaacggacagctatt                                      >  3:279502/1‑136 (MQ=21)                                                                                                                  cGCCACCTTCGCACAGTCATATTACTTTTTTTGATATCGACTCGGTTGGTGCGATGCGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGATCGTCTACCAATTCCGCCACCTTCGCCCAggttgaggtcttta                                    >  8:1235522/1‑137 (MQ=18)                                                                                                                   ggcacgaaataacaGTCAACAAACTTTTTTTGATTTCCCGTTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCaa                                  <  4:1328923/140‑1 (MQ=17)                                                                                                                     caccaacaaacaGGCATCAAACTTTATTGGATATACCCTATTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACg                                <  5:1552132/142‑1 (MQ=17)                                                                                                                           atcacaGGCATTATACTTTTTTTTATTTCCCGTGTTTTGGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttt                          <  6:41246/149‑1 (MQ=18)                                                                                                                            gCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACACTGACAACAATATAAACgacgtggt                         >  3:717708/1‑144 (MQ=18)                                                                                                                            gCACAGTCAACTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACATTCGCAAACTACGAGAAATATaaaagaattttt                         >  4:323489/1‑142 (MQ=17)                                                                                                                             aacatgaaTCTAAATTTTTTTGATATAGCATATTTTGTTGCG‑TGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                       <  3:118839/144‑1 (MQ=18)                                                                                                                              acatggatattACTTTATTTGATTTCCCGTATTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                       <  8:711086/142‑1 (MQ=18)                                                                                                                              acatgcatgaaacattttttGATTTCCCCACTTTTGGTGCGAGGGGGGGGACTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                       <  6:564957/137‑1 (MQ=18)                                                                                                                                aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACACACGGAGGAATATAAAAGTGTTTTTTGGt                     >  7:630829/1‑151 (MQ=18)                                                                                                                                       aaaattttattGATATACACTATTTTGGTGCTAGGGGGGGGACTTTACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggg              <  7:417491/147‑1 (MQ=18)                                                                                                                                        aaCTTTTTTTGATCTCGCCAGATTTGGTGCGAGGGGTGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg             <  5:843932/150‑1 (MQ=255)                                                                                                                                               tttGTTATGGCCTCGTTGGGGGCTTGGGGGGGCACTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGAc       <  7:85496/150‑1 (MQ=21)                                                                                                                                                ttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGTGAGCAATATAAACGTGGTTTTTGGTGCGAgggggggtcgtt     >  7:1328144/1‑146 (MQ=255)                                                                                                                                                   atatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGAGTGAGAGCAATATAAAAGAGGTTTTTCGTCCGAGGCGGGGGACtttcc  >  3:507713/1‑148 (MQ=18)                                                                                                                                                        |                                                                                                                                                  AAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAA  >  NC_000913/4606042‑4606337

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