Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,133,245	0	A	G	41.4%	85.8 / 17.3	58	G210G (GGA→GGG)	metF	5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (8/26);  new base G (24/0);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.40e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATG  >  NC_000913/4133095‑4133395                                                                                                                                                         |                                                                                                                                                       aaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGa                                                                                                                                                        <  1:1034524/151‑1 (MQ=255)   gcgcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGGAATTA‑TTCCCGGGAt                                                                                                                                                      >  1:992706/1‑151 (MQ=255)    cgcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAAt                                                                                                                                                      <  1:1038396/150‑1 (MQ=255)     gcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAAttt                                                                                                                                                    <  1:824966/151‑1 (MQ=255)        cAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTTG‑TTTGGGAAATA‑TTTCGGGGAtttt                                                                                                                                                   >  2:913011/1‑149 (MQ=255)          ggCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCg                                                                                                                                               <  1:249030/151‑1 (MQ=255)            cGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGATTTT‑TGTGGGAATTA‑TTTCGGGGAtttt                                                                                                                                                   >  2:891911/1‑145 (MQ=255)              gATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGGAATTA‑TTTCCGGGATTTTTCCggggt                                                                                                                                           >  1:841669/1‑148 (MQ=255)              gATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTAt                                                                                                                                           <  1:908288/151‑1 (MQ=255)               aTTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGATTGA‑TTTGGAAATTATTTCCGGGAAtttttc                                                                                                                                                 >  1:11110/1‑143 (MQ=255)                   gCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGT‑TGTGGGAATTA‑TTTCGGGGATTTTTCCGGGATCTaa                                                                                                                                       >  1:468025/1‑150 (MQ=255)                    cTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGG‑TGTGGGAATTA‑TTTCCGGGATTTTTCCGGGGTCTAAc                                                                                                                                      >  2:662506/1‑150 (MQ=255)                      tAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACtt                                                                                                                                    <  2:603573/150‑1 (MQ=255)                      tAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATAGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACt                                                                                                                                     <  1:556281/149‑1 (MQ=255)                       aaTCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACttt                                                                                                                                   <  2:30635/150‑1 (MQ=255)                              aaCGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGT‑TGTGGGAAATA‑TTTCGGGAAATTTGTCGGGATCTAACTTTAAAAAg                                                                                                                             >  2:839578/1‑149 (MQ=255)                                    aaGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGT‑TGTGGGAAATA‑TTTCGGGGAATTTGCCCGTGTGTAACTTTAAACAGGCgaa                                                                                                                        >  2:213987/1‑148 (MQ=255)                                    aaGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTTCGGGAATTTTTCCGGGATCTAACTTTTAAAAGGCgaag                                                                                                                       >  2:809116/1‑149 (MQ=255)                                        ggATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTGG‑TGTGGGAATTT‑TTTCGGGGATTTTTTCCGCATCTCACTTTAAAAAGGGGAAAAAAatt                                                                                                                 >  1:544070/1‑151 (MQ=255)                                         gATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGTGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAAttt                                                                                                                 <  2:327067/150‑1 (MQ=255)                                            gCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCg                                                                                                             <  1:845559/151‑1 (MQ=255)                                             ccGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTTG‑TGGGGGAATTA‑TTTCCGGGATTTTTCCGGGCTCTTACTTTTAACAAGCGAAGAAAAATGCCg                                                                                                             >  1:449045/1‑150 (MQ=255)                                             ccGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGa                                                                                                            <  2:256285/151‑1 (MQ=255)                                               ggAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCGTTTA‑TGTGGATATTA‑TTTCCGGGATTTTGCCGGGATCTAACATTAAAAAAGCGTAGAAAatt                                                                                                                 >  2:465399/1‑144 (MQ=255)                                               ggAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGGAATTA‑TTCCGGGGAATTTGCCGGTATCTAACTTTTAAAAAGCGAAGAAAATTTcc                                                                                                              >  1:983306/1‑147 (MQ=255)                                                 aGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATg                                                                                                        <  1:645146/151‑1 (MQ=255)                                                   ccAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTGA‑TGTGGAAATTA‑TTCCGGGGAATTTTCCGGGATCTAAATTTAAAAAGGCGAAGAAAATTGCCGGTATGa                                                                                                       >  1:511839/1‑150 (MQ=255)                                                    cAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTTG‑TGTGGGAATTT‑TTTCGGGGATTTTTCCCGTATCTAACTTTAAACAAGGGGAGAAAATTTTCGATTTGAcc                                                                                                     >  2:609016/1‑151 (MQ=255)                                                       ccGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAAc                                                                                                  <  2:597305/151‑1 (MQ=255)                                                             gATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTTCGGGGATTTTTCCGGGATCTAACTTTTAAAAGGCGAAGAAATTTTGCGGTATGAACAACGTg                                                                                               >  2:350010/1‑148 (MQ=255)                                                              aTTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCGTTTATTGGGGAAATTT‑TTCGGGGAATTTT‑TCGGGGTCTAACTTTTAAAACGCGAAGAAATTTGCCGAAATGAcc                                                                                                     >  2:717530/1‑141 (MQ=255)                                                                 aCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGAATGA‑TGTGGGAATTT‑TTTCGGGGATTTTTCCGGGATTTAACTTTTAAAAGGCGGAGAAAATTGACGATATGGACAACGAGCGTAt                                                                                          >  2:870639/1‑149 (MQ=255)                                                                          ttcttcGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGGAATTA‑TTTCCGGGGTTTTTCCGGTATCTAAATTTTAACAGGCGAAGAAAATTGCCGATATGACCAACGCGCGTTTTCCGGCGTgg                                                                                >  1:758209/1‑150 (MQ=255)                                                                                   gTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCaaa                                                                      <  1:744818/151‑1 (MQ=255)                                                                                    tCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGGAATTA‑TTTTGGGGATTTTTCCGGGATCTTACTTTTAAAAGGCGAAGAAATTTGTCGATATGACCAAAGGGCGTATTTCGGCGGGGATGGGGCaaa                                                                      >  2:276118/1‑150 (MQ=255)                                                                                     cGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGATTGA‑TGTGGAAAATA‑TTTCGGGGATTTTTCCGGGATCTAACTTTTAACAGGCGAAGAAATTTTCCGATATTACCAACGTGCGTATTTCCGCGTGGGTGGCGGaaa                                                                      >  1:701644/1‑149 (MQ=255)                                                                                     cGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGGAATTA‑TTTCGGGGATTTTTCCGGGATCTAACTTTTAACAAGCGAAGAAAATTTCCGATATGACCAAAGCGCGTATTTCGGGGTGGGTGGGGCaaa                                                                      >  2:829092/1‑149 (MQ=255)                                                                                     cGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAAt                                                                     <  1:331559/150‑1 (MQ=255)                                                                                      gAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGt                                                                   <  1:654768/151‑1 (MQ=255)                                                                                       aaaGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGtt                                                                  <  1:210016/151‑1 (MQ=255)                                                                                          gCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TTTGGGAATTA‑TTCCGGGGATTTTTGCGGGATCTAACTTTTAACAGGGGGAGAAAATTTCCGAGATGAAAAACAAGCGTTTTCCCGCGTGGATGGGGCAAAAGt                                                                   >  1:377704/1‑147 (MQ=255)                                                                                           cTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGGAATTA‑TTCCGGGGATTTTTCCGGGATTTAACTTTTAACAGGCGAAGAAAATTTCCGATATGACCAACGTGCGGATTCCGGCGGGGGTGGCGCCAAAGtt                                                                  >  1:967914/1‑147 (MQ=255)                                                                                               cTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGCATGA‑TGTGGGAATTT‑TTTCGGGGATTTTTTCCGGATTTAACTTTTAAAAGGCGAAGAAAATTGTCCATATTACCAAAGTGAGTATTTCGGGGCGGGTGGGGCaaaaatttgaccg                                                            >  2:795925/1‑139 (MQ=255)                                                                                                 gCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCt                                                         <  1:967797/150‑1 (MQ=255)                                                                                                    ttttCGTGACCGCTGTGTATCGGCGGGCATTTG‑AGTGGAAATTA‑TTCCCGGGATTTTTCCGGTATCTAAATTTTAAAAAGCGGAGAAAATTGTCGATATGACCAAAACGCGTATTTCGGCGGGGGTGGCGCAAAAGTTCGACGGTCttg                                                       >  1:166384/1‑147 (MQ=255)                                                                                                        cGTGACCGCTGTGTATCGGCGGGCCTTTA‑TGTGGAAATTT‑TTTCGGGGATTTTTCCGGTATTTAACTTTTAACAGGCGAAGAAAATTGCCGATATGACCAACCTGCGTATTCCGGCGGGGGTGGCGCAAATGTTCGGCGGGCTGggtgatg                                                 >  1:175050/1‑151 (MQ=255)                                                                                                             ccGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGcc                                            <  2:968859/151‑1 (MQ=255)                                                                                                                  gtgtATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGaaa                                        <  1:606926/150‑1 (MQ=255)                                                                                                                  gtgtATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAAc                                       <  1:717628/151‑1 (MQ=255)                                                                                                                               cATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGtt                          >  2:666626/1‑151 (MQ=255)                                                                                                                                aTTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGtt                          <  2:704395/150‑1 (MQ=255)                                                                                                                                   gA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGgc                       <  2:74742/150‑1 (MQ=255)                                                                                                                                            aaTTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATAtt              <  1:680420/151‑1 (MQ=255)                                                                                                                                             atta‑ttCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTg             <  1:870752/151‑1 (MQ=255)                                                                                                                                             atta‑ttCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTg             <  1:1054869/151‑1 (MQ=255)                                                                                                                                               ta‑ttCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGcc           >  2:178016/1‑151 (MQ=255)                                                                                                                                                a‑ttCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGcc           >  1:648012/1‑150 (MQ=255)                                                                                                                                                a‑ttCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGcc           <  2:888690/150‑1 (MQ=255)                                                                                                                                                       ggAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGat     >  1:668381/1‑150 (MQ=255)                                                                                                                                                         aaTTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATg  >  1:74280/1‑151 (MQ=255)                                                                                                                                                         aaTTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATg  <  1:1101426/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                       AAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGA‑TGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATG  >  NC_000913/4133095‑4133395

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