Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,397,830	0	T	C	45.9%	67.3 / 22.6	61	I265V (ATT→GTT)	uspE	stress‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (23/10);  new base C (0/28);  total (23/38)
Fisher's exact test for biased strand distribution p-value = 2.80e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACGC  >  NC_000913/1397684‑1397979                                                                                                                                                      |                                                                                                                                                      ttcGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑Tcccc                                                                                                                                                   >  2:771259/1‑151 (MQ=255) ttcGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGCGCCCAGAACCACAA‑Tcccc                                                                                                                                                   >  2:780992/1‑151 (MQ=255)   cGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCg                                                                                                                                                  <  2:433515/150‑1 (MQ=255)    gTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGcc                                                                                                                                                <  2:60689/151‑1 (MQ=255)      cATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTg                                                                                                                                              <  2:560672/151‑1 (MQ=255)               cacggggttttggATACTGTTCAGGTTAAAGGACCAGCAGGTCGCAGCAAAGATGTTCAACCACCTGTCCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAAGCCCGGTCCCTCCCCCCGTGCCCACACCCCCAC‑CCCCCGCCTGTAAGTGCTc                                                                                                                                     <  2:476793/139‑1 (MQ=255)                aaCGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGc                                                                                                                                                 >  1:704814/1‑138 (MQ=255)                aaCGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGc                                                                                                                                                 <  2:704765/138‑1 (MQ=255)                   gggTTTCTTATACTGGTCAGTTTTAATGCCCAGCAGGTCGCAGCGAAGATGAACAATCCCCTGTCCCGCC‑GTGTTGCCGAGGAAGGCTGCTG‑AAATCCCGGTGCGTCCTACCGGCCCCAAACCCACAA‑CCCCCGCCTGTAAGTGCTCCGcc                                                                                                                                 <  1:712075/151‑1 (MQ=255)                    ggTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGcc                                                                                                                                 >  2:885167/1‑150 (MQ=255)                    ggTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGcc                                                                                                                                 >  2:485811/1‑150 (MQ=255)                            aTACTGGTCAGGTTTAATGCCCAGCAGGTCGCAGCGAAGATGATCAATCCCCTGTCCCGCC‑GTTTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAg                                                                                                                          <  1:478904/149‑1 (MQ=255)                                 ggTCAGGTTAAATGACCAGCAGGTCGCAGCAAAGATGACCAATCCCCTGTCCCGCC‑GTGTTGCCGAGGAATGCTGCGG‑AAATCCCGGTGCGCCCTACCGTGCCAAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATc                                                                                                                     <  2:237731/149‑1 (MQ=255)                                      ggTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCtctt                                                                                                              >  2:648406/1‑151 (MQ=255)                                       gTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACctct                                                                                                               >  1:8657/1‑149 (MQ=255)                                        tttaagacccAGCAGGCCGCAGAAAAAATGGTCACTCCCTTTTCCCGCC‑GTGTTGCAGAGAAATGCGGCTG‑AATTCCCGGACCGTCCTACCGGCCCCAAAACCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACct                                                                                                                 <  2:16158/138‑1 (MQ=255)                                           aaTGCCCAGCAGGTCGCAGCGAAGATGATAAATCCCCTGTCCCGCC‑GTTTGGCCGAGGAATGCTGCTA‑AAATCCCGGGGCGCCCTACCGTGCCCAGAACCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTg                                                                                                           <  2:118962/149‑1 (MQ=255)                                           aaTGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGc                                                                                                         >  2:305556/1‑151 (MQ=255)                                             tGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGc                                                                                                         >  1:890675/1‑149 (MQ=255)                                                      gttCGCACCGAAAGTGAACAATCCCCTGTCCCGCC‑CTGTTCCCGAGGAATCCGGCAA‑AAATCCCGGTGCGCCCTCCCGTGCCCAAAACCCAAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCtttt                                                                                                <  1:567156/147‑1 (MQ=255)                                                      ggTCGCAGCGAGGATGATCAACCCCCTG‑TCCGCCGGTGTGGCCGAGGAATGCGGCTG‑AAATCCGGGTGGGCCCTCCCGGCCCCAGCCCCACAC‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTc                                                                                              <  1:157858/151‑1 (MQ=255)                                                      ggTCGCAGCGAAGATGATCAATCCCCTGTCCCGCC‑GTGTGGCCGAGAAATGCTGCTG‑AAATCCCGGTGCGTCCTACCGTGCCCAGAACCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCttt                                                                                                 <  2:423532/148‑1 (MQ=255)                                                       gTCGCAGCGAAGATGATCAATCCCCTGTCCCGCC‑GTGTTGCGGAGGAATGCTGCTA‑AAATCCCGGTGCGCCCTCCCGTGCCCAGAACCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCt                                                                                             <  2:575116/151‑1 (MQ=255)                                                         cggaGCGAAGAGAATCAATCACCTGTTCCGCC‑TTTTTGCCGAGGAATGCTGCTG‑AAATGCCGGTGCTTCCAACCGTCCCCAAACCCACAC‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCt                                                                                             <  1:905261/146‑1 (MQ=255)                                                           cAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGt                                                                                         >  1:50010/1‑151 (MQ=255)                                                            aGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTg                                                                                        >  2:292391/1‑151 (MQ=255)                                                            aGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTg                                                                                        >  2:186817/1‑151 (MQ=255)                                                                     gatCAACCCCCTGTCCCGCC‑GTGTTGCCGAGGAATGCTGCTA‑AAATCCCGGTGCGTCCTACCGTGCCCAGACCCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCAt                                                                                  <  2:12409/148‑1 (MQ=255)                                                                       tcaatcaCCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGTCGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGttt                                                                              >  2:228906/1‑150 (MQ=255)                                                                         aatccCCTGTCCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATCCGGGTGCGCCCTCCCGTCCCAAGACCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGt                                                                                <  2:40851/146‑1 (MQ=255)                                                                         aatcaCTTTTTCCGCC‑TTGTTGCGCAGTAATCCTGTCA‑AAATGCGGTTGCTCCCTACCGTCCCTATCACCACAC‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCAt                                                                                  <  2:108174/144‑1 (MQ=255)                                                                         aatcaCCTGTCCCGCC‑GTTTTGCCGAGAAATCCTGCTA‑AAATCCCGGTGCTCCCTACCGTGCCCAGAACCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGttt                                                                              <  1:610166/148‑1 (MQ=255)                                                                           tcaCCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCAtt                                                                         >  1:107502/1‑151 (MQ=255)                                                                            ccccTGTCCCGCC‑GTGTGGCCGAGCAATGCTGCTA‑AAAGACGGGGGCTTCCCCCGGTGCCAAAAACCACAC‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCa                                                                           <  1:760124/146‑1 (MQ=255)                                                                             cccTGTCCCGCG‑GTGTTG‑CGCAGAAAGCTTCCGAAAAGACCGTTGCGCCCTACCGTCCCCAGCCCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCAt                                                                          <  1:187611/147‑1 (MQ=255)                                                                              ccTTTTCCGCC‑GTGTTGCCGAGAAATGCGGCTA‑AAATCCCGGTGCGTCCTACCGTCCCCAGACCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAAt                                                                      <  1:710949/151‑1 (MQ=255)                                                                              ccTGTCCCGCC‑GTGTTGCCGGGGAATGCTGCTA‑AAATACCGTTGCGTCCTACCGTGCCCAGACCCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAAt                                                                      <  2:413289/151‑1 (MQ=255)                                                                                tGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGc                                                                    >  1:253583/1‑151 (MQ=255)                                                                                 tgtcccccgtggtcgCCAAGAAATGCTGCTA‑AAATACCGGGGCGCCCTACCGCCCCAAGACCCAAAC‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCAtt                                                                         <  2:496886/132‑1 (MQ=255)                                                                                  tcccgcc‑gTTTTGCCGAGAAATGCTGCTG‑AAATCCGGGGGCGCCCTCCCGTCCCCAGACCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGc                                                                    <  2:831499/147‑1 (MQ=255)                                                                                       cc‑gTTTTCCCGGGAAATGCTATTA‑AAATGCCGGTGCTTCCCACCTCGCCCAGAACCAAAT‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAAtt                                                              <  1:719902/150‑1 (MQ=255)                                                                                           tgtTGCCGAGGAATGCTGCTA‑AAATACCGGTGCGTCCTCCCGTCCCCAGAACCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCt                                                           <  1:186842/150‑1 (MQ=255)                                                                                             ttcccGAGAAATGCTGCTA‑AAATCCCGGGGCGCCCTCCCGTCCCCAAACCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCt                                                           <  2:55865/145‑1 (MQ=255)                                                                                                ccGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGc                                                      >  1:14494/1‑150 (MQ=255)                                                                                                  gagaAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGACCCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAg                                                    <  1:387837/150‑1 (MQ=255)                                                                                                    ggAATGCTGTAG‑AAATACCGTTGGTTCCTACGGTCCCCAAAACCACAC‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAggg                                                  <  2:432863/150‑1 (MQ=255)                                                                                                    gaaaTGCTGCTG‑AAATCCCGGTGCGTCCTCCCGTGCCCAGACCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAgg                                                   <  2:773315/147‑1 (MQ=255)                                                                                                     gAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCt                                                >  1:805663/1‑151 (MQ=255)                                                                                                       aTGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCttt                                              >  2:903137/1‑151 (MQ=255)                                                                                                          ctgcta‑AAATCCCGGGGCGTCCTCCCGTGCCCAGACCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCa                                            <  2:253557/150‑1 (MQ=255)                                                                                                                  aaTACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGc                                    >  2:625422/1‑151 (MQ=255)                                                                                                                      ccGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAAt                                >  1:226026/1‑151 (MQ=255)                                                                                                                            cGTCCTCCCGTCCCCAAACCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGtt                             <  2:805598/148‑1 (MQ=255)                                                                                                                            cGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGcc                          >  1:636582/1‑151 (MQ=255)                                                                                                                                   cggtcccaaaacccACAT‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACg                      <  1:343511/136‑1 (MQ=255)                                                                                                                                        cccAGACCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCAt              <  1:376235/151‑1 (MQ=255)                                                                                                                                         ccAGACCCACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATc             <  1:43427/151‑1 (MQ=255)                                                                                                                                          cAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCg            >  1:366127/1‑151 (MQ=255)                                                                                                                                            aaaCCCCAATTCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCg            <  2:566840/149‑1 (MQ=255)                                                                                                                                               ccACAA‑CCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTtata       <  1:563706/151‑1 (MQ=255)                                                                                                                                                    a‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACGc  >  1:494399/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCC‑GTGTTGCCGAGGAATGCTGCTG‑AAATACCGGTGCGTCCTACCGTGCCCAGAACCACAA‑TCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACGC  >  NC_000913/1397684‑1397979

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