Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,486,580	0	A	G	45.6%	79.3 / 12.9	57	G121G (GGA→GGG)	lptF	lipopolysaccharide export ABC permease
Reads supporting (aligned to +/- strand):  ref base A (9/22);  new base G (25/1);  total (34/23)
Fisher's exact test for biased strand distribution p-value = 1.06e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAT  >  NC_000913/4486433‑4486730                                                                                                                                                    |                                                                                                                                                        gctgATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGGGGGGcc                                                                                                                                                       >  2:900367/1‑150 (MQ=255) gctgATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGAcc                                                                                                                                                       <  1:985801/150‑1 (MQ=255)     atgaCGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAACGTGATGTGGGGGGGGCCgggg                                                                                                                                                   >  2:840281/1‑150 (MQ=255)        aCGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCAt                                                                                                                                               <  1:420484/151‑1 (MQ=255)                 aaaCTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAAAGGGATGTGGGGGGGGGCGCGGGCATcgcggcct                                                                                                                                       >  1:81533/1‑146 (MQ=255)                  aaCTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGTGGTGTGGGGGGGGCCGGGGGCATcgcggccac                                                                                                                                      >  2:752908/1‑145 (MQ=255)                        ataCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGTTCGCGACCCTGGTCATCACTTCATCAGTATg                                                                                                                                <  1:228877/150‑1 (MQ=255)                        ataCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATg                                                                                                                                <  1:982064/150‑1 (MQ=255)                         taCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGGGGTTTGGGCGGGGGCGTGGGCATcgcggcctccgggt                                                                                                                                 >  2:565629/1‑138 (MQ=255)                          aCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGaa                                                                                                                              <  1:491413/150‑1 (MQ=255)                               aaGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGGGATGGGGGGGGGGCCGGGGGGCTCGCGTCGACAGGGTGAAggggt                                                                                                                         >  1:476660/1‑146 (MQ=255)                                  tGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGGTGTGGGGGGGGCCGGGGGCATCGCGGCCTCAGGGTGGAggggt                                                                                                                         >  2:354257/1‑143 (MQ=255)                                   gAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGATGTGGGGGGGGCCGGGGGCATCGCGGCCTCCGGATGAAGAGtt                                                                                                                         >  2:901632/1‑146 (MQ=255)                                     aaTTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGGTGTGGGGGGGGGCGGGGGCATCGCGTCATCAGGGTGAAatggttgga                                                                                                                     >  2:304454/1‑142 (MQ=255)                                     aaTTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTTACCTGATGTGGGGGGGGCCGTGGGCCTCGCGGCATCCGGGTGAAGTGTTTGCAGa                                                                                                                   >  2:287908/1‑150 (MQ=255)                                     aaTTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAACGTGATGTGGGGGGGACCGTGGGCGTCCCCTCATAAGGagggaagggt                                                                                                                         >  1:621869/1‑135 (MQ=255)                                        tACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAACGGGATGTGGGGGGGGCCGGGGTCATCGCGTCATCAGGGTGAAGGGTTAGCAGaa                                                                                                                  >  2:644728/1‑148 (MQ=255)                                         aCGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGGGGGGCCGCGGGCATCGCGTCCTCAGGGTGAAGGGTTTGCAGaaggg                                                                                                               >  2:392767/1‑148 (MQ=255)                                         aCGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGaagcga                                                                                                              <  2:839196/151‑1 (MQ=255)                                           ggTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGaagcgaaa                                                                                                            <  2:642435/151‑1 (MQ=255)                                            gTAATGCATGCCTGCGGCCTGAGTAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGaagcgaaat                                                                                                           <  1:222673/151‑2 (MQ=255)                                            gTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGGGGTGTGGGGGGGGGCGGGGGCATCGCGGCATCAGGGTGAAGGGTTTGCAGaa                                                                                                                  >  2:285050/1‑144 (MQ=255)                                              aaTGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGaagcgaaagc                                                                                                          <  1:213731/150‑1 (MQ=255)                                               atgcatgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGTGGTGTGGGGGGGGCCGTGGGCATCGCGGCCTCAGGGTGAAGGGGTAGCAGaagaggaa                                                                                                            >  2:258725/1‑142 (MQ=255)                                               atgcatgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGCTGTGGGGGGGGCCGCGGGCATCCCGTCCTCCGGATGAAGGGGTAGCAGaagcggaa                                                                                                            >  2:769781/1‑147 (MQ=255)                                                 gcatgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAagcgaaagcgaa                                                                                                       <  1:673342/150‑1 (MQ=255)                                                  catgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGGCCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAAc                                                                                                      <  1:930167/150‑1 (MQ=255)                                                  catgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGGTGTGGGGGGGGCCGGGGGGATCGCGGCATCAGGAGGGAGGGTTAGCAGAagcgaaagagaaac                                                                                                     >  2:424803/1‑149 (MQ=255)                                                    tgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGATGTGGGGGGGGGCGGGGGCATCGCGCCCTCCGGGGGAAGAGGTAGCAGAAGCGAAAAAGAAccc                                                                                                    >  2:863905/1‑150 (MQ=255)                                                             ccTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCCGAAGCGAAAGCGAACCCTGGCATggcg                                                                                          <  2:889978/151‑1 (MQ=255)                                                              cTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTACGGGATGTGGGGGGGGGCGGGGGCATCGCGGCCTCAGGGTGAAGGGGTTGGAGAAACGGAAAAGAACCCCGGCATggcgg                                                                                         >  2:613882/1‑151 (MQ=255)                                                                gAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGGCCGGGGGCATCGCGGCATCAGGGTGAAGGGTTTGCAGAAACGAAAGCGAAACCCGGCCTgg                                                                                            >  2:900504/1‑146 (MQ=255)                                                                 aGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGGGGGGCCGGGGGGCTCGCGGCCTCAGGATGAAGGGTTTGCAGAAGCGAAAACGAAACCCGGAATGgcggcg                                                                                       >  2:779662/1‑150 (MQ=255)                                                                     aaGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGTGGTTTGGGGGGGGCCGTGGGGATCGCGGCCTCAGGGTGGAGGGGTTGGAGAAAAGAAAACGAACACCTGCATGGGGGGGCggg                                                                                   >  2:854027/1‑150 (MQ=255)                                                                        cGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑Ca                                                                              <  2:993759/151‑1 (MQ=255)                                                                          gTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAg                                                                            <  1:119269/151‑1 (MQ=255)                                                                               ggTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGc                                                                         <  1:165881/149‑1 (MQ=255)                                                                                    aaGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGATGTGGGGGGGACCGTGGGCATCGCGGCATCAGGGTGGAGGGGTTGCAGAAACGAAAGCGAACCCCCGCATGGGGGGGCTGGCG‑GCAGGGGAAtt                                                                     >  2:308006/1‑148 (MQ=255)                                                                                     aGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGGGGGGGCGGGGGCATCGCGGCATCAGGGTGAAGGGTTAGCAGAAGCGGAAGCGAAACCTGGCATGGCGGCGCCGGCG‑CCAGGGGaaattca                                                                  >  1:942237/1‑145 (MQ=255)                                                                                      gCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGc                                                                <  1:752741/151‑1 (MQ=255)                                                                                            aTGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGATGTGGGCGGGGCCGCGGGCATCGCGGCATCAGGAGGAAGTGTTAGCAGAAGAGAAAGCGAAACCCGGCATGGGGGGGCGGGGG‑CGAGGGGAATTCCAGCAAGCg                                                           >  2:38954/1‑150 (MQ=255)                                                                                             tGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAAAGTGCTGTGGGCGGGGCCCTGGGGATCGCGGCATCAGGGTGGGGGGGTAGCAGAAGGGAAAGAGAACCCCGGCCGGGCGGGGCTGGCG‑CAAGGGGAATTTCa                                                                  >  1:286603/1‑142 (MQ=255)                                                                                                 ccTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTaa                                                      <  1:982996/150‑1 (MQ=255)                                                                                                    tGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGCGATGTGGGGGGGGCCGTGGGCATCGCGTCGTCAGGGTGAAGGGGTTGCAGAAACGAAAACGAAACCCCGCATGGGGGCGCCGGCG‑CCAGGGGGAATTCAGCCAGCGACTAAAgg                                                   >  2:58711/1‑150 (MQ=255)                                                                                                         tATTCACGGCAATCGTAGCGGCGGTTAACGTGGTTTGGGGGGGACCCCGGTGATCCCGTCATCCGGGTGAAGGGGTAGCAGAAGGGGAAACGAACACCGGCCTGGCGGGGCCGGGGCCAGGGGGAATTCC‑GCAAACGACTaa                                                      >  1:502581/1‑142 (MQ=255)                                                                                                                ggCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTg                                      <  2:441941/151‑1 (MQ=255)                                                                                                                         aGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGaa                             <  2:286550/151‑1 (MQ=255)                                                                                                                           cggcggTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAg                           <  2:354197/151‑1 (MQ=255)                                                                                                                            ggcggTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGc                          >  2:165112/1‑151 (MQ=255)                                                                                                                                gTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTg                      >  2:845390/1‑151 (MQ=255)                                                                                                                                  tAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGAc                    >  2:137456/1‑151 (MQ=255)                                                                                                                                   aaCGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGAc                    <  2:81541/150‑1 (MQ=255)                                                                                                                                   aaCGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACg                   >  1:987750/1‑151 (MQ=255)                                                                                                                                      gTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCa                <  2:829367/151‑1 (MQ=255)                                                                                                                                        gATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGc              >  1:88629/1‑151 (MQ=255)                                                                                                                                               gCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCa       <  1:565513/151‑1 (MQ=255)                                                                                                                                                 gggACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCaa      >  2:881101/1‑150 (MQ=255)                                                                                                                                                    aCCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAt  <  2:27678/151‑1 (MQ=255)                                                                                                                                                    aCCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAt  <  2:569369/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        GCTGATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCAAAGAT  >  NC_000913/4486433‑4486730

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