Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,478,324	0	T	G	43.5%	105.1 / 10.3	63	G33G (GGT→GGG)	ynbB	putative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (7/28);  new base G (27/0);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 2.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACA  >  NC_000913/1478175‑1478470                                                                                                                                                        |                                                                                                                                                    gACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTg                                                                                                                                                    <  2:374341/151‑1 (MQ=255) gACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTg                                                                                                                                                    <  2:371481/151‑1 (MQ=255)     gCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAgg                                                                                                                                                 <  1:280020/150‑1 (MQ=255)        aTTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTCC‑GGAGAGAAAAGGGGGGGAgg                                                                                                                                                 >  1:77417/1‑147 (MQ=255)            aTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGGGAGAAAAAGGGGGGGGGGTACAtt                                                                                                                                          >  2:798567/1‑150 (MQ=255)            aTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑GTACC‑GGGGGGGAAAGGGGGGGGGGGTACAtt                                                                                                                                          >  2:645584/1‑150 (MQ=255)               gTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTATC‑GGAGAGAAAAGGGGGGGGGGTTACAtt                                                                                                                                          >  2:856683/1‑147 (MQ=255)               gTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGGGAGAAAAGGGGGGGAGGGTACAtt                                                                                                                                          >  2:494813/1‑147 (MQ=255)                  gTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑CGAGGGAAAAGGGGGGGAGGGGAAATTACGGA‑tt                                                                                                                                  >  1:64319/1‑151 (MQ=255)                    cGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTGCATTACGTA‑t                                                                                                                                   <  2:515048/148‑1 (MQ=255)                    cGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCg                                                                                                                                <  1:853197/151‑1 (MQ=255)                     ggAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGt                                                                                                                               <  1:169577/151‑1 (MQ=255)                     ggAGTGTGCTTCTCATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTCC‑GGGGAGAAAAGGGGGGGAGGTTATATTACGTA‑ttt                                                                                                                                 >  2:403451/1‑148 (MQ=255)                         tgtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑ATACC‑GGGGAGAAAAGGGGGGGTGGGTACAtttccgaatccttcc                                                                                                                             >  2:403326/1‑137 (MQ=255)                               tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTATC‑GGGGGGAAAAGGGGGGGGGGTTTTATTTCGTA‑TTTGTACCTgggggg                                                                                                                     >  2:838946/1‑151 (MQ=255)                                   aTGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCa                                                                                                                  <  1:263104/150‑1 (MQ=255)                                   aTGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCAt                                                                                                                 <  1:667864/151‑1 (MQ=255)                                    tGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTAC‑GGAGAGAAAAGGGGGGGAGGGTAAATTACGTA‑TTCGTACCTGGTgggg                                                                                                                    >  2:112915/1‑147 (MQ=255)                                       gcgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTCC‑GGGGGGGAAAGGGGGGGGGGGTTCAATTCGCA‑TTTGTTCCTTGTGGGGGCTCGTc                                                                                                             >  2:712407/1‑151 (MQ=255)                                       gcgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTc                                                                                                             <  2:786324/151‑1 (MQ=255)                                           tGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑CGAGAGAAAAAGGGGTGAGGGTACATTACGTA‑TTCGTACCTGGGGGGGCATCGGCAttt                                                                                                         >  1:359150/1‑151 (MQ=255)                                            gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCAttt                                                                                                         <  1:208300/150‑1 (MQ=255)                                              aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTCC‑GGGGAGAAAAGGGGGGGGGGTTAAATTACGTA‑TTTGTACCTGGTGGGGGATCGTCAtttttt                                                                                                      >  1:573619/1‑151 (MQ=255)                                                 ataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTCCTGGTGGGAAAGGGGGGGGGGGTT‑TATTTCGTA‑TTTCTTTCTGGGGGGGGATCGTCATTTGttttt                                                                                                   >  2:712462/1‑151 (MQ=255)                                                   aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGGGAAAAAGGGGGGGGGGTTCAATACGTT‑TTTGTACCTTGTGGGGGAACGTCATTTGtttttt                                                                                                  >  1:268967/1‑149 (MQ=255)                                                   aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGGAAAGGGGGGGGGGTTACATTACGTA‑TTCGTAACTGGTGGGGGATCGTCATTTGttttt                                                                                                   >  1:674201/1‑149 (MQ=255)                                                       cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGGGAGAAAAGGGGGGGGGGTTACATTAAGTA‑TTTGTACCTGGGGGGGGATCGGCATTTGTTTTTTAATgg                                                                                             >  1:38138/1‑151 (MQ=255)                                                       cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATg                                                                                              <  1:504342/150‑1 (MQ=255)                                                            aaaaTCTCTGGCAACGCTTTTCGCGTTGTTAATTGTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGAtt                                                                                        <  1:231831/151‑1 (MQ=255)                                                               aTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTAAC‑GGGGAGAAAAAGGGGGGGGGGTACATTACGGA‑TTTGTACCTGGTGGGGGATCGTCATTTTTTTTTTAATGGGGAttt                                                                                       >  1:322269/1‑149 (MQ=255)                                                                     ggCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTcc                                                                                <  1:623167/150‑1 (MQ=255)                                                                           acTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTGGC‑GGAGCGAAAAGGGGGGTTGGGTTCATTATGAA‑TTTGACACTGGTGGTGGAGCGTCATTTTTTTTTTTATTGTGATTTTAGGGTCACGCTgg                                                                         >  1:790902/1‑151 (MQ=255)                                                                            cTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACC‑GG‑GAGAAAAGGGGGGGAGGGTACATTACGGA‑TTCGTACCTGGGGGGGGATCGTCATTTGTTTTTTTATTGGGGTTTTAGGTCCACGc                                                                            >  1:296997/1‑147 (MQ=255)                                                                                      tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTAC‑GGAGAGAAAAGGGGGGGAGGGTATATTACGTA‑TTTGCACCTGGTGGGGGGACGTCATTTGTTTTTCAATGGTTATTTCAGGGCCACCCCGGGtgacggt                                                                 >  1:14473/1‑146 (MQ=255)                                                                                      tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTATC‑CGGGAGAAAAGGGGGGGGGGTTACATTACGTA‑TTCTGACCTGGTGGGGGAACGTCAATTGTTTTTTTATAGTGATTTTTGGTGCCCGCTGGGtgacgtt                                                                 >  2:580964/1‑148 (MQ=255)                                                                                      tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGGGAGAAAAGGGGGGGGGGGTTCATTAAGTA‑TTTGTTCCTGGTGGGGGGTCGTCCTTTTTTTTTTAATGGTGATTTTAGGGCCACGCTGGGtgacggt                                                                 >  1:442892/1‑146 (MQ=255)                                                                                      tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑GTTCC‑GGAGAGAAAAGGGGGGGAGGGTAACTTTCGTATTTTTTACCGGGGGGGGC‑TCGTCGTTTGTTTTTTTATGGTGAATTTAGGGTCACGGCGGAtgacg                                                                   >  1:510367/1‑146 (MQ=255)                                                                                       gttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttgac                                                              <  1:617637/150‑1 (MQ=255)                                                                                            ttttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTTC‑GGAGAGAAAAGGGGGAGGGGGTTCAGTTCATT‑TTTGTATCCGGTGGGGGCTCGTTCTTTGTTTTTTTATGGTGTTTTTAGTTCCCCCGCGGTTGACGTTGGAGttttt                                                        >  2:125820/1‑151 (MQ=255)                                                                                                     cGCTGATAAACCGTTTTCTGCTATGGCGG‑TTTCC‑GGGGAGAAAAAGGGGGGGGGGTACATTACGTA‑TTTGTACCTGGGGGGGGGTCGGCCTTTTTTTTTTTATGGGGGTTTTAGGGCCACGCTGGGTGACGGTGACGtttttt                                                       >  1:444388/1‑143 (MQ=255)                                                                                                        tGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGGGAGAAAAGGGGGGGGGGGTACATTACGTA‑TTTGTACCTGGGGGGGGATCGTCATTTGTTTTTTAATGGTGATTTTAGGGCCCCGCTGGGTGACGTTGAAGTTTTTTTTGCTGGTa                                              >  1:602258/1‑149 (MQ=255)                                                                                                         gATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAg                                            <  2:22023/150‑1 (MQ=255)                                                                                                         gATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAg                                            <  2:548885/150‑1 (MQ=255)                                                                                                            aaaCCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGttt                                         <  2:682603/150‑1 (MQ=255)                                                                                                            aaaCCGTTTTCTGCTATGGCGG‑TTACC‑AGGGGGAAAAGGGGGGGGGGGTACATTTCGTA‑TTTTTACCTGGTGGGGGAGCCTCATTTTTTTTTCCATGGGGATTTTAGGTCCACCCTGGGTTAGGGTGACGGTTTTTgtgcgggtgaggttt                                        >  1:206167/1‑140 (MQ=255)                                                                                                               ccGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTg                                     <  2:786789/151‑1 (MQ=255)                                                                                                               ccGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTg                                     <  2:73915/151‑1 (MQ=255)                                                                                                                 gTTTTCTGCTATGGCGG‑TTACC‑CGAGAGAAAAGGGGGGGGGGTTTCATTTAGGG‑TTCGTACCTGGGGGGGGATCGTCCTTTGTTTTTTAATGGTGATTTTAGGGGCACGCTGGGTGACGGTGACGGTTTTTTCGCTGTTTAGtttttt                                      >  2:244431/1‑146 (MQ=255)                                                                                                                    ttCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCAt                                 <  2:211362/150‑1 (MQ=255)                                                                                                                     tCTGCTATGGCGG‑TTACC‑GGTGAGAAAAGGGGGTGAGGGTTTCTTTCGTC‑TTTCTTACTCGGGGGGCATCGGGCTTTTTTTTTCAATTGGGGTTTCTGGTCCACGCTGGGTGGCGGTTACGTTTTTTTCGCTGGTATGtttttttgcctt                                >  2:147052/1‑142 (MQ=255)                                                                                                                        gCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGaa                             <  2:77416/150‑1 (MQ=255)                                                                                                                           aTGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAg                           <  2:98825/149‑1 (MQ=255)                                                                                                                                 g‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAatat                      >  2:661415/1‑148 (MQ=255)                                                                                                                                     aCC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACg                >  2:552175/1‑151 (MQ=255)                                                                                                                                            gagaAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtata          <  2:666951/151‑1 (MQ=255)                                                                                                                                             agaAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTatat         <  1:721512/151‑1 (MQ=255)                                                                                                                                               aaaaGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCt       >  1:323510/1‑151 (MQ=255)                                                                                                                                                aaaGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCGTCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTg      >  1:66432/1‑151 (MQ=255)                                                                                                                                                 aaGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGt     >  2:773711/1‑151 (MQ=255)                                                                                                                                                   gggggTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAc   <  1:602662/151‑1 (MQ=255)                                                                                                                                                    ggggTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACa  <  2:81335/151‑1 (MQ=255)                                                                                                                                                    ggggTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACa  <  2:444384/151‑1 (MQ=255)                                                                                                                                                    ggggTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACa  <  2:342322/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                    GACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGG‑TTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTA‑TTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACA  >  NC_000913/1478175‑1478470

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