Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,431	0	A	C	48.5%	66.4 / 30.2	67	G34G (GGT→GGG)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (27/7);  new base C (0/32);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 2.28e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACG  >  NC_000913/2248282‑2248579                                                                                                                                                       |                                                                                                                                                      gTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAc                                                                                                                                                      <  2:340853/151‑1 (MQ=255)   ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAcc                                                                                                                                                     <  2:857664/150‑1 (MQ=255)    gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑Tg                                                                                                                                                  <  1:470603/151‑1 (MQ=255)      aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGt                                                                                                                                                 >  2:976850/1‑150 (MQ=255)        aCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTAcc                                                                                                                                              <  1:837980/151‑1 (MQ=255)                   aCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGc                                                                                                                                   >  2:948440/1‑151 (MQ=255)                        cATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGAACCGATGGAACCGCCAATg                                                                                                                              >  1:124848/1‑151 (MQ=255)                          tatatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGc                                                                                                                            >  2:787008/1‑151 (MQ=255)                           atatGCAGCCAGTCCCCCGAGACGGGTCATCAGGAAGAAACCCACCAGGCCAATCAGCATATACGAGAGAAACGCCCCCCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAAAAAAAAGACC‑TGTACCGATGGAACCGCCAATGGc                                                                                                                            <  1:807283/150‑1 (MQ=255)                            tatGCAGCAAGTTCACCGAGACTGGCCATCAGGAAGTAAAACATCCGGCCAATCAGCATATCCGAGAGCAATGCCCCGCGCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCa                                                                                                                           >  1:284851/1‑150 (MQ=255)                                  ccaattcccccGAGACTGGTCACCAGAAAGAAACCCAACAGCCAAACAAGCATACCCGAGAGCAATCCCCCCCCCGGC‑CCGGCCGGAGAAATCGTGGCGCCAGAGAAAAAAAAAAACCT‑TTTACCGATGAAACCGCCAATGGCAatcatc                                                                                                                     <  1:951085/141‑1 (MQ=255)                                   ccAGTCCCCCGAGATGGGCCACCAGAAAGAAACCCATCAGCCAAATCAGCATAAACGAGAGAAATCCCCCGCCCGGC‑CCTGCCTGAAAAATCTTGGCGCCAGAGGCAAAAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAatcatc                                                                                                                     <  1:18208/149‑1 (MQ=255)                                       tccACCGAGACTGGTCATCAGAAATAAACCCATCAGCCCAATAAGAATATACGAAACCAAGCCCCCGCCCGGA‑CCTGCCGGAGAAATCGTTCCGCCAGAGACAAAAAAAAGCCC‑TTTACCGATGGAACCGCCAATGGCAATCATCGTc                                                                                                                  <  1:466603/146‑1 (MQ=255)                                          aCCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTg                                                                                                             >  1:127820/1‑150 (MQ=255)                                           ccGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAc                                                                                                           >  1:695334/1‑151 (MQ=255)                                             gagaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                                         >  1:376534/1‑151 (MQ=255)                                               gaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACgcgc                                                                                                       >  1:437862/1‑151 (MQ=255)                                                acgggTCATCAGAAAGAAACCCATCGGCCAAACCAGCATATCCGAGAGCAAGCCCCCCCCCGGC‑CCTGCCTGAAAAATCGTGGCGCCGGAGCAAAAAAAAAGCCC‑TGTCCCGTGGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                                         <  2:309395/145‑1 (MQ=255)                                                    gTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                                   >  2:1113354/1‑150 (MQ=255)                                                          aGGAAGTAACCCATCAGCCCAAGCAGAATATACGAGACAAACCCCCCCCCCGCC‑CCTGCCGGAAAAATCGTGGCGCAAGAGGCAAAAAAAAGCCC‑TGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTc                                                                                              <  1:348666/149‑1 (MQ=255)                                                           ggAAGAAACCCACCAGCCAAAACAGAATAAACGAAAGAAAT‑CCCCCCCCGGACCCTGCCGGAAAAATCGTTGCGCCGGGGAAAAAAAAAAGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAtt                                                                                               <  2:752868/147‑1 (MQ=255)                                                             aaGAAACCCATCAGACCAATCACCACAGACGAGAGCAAGGCCCCGCCCGGA‑CCGCCCGAAGAAACCGTGGCGCAAGGGGAAACAAAAAGCCC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACgg                                                                                          <  1:337172/150‑1 (MQ=255)                                                              aGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCg                                                                                        >  1:773592/1‑151 (MQ=255)                                                                  acccATCAGGCCAATCAGCATATACGAGAGCAAGGCCCCGCCCGGA‑CCTGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCC‑TGTCCGGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                                      <  2:888904/147‑1 (MQ=255)                                                                    ccATCAGGCCAATCAGCATATACGAGAGCAATGCTCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGcc                                                                                  >  1:191652/1‑151 (MQ=255)                                                                      aTCAGCCAAACCAGCATATACGAGAGCAATGCCCCCCCCGGA‑CCTGCCTGAAAAATCGTTGCCCCAGAGGCAAAAAAAAGCCC‑TGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGccc                                                                                 <  2:622642/150‑1 (MQ=255)                                                                         aGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgcg                                                                             >  2:1015061/1‑151 (MQ=255)                                                                          gcccAATCAGCATATACGAGAGCAATCCCCCCCCCGGC‑CCTCCCTGAAAATTCTTTGCGCCAGAGCCAACAAAAAGCCC‑TGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCgg                                                                               <  2:263028/146‑1 (MQ=255)                                                                              aacaaGCATATACGAGAGCACCCCCCCGCCCGGC‑CCTGCCTGAAAAACCGTGGCGCCAGAGCAAAAAAAAAGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCt                                                                           <  1:109555/144‑1 (MQ=255)                                                                              aaTCAAAATAAACGAGAGCAATCCCCCCCCGGCC‑CCTGCCTGAGAAATCGTGGCCCCAGAGCCAACAAAAAGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCtt                                                                          <  2:729982/149‑1 (MQ=255)                                                                               aTCAGAAAAAACGAGAGCAATCCCCCGCCGGGC‑CCTGCTGAAAAAATCGTGGCGCCGGAGGCAAAAAAAAGCCC‑GTTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTc                                                                         <  1:100395/149‑1 (MQ=255)                                                                                 cAGAATATACGAGAGCAATCCCCCGCCCGGA‑CCTGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCC‑TGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                      <  2:672885/150‑1 (MQ=255)                                                                                  agaatatAGGAGAGCAAGCCCCCGCCCGGC‑CCGGCCTAAAAAATCGTTCCGCCAGAGCCAAAAAAAAGCCC‑GTTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                      <  2:582090/146‑1 (MQ=255)                                                                                  aacATATACGAGAGCAATGCCCCGCCCGGC‑CCGGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCC‑TGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTgg                                                                    <  2:887041/149‑1 (MQ=255)                                                                                  aacATAAAGGAGAGCAACCCCCCGCCCGGC‑CCTGCCGGAAAAATCTTGGCCCCAGAGGCAAAAAAAAGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtttg                                                                     <  1:532939/148‑1 (MQ=255)                                                                                  aGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTgg                                                                    >  1:843966/1‑151 (MQ=255)                                                                                    cATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                                  >  1:536387/1‑151 (MQ=255)                                                                                     atatCCGAGACCAATGCCCCGCCCGGC‑CCTGCCAGAAAAACCGTCGCGCCAGAGGCAACAAAAAGCCC‑TGCACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTgg                                                                    <  1:256631/148‑1 (MQ=255)                                                                                       aaccgaaaccAACGCCCCGCCCGGC‑CCTGCCGGAAAAATCGTGGCGCCAGGGGCAAAAAAAACCCC‑GGTACCAATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTa                                                               <  2:190224/142‑1 (MQ=255)                                                                                       aaaCGAGAGCAATCCCCCCCCCGGC‑CCTGCCTGAAAAATCGTTGCGCCAGAGGCAACAAAAAGCCC‑TGTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGt                                                                   <  2:560633/145‑1 (MQ=255)                                                                                         aCGAGAGCATTGCCCCGCCCGCC‑CCTGCCTGAAAAATCGTTGCGCCAGAGGCAACAAAAAGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTa                                                               <  2:496316/149‑1 (MQ=255)                                                                                         aCGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                                             >  1:7648/1‑151 (MQ=255)                                                                                              tgtaaTTCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCACTCATCGTCAGGTGACGCGCCTTTAATTGACGGCGTAAGCCCGGCGCTTCTGTGGTTTAAGttt                                                           >  1:244526/4‑148 (MQ=255)                                                                                                  atcccccGCCCGGC‑CCTGCCGGAGAAAGCGTTGCGCCAGAGGCAAAAAAAACCCC‑TGTCCCTATGGACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                      <  1:950216/146‑1 (MQ=255)                                                                                                     ccccGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAACAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGCGACGCGCCTTTAATTCACGGCGTAAGCCCGACGCTTCTGTGGTTTTAGTTTCGGAAAcaa                                                  >  1:414896/1‑148 (MQ=255)                                                                                                        cgcccgGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATAc                                               >  1:229832/1‑150 (MQ=255)                                                                                                         gcccgGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACa                                              >  2:433920/1‑150 (MQ=255)                                                                                                           cccgccccgccgtaagaAAGCGTTGCGCCAGAGAAAAAAAAAACACC‑TTCCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaa                                             <  1:150842/136‑1 (MQ=255)                                                                                                              ggccctccctaAAAAATCGTTGCGCCAGAGCCAAAAAAAAGCCC‑TGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaaaa                                          <  2:681871/147‑1 (MQ=255)                                                                                                              gA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                        >  1:390377/1‑151 (MQ=255)                                                                                                                  ctgcctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCt                                      >  2:632846/1‑150 (MQ=255)                                                                                                                       tgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                                >  1:819696/1‑151 (MQ=255)                                                                                                                         aaaaaTCGTGGCGCCAGAGGCAAAAAAAAGCCC‑TGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                                <  2:1020809/147‑1 (MQ=255)                                                                                                                           aaaTCGTTGCGCCGGGGGAAAAAAAAAGCCC‑TGTCCCAAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                                <  2:187591/147‑1 (MQ=255)                                                                                                                           aaaTCGTGGCGCAAGAGGCAAAAAAAAGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa                            <  1:663462/151‑1 (MQ=255)                                                                                                                            aaTCTTTGCGCCAGAGGCAAAAAAAAGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaa                           <  2:481664/151‑1 (MQ=255)                                                                                                                             aTCGTGGCGCCAGAGGCAAAAAAAAGCCC‑GGTACGGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaa                          <  1:942910/151‑1 (MQ=255)                                                                                                                                    cgcAAGAGGCAAAAAAAAGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAg                    <  2:814833/150‑1 (MQ=255)                                                                                                                                         gagGCAAAAAAAAGCCC‑TGCCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcgc                <  2:757227/149‑1 (MQ=255)                                                                                                                                          gagaaaaaaaaaGGCCTTGTACGAATGGATCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAgg                   <  1:146299/139‑1 (MQ=255)                                                                                                                                            gCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTg           >  2:898222/1‑151 (MQ=255)                                                                                                                                              aacaaAAAGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGAtt            <  2:446833/148‑1 (MQ=255)                                                                                                                                              aaaaaaaaGCCT‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGAtt            <  1:463645/145‑1 (MQ=255)                                                                                                                                               aaaaaaaGCCC‑TGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAg        <  1:948313/149‑1 (MQ=255)                                                                                                                                                  aaaaGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGc       >  1:928931/1‑149 (MQ=255)                                                                                                                                                     aGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACg  >  2:583986/1‑151 (MQ=255)                                                                                                                                                     aGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACg  >  1:709788/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      GTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGA‑CCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACC‑TGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACG  >  NC_000913/2248282‑2248579

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