Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,486,580	0	A	G	41.4%	94.8 / 11.3	58	G121G (GGA→GGG)	lptF	lipopolysaccharide export ABC permease
Reads supporting (aligned to +/- strand):  ref base A (6/28);  new base G (24/0);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 4.79e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCA  >  NC_000913/4486433‑4486725                                                                                                                                                    |                                                                                                                                                    gctgATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGAcc                                                                                                                                                   <  1:186792/150‑1 (MQ=255)    gatgaCGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTg                                                                                                                                                <  2:892525/150‑1 (MQ=255)    gatgaCGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTg                                                                                                                                                <  2:610813/150‑1 (MQ=255)       gaCGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTc                                                                                                                                             <  1:480868/150‑1 (MQ=255)       gaCGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCa                                                                                                                                            <  2:388338/151‑1 (MQ=255)         cGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATc                                                                                                                                          <  2:231870/151‑1 (MQ=255)                 aaaCTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGCGATGTGGGGGGGGCCGCGGGCATagcggcctc                                                                                                                                  >  2:351415/1‑146 (MQ=255)                   aCTGTGTACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAg                                                                                                                                <  1:838311/151‑1 (MQ=255)                       tataCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTACGTGATGGGGGCGGGGCCGTGGGCATCGCGGCATCAGGGT‑g                                                                                                                           >  1:128872/1‑151 (MQ=255)                       tataCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGATGTGGGGGGGGGCGGGGGCATcgcggcctccgggtg                                                                                                                            >  2:441004/1‑140 (MQ=255)                             gAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAgtgt                                                                                                                     <  1:587381/151‑1 (MQ=255)                             gAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGCTTACGTGATGTGGGGGGGGCCGTGGGCATCGCGTCATCCgggggtagag                                                                                                                       >  1:711375/1‑142 (MQ=255)                              aaaGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTACGTGATGTGGGGGGGGGCGGGGGGGTCGCGGCCTCCGGAT‑GAAggggt                                                                                                                    >  2:689727/1‑147 (MQ=255)                              aaaGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGtt                                                                                                                    <  2:317275/151‑1 (MQ=255)                               aaGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTACGAGAGGTGGGCGGGGGCGTGGGGAACCCGTCTTCAGGagggagaggt                                                                                                                     >  2:234733/1‑141 (MQ=255)                                aGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAg                                                                                                                  <  1:78469/151‑1 (MQ=255)                                  tGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTGAAAAGGGTGTTGGGGGGGCCGGGGGCATcgcgccgtctgg                                                                                                                               >  2:304564/1‑129 (MQ=255)                                   gAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCa                                                                                                                <  1:672616/150‑1 (MQ=255)                                    aaaTTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGGGGTGGGGGGGGGGGCGGGGGCATCGCGTCCTCAGGGT‑Gaa                                                                                                                         >  1:24114/1‑140 (MQ=255)                                      aTTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGTGGTGTGGGGGGGGACGGGGGCATCGCGGCATCAGGGT‑GGAGGGGTAGGAGaa                                                                                                             >  2:926403/1‑150 (MQ=255)                                        tACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGGGATGTGGGGGGGGGCGCGGGCATCGCGGCACCTGGGG‑GAAGGGTTTGCaaaa                                                                                                             >  2:99229/1‑148 (MQ=255)                                        tACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGGTGTGGGGGGGGCCGGGGGTCTCGCGGCATCAGGAT‑GGAggggttggagaa                                                                                                             >  2:813788/1‑137 (MQ=255)                                           ggTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGATGTGGGCGGGACCGGGGGCATCGCGGCATCAGGAT‑GGAGGGTTTGCAGaagcg                                                                                                          >  2:142860/1‑148 (MQ=255)                                             tAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGGGATGTGGGGGGGGCCGGGGGCTTTCGGTCCTCCGGGT‑GAAGGGTTATGAGaaggg                                                                                                          >  2:830170/1‑144 (MQ=255)                                             tAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGATGTGGGGGGGGCCGGGGGCATCGCGGCATCAGGGT‑GGAGGGGTTGCAGaaacgaaa                                                                                                       >  2:278071/1‑149 (MQ=255)                                                         gCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGGTGTGGGCGGGGGCCGGGGCATCGCGGCATCAGGGT‑GAAGGGGTTGGAGaagggaaaaagaaacccgg                                                                                            >  2:897424/1‑137 (MQ=255)                                                          cGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATg                                                                                        <  1:430198/151‑1 (MQ=255)                                                          cGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAACGGGGTGTGGGGGGGGGCGGGGGGATCGCGGCCTCAGGGG‑GGAGGGGTTGCAGAAGCGCAAGAGAAccc                                                                                               >  1:734452/1‑144 (MQ=255)                                                               tGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGGGATGTGGGGGGGGGCGGGGGCATCGCGGCCTCCGGAT‑GAAGGGGTAGCAGAagcgaaagcgaaaacccgcctggcgg                                                                                    >  1:406585/1‑136 (MQ=255)                                                               tGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGATGTGGGGGGGGGCCGGGGGATCGCGGCCTCAGGGT‑GGAGTGGTTGCAGAAACGGAAACGAAACCCgggctggggg                                                                                    >  2:358149/1‑142 (MQ=255)                                                                gAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATggcggc                                                                                   <  1:59891/150‑1 (MQ=255)                                                                   cAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAATGGTGTGGGGGGGGGCGTGGTCATCGCGGCCTCCAGAG‑GAAAGGGTTGCGGAAGCGAAAAAGAACCCCGGCCCGGCGggg‑cgg                                                                              >  1:285104/1‑149 (MQ=255)                                                                     aaGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGgc                                                                            <  2:251437/151‑1 (MQ=255)                                                                        cGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACATGGTGGGGGCGGGGACGGGGGGATCGCCGCCTCAGGGTGGAAGGGTT‑GGAGAAACGAAAACCAACCCCCGGCTGGCGGGG‑CTgg                                                                             >  2:2839/1‑147 (MQ=255)                                                                          gTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGTTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCa                                                                         >  2:931242/1‑149 (MQ=255)                                                                                gTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAAt                                                                 <  2:107792/151‑1 (MQ=255)                                                                                gTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGATGTGGGGGGGGCCGGGGGCATCGCGGCATCAGGAT‑GGAGTGTTAGCAGAAGCGAAAACGAACCCCGGCCTGGCGGGG‑CCgggggcaggg                                                                     >  2:80437/1‑139 (MQ=255)                                                                                  gAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTc                                                               <  2:832546/151‑1 (MQ=255)                                                                                    aaGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTATGTGGTGTGGGGGGGAGCGTGGGCATCGCGTCATCAGGGT‑GAAGTGTTTGCAGAAGGGAAAGAGAACCCCGGCGTGGCGGGG‑CCGGCGGACGGGGAAATcc                                                              >  1:98897/1‑150 (MQ=255)                                                                                     aGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGGTGTGGGGGGGGCCGGGGGCATCGCGTCATCAGGGT‑GAAGGGTTAGCAGAAGCGAAAGCGAAACCCGGCATGGGGGGG‑CCGGGGCCAAGGGAAtttc                                                              >  2:601962/1‑147 (MQ=255)                                                                                           aaTGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCg                                                      <  1:518410/151‑1 (MQ=255)                                                                                                  cTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGATGTGGGCGGGGCCGGGGTTATCGCGGCCTCAGGAT‑GAAGGGTTTGCAGAAGCGAAAACGAAACCCCGCATGGCGGGG‑CTGGCGCCAGGGGAATTTCAACAAACGACt                                                   >  2:463438/1‑147 (MQ=255)                                                                                                    tGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATgg                                              <  2:98895/150‑1 (MQ=255)                                                                                                    tGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGGTGGGGGGGGGGCCGGGGGCATCGCGGCATCAGGGT‑GGAGGGTTAGCAGAAACGAAAAAGAAACCCGGCCTGGCGGGG‑CCGGGGCCAGGGGAAATTCAACAAACGACTaaatg                                              >  2:249198/1‑147 (MQ=255)                                                                                                      cGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGATGTGGGCGGGACCGGGGGCATCGCGTCCTCAGGGT‑GGAGGGGTTGGAGAAGCGAAAACGAAACCCGGCCTGGCGGGGCCGGGGGCAGGGGGAATT‑CAGCCAGCGACCAATgg                                              >  1:951033/1‑148 (MQ=255)                                                                                                          aTTCACGGCAATCGTAGCGGCGGTTTAAGTGATGTGGGGGGGGCCCTGGGCATCGCGGCCTCAGGGT‑GAAGGGGTTGGAGAAGGGAAAAAGAAACCCGGCCTGGCGGCG‑CCGGCGCCAGGGGAATTTCCGCAAGCGACAAATgggcgccc                                        >  1:840691/1‑144 (MQ=255)                                                                                                            tCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCgg                                      <  1:805588/150‑1 (MQ=255)                                                                                                                ggCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTg                                 <  1:926386/151‑1 (MQ=255)                                                                                                                          gcggcgGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCg                          <  1:705844/148‑1 (MQ=255)                                                                                                                               ggTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGtt                  <  1:443724/151‑1 (MQ=255)                                                                                                                                gTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTg                 <  1:181407/151‑1 (MQ=255)                                                                                                                                     cGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGc            >  2:861494/1‑151 (MQ=255)                                                                                                                                     cGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGc            <  2:513148/151‑1 (MQ=255)                                                                                                                                     cGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGc            <  1:60756/151‑1 (MQ=255)                                                                                                                                     cGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGc            >  2:19516/1‑151 (MQ=255)                                                                                                                                           gtgGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGa       <  2:681655/150‑1 (MQ=255)                                                                                                                                            tgGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGAt      <  2:759335/150‑1 (MQ=255)                                                                                                                                               gCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCa  <  1:440994/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCTGATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGAT‑GAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCG‑CTGGCGCAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCA  >  NC_000913/4486433‑4486725

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