Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,578,143	0	A	C	44.1%	89.2 / 25.0	59	L79V (TTA→GTA)	maeB	malic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (30/3);  new base C (0/26);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 9.23e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTT  >  NC_000913/2577996‑2578288                                                                                                                                                         |                                                                                                                                                  gttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAAccc                                                                                                                                                >  1:637148/1‑151 (MQ=255)    cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAAcccc                                                                                                                                               >  1:535890/1‑149 (MQ=255)         gcggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCAccg                                                                                                                                        >  1:906237/1‑151 (MQ=255)            gcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCAccgccg                                                                                                                                     <  1:18016/151‑1 (MQ=255)                        aTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTag                                                                                                                            >  2:244707/1‑148 (MQ=255)                          aattttGTCCGGGTCGAGTTCCTCATCTTCAATGTAAAATACATAAACCCCGGCAATTTTCTTAAACAGACCGCCCTTGCCTTCCACCACCGTTTT‑C‑CC‑GGCA‑GGCGCGCGATTGTT‑GCCCACCCCGACCACCGCCGTACCGTTAGAGATc                                                                                                                       <  2:953254/148‑1 (MQ=255)                                  cggggTCGAGTTCATCATCTTCAATGAAAAATCCCACACCCCCGGCAAATTTCTAAAAAAGACCCCCCTTCCCTCCCACCCCCGGTTT‑G‑CC‑TGCC‑GGCCCGCCAATGTT‑GCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGcc                                                                                                                 <  2:963137/147‑1 (MQ=255)                                  ccGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCATTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGccac                                                                                                               >  1:31900/1‑151 (MQ=255)                                       tCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGt                                                                                                          >  1:331091/1‑151 (MQ=255)                                         gttttCGTCAACTTAAATGTAAAATACATCAACCCCGGCAATTTTTTAAAACAGACCGCCCTTGCTTCCCATCCCCGTTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                         <  1:9482/147‑1 (MQ=255)                                         gtgttccacAACTTCAATGTAAAATCACCCAACCCCGCCGATTTTCTTAACCAAACCCCCCTTCCTTCCCACCCCGGGTTT‑C‑CC‑TCCC‑ACCGCCCAAAGGTT‑GCCAACCCCCGACACCGCCGTACCGTTAGAGATCACCGCCACCAgg                                                                                                           <  2:749073/140‑1 (MQ=255)                                         gatttCGTAAACTCAAATGTAAAATACATAACCCCCGGGGAATTTCTAAAACAGACCGCCCTTCCCTCCCACCCCCGGTTC‑G‑CC‑TGCC‑GGCGCGCAAGTGTT‑GCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                         <  1:753407/147‑1 (MQ=255)                                         gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                         >  1:394440/1‑150 (MQ=255)                                         gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                         >  2:582517/1‑150 (MQ=255)                                         gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                         >  2:652715/1‑150 (MQ=255)                                            ttCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTAcc                                                                                                      >  2:895020/1‑150 (MQ=255)                                                tCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCggg                                                                                                 >  1:297475/1‑151 (MQ=255)                                                tCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCggg                                                                                                 >  1:286324/1‑151 (MQ=255)                                                 cAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                                >  2:83667/1‑151 (MQ=255)                                                     tacagtGTAAAATACTACACCCCCGGCTAATTCCTCAAATAGAGTGCCTTCGCCCCCCATCCCCGTTTT‑C‑CC‑TGAC‑AGCCCGCCATTTTC‑GCCCACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCg                                                                                               <  2:714785/143‑1 (MQ=255)                                                        attgcaaaaTCCATCAACCCCGGAGATTTTCTAAACTAGAGCCCCCTTG‑CATCCAACCCCTGTTT‑G‑CC‑CCCCGAGCCCGCCATTTTT‑CCCTACCCCAGGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGta                                                                                            <  1:776880/142‑1 (MQ=255)                                                        aaTGTAAAATACATCAACCCCGGCAAATTTTTAAAACAGACCCCCCTTGCCTTCCACCCCGGTTTT‑G‑CC‑TGCC‑GGCGCGCAATTGTT‑GCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                           <  1:659273/149‑1 (MQ=255)                                                          tGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                       >  2:794678/1‑151 (MQ=255)                                                           gTCAAAAACATCAACCCCGGCGATTTTCTAAAACAGACCCCCCTCCCCTCCCATC‑CCTGTTT‑GCCC‑TGCC‑GGCGCCCCAGTGTC‑CCCAACCCCCAGCTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAtt                                                                                        <  1:386187/149‑1 (MQ=255)                                                              aaaTCCATCACCCCCGGGGATTTTCTAAAACAGACCCCCCTTCCCCTCAAGC‑GCGGTTT‑G‑CCGGGCC‑AGCGCACCAATGTC‑CCCTACCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAg                                                                                   <  2:584717/151‑1 (MQ=255)                                                              aaaTACATCATTCCGGGCGATTTTCTAAAACAGACCCCCCTTCCCCCCCCTCGCGGTTTT‑G‑CC‑TCCC‑GGCGCCCAAGTGTT‑GCCTCCCCCCCGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTa                                                                                    <  2:425091/150‑1 (MQ=255)                                                              aaaTACATCAATCCCGGCGATTTTCTTAAACAGACCGCCCTTGCCTTCCACCCCGGTTTT‑G‑CC‑GGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                       <  1:116166/147‑1 (MQ=255)                                                              aaaTACAACAGCCCCGGCAATTTTCTAAAAAAGACCCCCCTTCCCTCCCTCCACCGTTTT‑G‑CC‑GGCC‑GGCCCCCCATTGTT‑GCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTa                                                                                    <  2:659138/150‑1 (MQ=255)                                                                aTACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGc                                                                                 >  2:702916/1‑151 (MQ=255)                                                                    aTCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                             >  1:558962/1‑151 (MQ=255)                                                                       atccccGGCAATTTTTTAAAACAGACCCCCCTTCCTTCCCATCCCGTGTTT‑C‑CC‑GGCG‑GGCGC‑CCAATTTTCCCCCACCCCCAGACCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                             <  2:977068/145‑1 (MQ=255)                                                                       aaccccggaaattttCTAAAACAAACCCCCCTTCCCTCCCACCCCCGTTTT‑G‑CC‑TGCC‑AGCGCCCAAATTTT‑CCCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtt                                                                               <  1:877413/134‑1 (MQ=255)                                                                        aTCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                         >  1:949501/1‑151 (MQ=255)                                                                          ccgtgCGATTTTCTAAAACAGCCCCCCCTCCCCTCCCACCCCCGTTTT‑C‑CC‑GCCC‑AGCGCCCCATTGTT‑CCCAACCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg                                                                       <  2:592829/147‑1 (MQ=255)                                                                          cccGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg                                                                       >  1:279214/1‑151 (MQ=255)                                                                          cccGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg                                                                       >  2:422554/1‑151 (MQ=255)                                                                          cccGGCAATTTTTTAAAACAGAACGCCTTGGCCTCCCACCCCCGGTTT‑G‑CC‑TCCC‑AGCGCCCCATTGTT‑GCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                              <  1:62157/144‑1 (MQ=255)                                                                                gAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                 >  1:254324/1‑151 (MQ=255)                                                                                gAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                 >  1:373988/1‑151 (MQ=255)                                                                                gAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                 >  1:604510/1‑151 (MQ=255)                                                                                 aaTTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                                >  2:504878/1‑151 (MQ=255)                                                                                 aaTTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                                >  1:83961/1‑151 (MQ=255)                                                                                      cTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAtt                                                           >  2:159670/1‑151 (MQ=255)                                                                                        aaaaCAGACCCCCCTTCCCTCCCA‑CCCCGGTTTCG‑CC‑TGCC‑AGCGCGCCAATGTG‑GCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAtt                                                           <  1:391361/148‑1 (MQ=255)                                                                                         aaaCAGACCGCCTTTGCCTCCCACCCCCGGTTT‑G‑CC‑TGCC‑AGCGCCCCAATGTT‑GCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                            <  1:282498/147‑1 (MQ=255)                                                                                         aaaCAGACCGCCCTTGCCTCCCACCCCCGGTTT‑G‑CC‑TGCC‑AGCGCCCCAATGTT‑GCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                         <  1:866158/150‑1 (MQ=255)                                                                                              gAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                    >  2:1011450/1‑150 (MQ=255)                                                                                               aaCGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                                   >  2:261659/1‑150 (MQ=255)                                                                                                   cccTTCCCTCCCACCCCGGTTTT‑G‑CC‑GCCC‑GGCGCCCAAATTTG‑CCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                               <  1:237843/150‑1 (MQ=255)                                                                                                      ttGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                           >  1:828851/1‑151 (MQ=255)                                                                                                          ctcccaccccgggTTC‑G‑CC‑CGCC‑AGCGCGCAAGTGTT‑CCCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                           <  2:777533/136‑1 (MQ=255)                                                                                                            cccacccccGGTTT‑C‑CC‑TGCC‑GCCGCCCCAATTTG‑CCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGc                                     <  2:788771/144‑1 (MQ=255)                                                                                                              cacccccGTTTT‑G‑CC‑GGCC‑GGCCCCCCAATGTT‑GCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCt                                   <  2:925377/146‑1 (MQ=255)                                                                                                                      gTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGa                              >  2:913358/1‑148 (MQ=255)                                                                                                                                         cgcCAGTGTG‑CCCTCACCCCGCCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGa              <  2:682563/149‑1 (MQ=255)                                                                                                                                              atttttcccTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg         <  2:279042/143‑1 (MQ=255)                                                                                                                                               tGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCtg      >  1:618904/1‑151 (MQ=255)                                                                                                                                                    cccTACCCCCACCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCtgt     <  1:619423/147‑1 (MQ=255)                                                                                                                                                     ccTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGtt  <  1:203982/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                  GTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑G‑CC‑TGCC‑AGCGCGCCAATGTT‑GCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTT  >  NC_000913/2577996‑2578288

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: ‑