Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,424	0	C	A	41.9%	70.2 / 38.1	62	V37F (GTT→TTT)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base C (26/10);  new base A (0/26);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.30e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTG  >  NC_000913/2248274‑2248570                                                                                                                                                        |                                                                                                                                                   gACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGcaac                                                                                                                                                    <  2:340448/151‑1 (MQ=255) gACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGcaac                                                                                                                                                    <  1:297106/151‑1 (MQ=255)    cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGcaacaa                                                                                                                                                  >  2:533016/1‑150 (MQ=255)     gTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACaaa                                                                                                                                                 >  2:826969/1‑150 (MQ=255)      tAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACaaaa                                                                                                                                                >  1:696135/1‑150 (MQ=255)         gTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAc                                                                                                                                            <  2:340853/151‑1 (MQ=255)           ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAcc                                                                                                                                           <  2:857664/150‑1 (MQ=255)            gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTg                                                                                                                                         <  1:470603/151‑1 (MQ=255)              aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGt                                                                                                                                        >  2:976850/1‑150 (MQ=255)                aCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTAcc                                                                                                                                     <  1:837980/151‑1 (MQ=255)                           aCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGc                                                                                                                          >  2:948440/1‑151 (MQ=255)                                cATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGAACCGATGGAACCGCCAATg                                                                                                                     >  1:124848/1‑151 (MQ=255)                                  tatatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGc                                                                                                                   >  2:787008/1‑151 (MQ=255)                                   atatGCAGCCAGTCCCCCGAGACGGGTCATCAGGAAGAAACCCACCAGGCCAATCAGCATATACGAGAGAAACGCCCCCCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAAAAAAAAGACCTGTACCGATGGAACCGCCAATGGc                                                                                                                   <  1:807283/150‑1 (MQ=255)                                    tatGCAGCAAGTTCACCGAGACTGGCCATCAGGAAGTAAAACATCCGGCCAATCAGCATATCCGAGAGCAATGCCCCGCGCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCa                                                                                                                  >  1:284851/1‑150 (MQ=255)                                          ccaattcccccGAGACTGGTCACCAGAAAGAAACCCAACAGCCAAACAAGCATACCCGAGAGCAATCCCCCCCCCGGC‑CCGGCCGGAGAAATCGTGGCGCCAGAGAAAAAAAAAAACCTTTTACCGATGAAACCGCCAATGGCAatcatc                                                                                                            <  1:951085/141‑1 (MQ=255)                                           ccAGTCCCCCGAGATGGGCCACCAGAAAGAAACCCATCAGCCAAATCAGCATAAACGAGAGAAATCCCCCGCCCGGC‑CCTGCCTGAAAAATCTTGGCGCCAGAGGCAAAAAAAAGACCTGTACCGATGGAACCGCCAATGGCAatcatc                                                                                                            <  1:18208/149‑1 (MQ=255)                                               tccACCGAGACTGGTCATCAGAAATAAACCCATCAGCCCAATAAGAATATACGAAACCAAGCCCCCGCCCGGA‑CCTGCCGGAGAAATCGTTCCGCCAGAGACAAAAAAAAGCCCTTTACCGATGGAACCGCCAATGGCAATCATCGTc                                                                                                         <  1:466603/146‑1 (MQ=255)                                                  aCCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTg                                                                                                    >  1:127820/1‑150 (MQ=255)                                                   ccGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAc                                                                                                  >  1:695334/1‑151 (MQ=255)                                                     gagaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                                >  1:376534/1‑151 (MQ=255)                                                       gaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACgcgc                                                                                              >  1:437862/1‑151 (MQ=255)                                                        acgggTCATCAGAAAGAAACCCATCGGCCAAACCAGCATATCCGAGAGCAAGCCCCCCCCCGGC‑CCTGCCTGAAAAATCGTGGCGCCGGAGCAAAAAAAAAGCCCTGTCCCGTGGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                                <  2:309395/145‑1 (MQ=255)                                                            gTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                          >  2:1113354/1‑150 (MQ=255)                                                                  aGGAAGTAACCCATCAGCCCAAGCAGAATATACGAGACAAACCCCCCCCCCGCC‑CCTGCCGGAAAAATCGTGGCGCAAGAGGCAAAAAAAAGCCCTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTc                                                                                     <  1:348666/149‑1 (MQ=255)                                                                   ggAAGAAACCCACCAGCCAAAACAGAATAAACGAAAGAAAT‑CCCCCCCCGGACCCTGCCGGAAAAATCGTTGCGCCGGGGAAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAtt                                                                                      <  2:752868/147‑1 (MQ=255)                                                                     aaGAAACCCATCAGACCAATCACCACAGACGAGAGCAAGGCCCCGCCCGGA‑CCGCCCGAAGAAACCGTGGCGCAAGGGGAAACAAAAAGCCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACgg                                                                                 <  1:337172/150‑1 (MQ=255)                                                                      aGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCg                                                                               >  1:773592/1‑151 (MQ=255)                                                                          acccATCAGGCCAATCAGCATATACGAGAGCAAGGCCCCGCCCGGA‑CCTGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCGGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                             <  2:888904/147‑1 (MQ=255)                                                                            ccATCAGGCCAATCAGCATATACGAGAGCAATGCTCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGcc                                                                         >  1:191652/1‑151 (MQ=255)                                                                              aTCAGCCAAACCAGCATATACGAGAGCAATGCCCCCCCCGGA‑CCTGCCTGAAAAATCGTTGCCCCAGAGGCAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGccc                                                                        <  2:622642/150‑1 (MQ=255)                                                                                 aGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgcg                                                                    >  2:1015061/1‑151 (MQ=255)                                                                                  gcccAATCAGCATATACGAGAGCAATCCCCCCCCCGGC‑CCTCCCTGAAAATTCTTTGCGCCAGAGCCAACAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCgg                                                                      <  2:263028/146‑1 (MQ=255)                                                                                      aacaaGCATATACGAGAGCACCCCCCCGCCCGGC‑CCTGCCTGAAAAACCGTGGCGCCAGAGCAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCt                                                                  <  1:109555/144‑1 (MQ=255)                                                                                      aaTCAAAATAAACGAGAGCAATCCCCCCCCGGCC‑CCTGCCTGAGAAATCGTGGCCCCAGAGCCAACAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCtt                                                                 <  2:729982/149‑1 (MQ=255)                                                                                       aTCAGAAAAAACGAGAGCAATCCCCCGCCGGGC‑CCTGCTGAAAAAATCGTGGCGCCGGAGGCAAAAAAAAGCCCGTTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTc                                                                <  1:100395/149‑1 (MQ=255)                                                                                         cAGAATATACGAGAGCAATCCCCCGCCCGGA‑CCTGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                             <  2:672885/150‑1 (MQ=255)                                                                                          agaatatAGGAGAGCAAGCCCCCGCCCGGC‑CCGGCCTAAAAAATCGTTCCGCCAGAGCCAAAAAAAAGCCCGTTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                             <  2:582090/146‑1 (MQ=255)                                                                                          aacATATACGAGAGCAATGCCCCGCCCGGC‑CCGGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTgg                                                           <  2:887041/149‑1 (MQ=255)                                                                                          aacATAAAGGAGAGCAACCCCCCGCCCGGC‑CCTGCCGGAAAAATCTTGGCCCCAGAGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtttg                                                            <  1:532939/148‑1 (MQ=255)                                                                                          aGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTgg                                                           >  1:843966/1‑151 (MQ=255)                                                                                            cATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                         >  1:536387/1‑151 (MQ=255)                                                                                             atatCCGAGACCAATGCCCCGCCCGGC‑CCTGCCAGAAAAACCGTCGCGCCAGAGGCAACAAAAAGCCCTGCACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTgg                                                           <  1:256631/148‑1 (MQ=255)                                                                                               aaccgaaaccAACGCCCCGCCCGGC‑CCTGCCGGAAAAATCGTGGCGCCAGGGGCAAAAAAAACCCCGGTACCAATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTa                                                      <  2:190224/142‑1 (MQ=255)                                                                                               aaaCGAGAGCAATCCCCCCCCCGGC‑CCTGCCTGAAAAATCGTTGCGCCAGAGGCAACAAAAAGCCCTGTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGt                                                          <  2:560633/145‑1 (MQ=255)                                                                                                 aCGAGAGCATTGCCCCGCCCGCC‑CCTGCCTGAAAAATCGTTGCGCCAGAGGCAACAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTa                                                      <  2:496316/149‑1 (MQ=255)                                                                                                 aCGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                                    >  1:7648/1‑151 (MQ=255)                                                                                                      tgtaaTTCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCACTCATCGTCAGGTGACGCGCCTTTAATTGACGGCGTAAGCCCGGCGCTTCTGTGGTTTAAGttt                                                  >  1:244526/4‑148 (MQ=255)                                                                                                          atcccccGCCCGGC‑CCTGCCGGAGAAAGCGTTGCGCCAGAGGCAAAAAAAACCCCTGTCCCTATGGACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                             <  1:950216/146‑1 (MQ=255)                                                                                                             ccccGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAACAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGCGACGCGCCTTTAATTCACGGCGTAAGCCCGACGCTTCTGTGGTTTTAGTTTCGGAAAcaa                                         >  1:414896/1‑148 (MQ=255)                                                                                                                cgcccgGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATAc                                      >  1:229832/1‑150 (MQ=255)                                                                                                                 gcccgGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACa                                     >  2:433920/1‑150 (MQ=255)                                                                                                                   cccgccccgccgtaagaAAGCGTTGCGCCAGAGAAAAAAAAAACACCTTCCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaa                                    <  1:150842/136‑1 (MQ=255)                                                                                                                      ggccctccctaAAAAATCGTTGCGCCAGAGCCAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaaaa                                 <  2:681871/147‑1 (MQ=255)                                                                                                                      gA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                               >  1:390377/1‑151 (MQ=255)                                                                                                                          ctgcctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCt                             >  2:632846/1‑150 (MQ=255)                                                                                                                               tgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                       >  1:819696/1‑151 (MQ=255)                                                                                                                                 aaaaaTCGTGGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                       <  2:1020809/147‑1 (MQ=255)                                                                                                                                   aaaTCGTTGCGCCGGGGGAAAAAAAAAGCCCTGTCCCAAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                       <  2:187591/147‑1 (MQ=255)                                                                                                                                   aaaTCGTGGCGCAAGAGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa                   <  1:663462/151‑1 (MQ=255)                                                                                                                                    aaTCTTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaa                  <  2:481664/151‑1 (MQ=255)                                                                                                                                     aTCGTGGCGCCAGAGGCAAAAAAAAGCCCGGTACGGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaa                 <  1:942910/151‑1 (MQ=255)                                                                                                                                            cgcAAGAGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAg           <  2:814833/150‑1 (MQ=255)                                                                                                                                                 gagGCAAAAAAAAGCCCTGCCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcgc       <  2:757227/149‑1 (MQ=255)                                                                                                                                                    gCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTg  >  2:898222/1‑151 (MQ=255)                                                                                                                                                      aacaaAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGAtt   <  2:446833/148‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   GACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTG  >  NC_000913/2248274‑2248570

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