Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	659,940	0	T	A	44.3%	85.5 / 30.2	61	I93F (ATT→TTT)	lipA	lipoate synthase
Reads supporting (aligned to +/- strand):  ref base T (2/32);  new base A (27/0);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 2.57e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTG  >  NC_000913/659790‑660088                                                                                                                                                       |                                                                                                                                                         cgcgCAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAAt                                                                                                                                                          <  2:151481/151‑1 (MQ=255)   cgcAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAg                                                                                                                                                        <  2:483262/151‑1 (MQ=255)             caCGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGAt                                                                                                                                              <  2:1238613/151‑1 (MQ=255)               cGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCa                                                                                                                                            <  2:504643/151‑1 (MQ=255)                ggTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTGCAAAAAAAGCCGAgggac                                                                                                                                             >  1:307603/1‑146 (MQ=255)                    aaCGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAAc                                                                                                                                       <  1:318667/151‑1 (MQ=255)                       ggAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGtt                                                                                                                                    <  1:364131/151‑1 (MQ=255)                       ggAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAAAGCCGCGGAGCATaaaaatt                                                                                                                                    >  1:900977/1‑147 (MQ=255)                        gAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAgggagataaaagt                                                                                                                                     >  1:923084/1‑138 (MQ=255)                              aTAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAGGGAGATAAAAGTTggtttg                                                                                                                              >  1:567372/1‑145 (MQ=255)                                 aCCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGt                                                                                                                         <  2:299247/151‑1 (MQ=255)                                     caTAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGCACAAAAAGCCACGAGGGGGATAAACATTTCttgggc                                                                                                                            >  2:782063/1‑140 (MQ=255)                                             aGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACCAAAAAAGCCCCGAACCAAAAAAAATGCTGT‑GCCGTGGGTGGAGCaat                                                                                                             >  2:1077390/1‑149 (MQ=255)                                                 ccATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑c                                                                                                        <  2:856571/151‑1 (MQ=255)                                                  cATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑ca                                                                                                       <  2:259887/151‑2 (MQ=255)                                                       cGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGtt                                                                                                  <  2:155672/151‑1 (MQ=255)                                                        ggCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTa                                                                                                 <  2:705337/151‑1 (MQ=255)                                                        ggCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGGCAAAAAACGCCGAGGGTCATAAAAAATGCTGT‑GCCGCGGGTGGAGGATTCcc‑ccccggt                                                                                                  >  2:43439/1‑145 (MQ=255)                                                         gCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAg                                                                                                <  1:729399/151‑1 (MQ=255)                                                          cAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAATAACGCCGCGGGTCATAAAAATTGCTGT‑GGCGGGGGTGAAGGAATTCC‑CCAGGGTAgg                                                                                               >  1:132951/1‑151 (MQ=255)                                                           aaTGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAAAGCCGAGGAGCATAAAAAATGCTGT‑GCCGCGGGTGGAGAATTTCC‑CCAGGGTAggg                                                                                              >  1:934014/1‑151 (MQ=255)                                                              ggTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCag                                                                                           <  2:1030427/151‑1 (MQ=255)                                                               gTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCagg                                                                                          <  1:17112/151‑1 (MQ=255)                                                                 cTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGAAAAAAAAAGCCGCGGGGCCTAAAAGCTGCTGT‑GGCGCGGGGGAAGCCCTCCC‑CCCGGGTTGGGGagg                                                                                          >  1:364304/1‑149 (MQ=255)                                                                    cgcCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAGGGACATAAAAGGTTGTGT‑GCCGCGGGTGGAGCATTCCC‑CCCGGGTTGGGGaggggg                                                                                       >  1:616203/1‑149 (MQ=255)                                                                     gcCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGAAAAAAAGAGCCGGGGAACATAAACATTGCTGT‑GTCGGGGTTGTAGCATTTCC‑CCCGGGTAGGGCAGGGGGcg                                                                                     >  2:692181/1‑149 (MQ=255)                                                                       cAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAAAGCCGGGGATTATAAAAAGTGCTGT‑GGCGGGGGTTGAGCAATTCT‑CCAGGGTAGGGGcggagggct                                                                                    >  1:150910/1‑146 (MQ=255)                                                                       cAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAAAGCCCGGGAACATAAAACATGCTCT‑TCCCCGGGTGGAGGATTCCG‑CCCGGCTTGGGGaggggg                                                                                       >  2:846796/1‑146 (MQ=255)                                                                         gTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTc                                                                                <  1:84446/151‑1 (MQ=255)                                                                           ttCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAAAAAAGGCGAGGGATATAAAAGTTTGTGTGGCGGGGGTTGAAGATTTCCC‑CCGGGTTgggggggggggttcctc                                                                                 >  1:390319/1‑136 (MQ=255)                                                                             cACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGcag                                                                            <  1:1173790/151‑1 (MQ=255)                                                                               cTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTGCAAAAAAAGCCGGGGGGTATAAAAGATGCTGT‑GGCGGGGGTGGAGGATTCCC‑CCCGGGTTGGGGgggaggggttcccgcag                                                                            >  1:626619/1‑138 (MQ=255)                                                                                  gTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGcagacaga                                                                       <  2:1006649/151‑1 (MQ=255)                                                                                   ttCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAGGGGCATAAAAATTGCTGT‑GCCGTGGGTGAAGCAATCCG‑CCAGGGTTGGGGAGGGGGGTTCCTCGccgacaga                                                                       >  1:941356/1‑150 (MQ=255)                                                                                    tCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAAAAAGGCCCAGGGACATAAAAGATTCTGT‑GGCGGGGGTGGAGCATTTCC‑CCAGGGTTGGGGAGGGGGGGTCCCCGCCGACAGaa                                                                      >  1:1083512/1‑150 (MQ=255)                                                                                      aTTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGAGGGTCATAAAAGATGCTGT‑GCCGGGGGTGAAGCAATTCG‑CCCGGGTTGGGGAGGGGGGGTTCTCGCAGACAGaa                                                                      >  1:185354/1‑148 (MQ=255)                                                                                          gCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑c                                                             <  2:582748/151‑1 (MQ=255)                                                                                               aGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCAttt                                                        <  2:1219998/151‑1 (MQ=255)                                                                                                ggAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCAtttt                                                       <  1:712399/151‑1 (MQ=255)                                                                                                 gAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCAttt                                                        <  2:285977/149‑1 (MQ=255)                                                                                                  aGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGAACAAAAAGAGCCCCCGAGCATAAAAATTTCTTT‑GCCCCGGGTGGAGCAATTCC‑CCACGGTAGGGGAGGAGGGGTTCTCCCAGCCAAAAAGGA‑GG‑GCAattttt                                                     >  2:647887/1‑150 (MQ=255)                                                                                                   gCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTgc                                                    <  2:136134/151‑1 (MQ=255)                                                                                                       ccGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTgcgc                                                  <  2:267843/149‑1 (MQ=255)                                                                                                       ccGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCCGGGGACATAAAAAATGCTGT‑GCCGTGGGTGAAGCATTTCC‑CCCGGTTAGGGGAGGGGGGGTCCTCGCCGACAGAAAGGA‑GG‑GCAttttt                                                      >  1:1218460/1‑145 (MQ=255)                                                                                                         gggCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTg                                              <  2:654852/151‑1 (MQ=255)                                                                                                         gggCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCGGGGGGGATAAAACATGCGGT‑GCCGGGGGTGAAGGAATTCG‑CCCCGGTAGGGGAGGGGGGGTCCCCGCCGCCAGAATGGA‑GG‑CCCTTTTTTTGCAt                                                >  2:251222/1‑149 (MQ=255)                                                                                                             gACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAAAAACGCCCAGGGTCAATAAAATTTCTGT‑GCCGTGGGGGAAGGAATACC‑CCCGGGGAAGGGAGGGGGGGTCCCCGCCGACAAAAAGCA‑GG‑CCAtttttttgcctttcgg                                           >  1:799518/1‑139 (MQ=255)                                                                                                             gACCGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAACGCCCCGGGACATAAAAAATGGTGT‑GCCGGGGGTGGAGCAATCCGCCCAGGTT‑GGGGAGGGGGGGTCCCCCCCGACCAAAAGCAGGG‑CCAttttttggctttggg                                            >  2:1177049/1‑140 (MQ=255)                                                                                                              aCCGTGGGCAACATCACAGAACGGACAACGGCGGGGAACAAAAACGCCGAGGAACACAAAAACTGCTGT‑GACGTGGTGGGAGCAATCAG‑CCAGGGTAGGGGAGGGGGGGTTCTCGCCGACAGAATGCA‑GG‑GCAtttttttgcgtttcg                                            >  2:118815/1‑138 (MQ=255)                                                                                                               ccGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCt                                         <  2:1027416/150‑1 (MQ=255)                                                                                                                cGTGGGCAACATCACAGAACGGACAACGGCGGGGACAAAAAAAGCCGAGGGGCATAAAAATTGCTGT‑GCCGGGGGTGAAACATTTCC‑CCCGGGTTGGGGAGGGGGCCTTCTCCCCGACAGAAAGCA‑GG‑CCACTTTTTGGCAttttggggtt                                       >  1:890465/1‑143 (MQ=255)                                                                                                                         cATCACAGAACGGACAACGGCGGGGACAAAAAAAGCCCAGGATCATAAACAATGCTGC‑GGCGCGGGTGAAACATTCCC‑CCCCGGTTGGGGAGGGGGGGTCCTCGCCGACAGAAAGCA‑GG‑GCCTTTTTTTGCATTTCGGCGTTgtttacc                                >  2:32020/1‑143 (MQ=255)                                                                                                                           tCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTgg                             <  1:833495/150‑1 (MQ=255)                                                                                                                            cacaGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTGGAt                           <  2:338347/151‑1 (MQ=255)                                                                                                                            cacaGAACGGACAACGGCGGGTACAAAAAACGCCCGGGGTCATAAAAATTGCTGT‑GGCGTGGGTGAAGCATTTCG‑CCACGGTAGGGGAGGGGGCTTCCTCCCCGGAAGAAAGCA‑GGCCCATTTTTTC‑CATTTGGGCTTTTAttgcctcggt                           >  2:1104764/1‑142 (MQ=255)                                                                                                                                 aaCGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACgtg                      <  1:759638/151‑1 (MQ=255)                                                                                                                                     gACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTaga                  <  1:88636/151‑1 (MQ=255)                                                                                                                                         aCGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTcc              <  1:855845/151‑1 (MQ=255)                                                                                                                                                gTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGaa       <  2:648599/151‑1 (MQ=255)                                                                                                                                                  aCAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCTGCTTTGATGCCCTGGATACGGGTAGAGTCCGCTGGGAGc     >  2:678044/1‑151 (MQ=255)                                                                                                                                                   cAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCt    <  1:331783/151‑1 (MQ=255)                                                                                                                                                     aaTAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTg  <  1:517400/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                         CGCGCAGGTCATCACGGTCAACGGAGGTGATAACCACATAACGCAGCGCCATATCGGCAATGGTCTGCGCCAGTTTCACTGGTTCATTGGCATCAGGAGCTACCGGGCGACCGTGGGCAACATCACAGAACGGACAACGGCGGGTACAAATAGCGCCGAGGATCATAAACGTTGCTGT‑GCCGTGGTTGAAGCATTCCG‑CCAGGTTAGGGCAGGAGGCTTCCTCGCAGACAGAATGCA‑GG‑CCATTTTTGCGCATTGCGGCTTTGATGCCCTGGATACGTGTAGAGTCCGCTGGAAGCTTG  >  NC_000913/659790‑660088

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