Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,434,882	0	T	G	48.4%	93.5 / 14.6	63	Y252S (TAC→TCC)	truA	tRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (4/28);  new base G (30/0);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.04e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  NC_000913/2434739‑2435030                                                                                                                                                 |                                                                                                                                                       cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGGCAGGGGAATCCAc                                                                                                                                                 >  1:307476/1‑151 (MQ=255) cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATA‑CCGGTCAGGGGAATTCAc                                                                                                                                                 >  1:485585/1‑151 (MQ=255) cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATA‑CCGGTCAGGGGAAAcc                                                                                                                                                   >  1:370417/1‑149 (MQ=255)  atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACAAA‑CCGGGCAGGGGAATTCAcc                                                                                                                                                >  1:824424/1‑151 (MQ=255)      ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGGCAGGGGAATCCAcc                                                                                                                                                >  1:328037/1‑147 (MQ=255)              gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATa                                                                                                                                     <  2:944530/150‑1 (MQ=255)               ccTTGTCCCCTTCACTTTCCAGGTTGTACCTACGTTGGCTCAAGCCTGAAAGTTCCCGGGCATTATTCCCTTTAGCCCCCCAGAAATAGCGGGCCCATT‑GCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATa                                                                                                                                     <  1:569059/148‑1 (MQ=255)                cTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTTGAAGAACAAA‑ACGGTCAGGGGAAAACACCCCGACCAGATA‑Ca                                                                                                                                  >  2:232305/1‑150 (MQ=255)                   gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑Gcc                                                                                                                              <  2:327516/150‑1 (MQ=255)                         gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑ACCGGCAGGGGAATTCACCGCGACCCGAAA‑CA‑ACCCtt                                                                                                                           >  1:833439/1‑147 (MQ=255)                           cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCt                                                                                                                      <  2:361158/150‑1 (MQ=255)                                 tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGAAGAACATA‑ACCGTCAGGGGAAGCCAACCCGACCCGAAA‑CA‑AGCCCTCCGCtttt                                                                                                                   >  2:448345/1‑147 (MQ=255)                                  cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCg                                                                                                               <  2:282251/150‑1 (MQ=255)                                    ggCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTtg                                                                                                            <  1:706273/151‑1 (MQ=255)                                     gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATGGGCGGTTTTTGAAAATCATA‑CCGGCCCGGGGAAACCCCCCCGGCCAACAA‑CA‑GCCCCTCCCCtttt                                                                                                                   >  2:453092/1‑143 (MQ=255)                                      cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATA‑CCCGGCCGGGGAAACCCCCCCGACCAGATACCA‑CCCCTTCC‑CTTTTTcc                                                                                                                >  2:959262/1‑145 (MQ=255)                                        gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTtgctg                                                                                                         <  2:846466/150‑1 (MQ=255)                                          aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑ACCGGCCGGGGAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTTCCGTTTCTgcc                                                                                                       >  1:209658/1‑150 (MQ=255)                                              aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCa                                                                                                   <  2:586648/150‑1 (MQ=255)                                              aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAACGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTgccgcc                                                                                                    <  1:933933/149‑1 (MQ=255)                                               cGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGc                                                                                                 <  1:627316/151‑1 (MQ=255)                                                gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCg                                                                                                <  2:993834/151‑1 (MQ=255)                                                 ttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATCATA‑CCCGTCAGGGGAATCCACCGCGACCAGATA‑AA‑ACCCCTCCGCTTTTTCCGGttctggcccc                                                                                                    >  2:515951/1‑141 (MQ=255)                                                      tACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGGCAGGGGAAACCACCGCGACCAGATA‑CA‑ACCCTTCCCCTTTTTCCGTTGCTACCCCCCGCGTTCTGt                                                                                          >  1:525189/1‑151 (MQ=255)                                                         aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTGTTGCCGTTGCTGCCGCCAGCGTTCTGTCCt                                                                                       <  2:877134/151‑1 (MQ=255)                                                           ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                                     <  1:864657/151‑1 (MQ=255)                                                           ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                                     <  1:1032584/151‑1 (MQ=255)                                                              gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCAGAAA‑AA‑GCCCTTCCGCCTTTTCCGTTTGTGCCCCCAGCGTTCTGTCCttt                                                                                     >  2:330512/1‑148 (MQ=255)                                                                       cGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGCAGATAATA‑CCCGTCACGGGAATTCACCGCGACCCCATA‑TA‑AGCCTTCCCCCTTTGCCCGTGCTGCCCCCAACGTTCTGTCCTTTTCCGCCa                                                                              >  2:531815/1‑146 (MQ=255)                                                                         aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCCTA‑ACGGGCAGGGGAATCCACCGCGACCCGATA‑CA‑GCCCTTCCCCCTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCCTTTTCGCCAGCAActc                                                                       >  1:206801/1‑151 (MQ=255)                                                                          gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGAGAAAATACCCGGCCAGGGAAAT‑CACCGCGACCCAATA‑CAGGCCCTTCC‑CTTTTTCCGGTGGCGCCCCCAGCGGTCTGTTCCTTtccgccccccacac                                                                       >  2:8018/1‑142 (MQ=255)                                                                           cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATA‑CCCGGCCGGGGAATCCACCGCGACCAGAAA‑CA‑ACCCCTCCGCTTTTGCCGCTGCTGCCGCCCGCGGTCTGTTCTTTTccgcccgcca                                                                          >  2:87784/1‑144 (MQ=255)                                                                           cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGTCAGGGGAATTCACCCCGACCAGATA‑AA‑ACCCTTCCGCCTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCAGCAAct                                                                        >  2:807790/1‑148 (MQ=255)                                                                            attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCCGGCAGGGGAATACACCCCCACCAGATA‑CA‑ACCCCTCCCCCTTTGCCGTTGGTGCCCCCCGCGGTTTGTTCCTTGCCCCCCGccaacc                                                                       >  2:134203/1‑143 (MQ=255)                                                                              tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAACATA‑CCGGTCAAGGGAATACACCCCGACCCGGAA‑CA‑GCCCTTCCGCTTTTTGCGTTTCTGCCGCCCGCGTTCTGTCCTTTTGCGCCAGc                                                                            >  1:114201/1‑141 (MQ=255)                                                                                ttCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCCGTCCGGGGAATCCACCCCGACCAGATA‑AA‑GCCCCTCCGCTTTTGCCCGTGCTGCCGCCAGCGGTCTGTCCCTTGCCGCCCGCAACCCTGCt                                                                   >  2:1027593/1‑148 (MQ=255)                                                                                    cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTTGGGGATCATA‑CCGGGCAGGGGAATCCACCCCGCCCAAATA‑CA‑GCCCCTCCCGTTTTTCCCCTGCTTCCCCCCACGTTCTTTTCTTTTGCGCCCCCGACACCGGTATCCAGc                                                            >  2:1024616/1‑151 (MQ=255)                                                                                       tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATA‑CCGGGCAGGGGAATCCACCGCGGCCAGATA‑CA‑GCCCCTCCGCTTTTGCCGTTGCTGCCCCCAGCGGTCTGTCCTTTTCCGCCAGCCAATCTGCTAt                                                                 >  1:466015/1‑143 (MQ=255)                                                                                          tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCAGATA‑CA‑GCCCCTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTGCCGCCCGCAACTCTGCTATCCAGCtctc                                                        >  1:1061458/1‑149 (MQ=255)                                                                                           ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGGCCGGGGAATCCACCGCGACCAGATA‑CA‑ACCCTTCCGCCTTTTCCGTTTCTGCCCCCAACGTTCTGTCCTTTGCCCCCAGCAACTCTGCTATCCCGCTCTCCgg                                                     >  1:980133/1‑151 (MQ=255)                                                                                           ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGAACATA‑CCGGTCCGGGGAATCCACCCCGACCAGATA‑CA‑GACCCTCCGCTTTTTCCGGTGCTGCCCCCCGCGTTCTGTCCTTTTCCGCCCGCCACTCTGCTATCCAGCTCTcc                                                       >  2:798854/1‑149 (MQ=255)                                                                                              ccAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGGCAGGGGAATCCACCGCGACCAGATA‑CA‑GCCCTTTCCCTTTTTCCGTTTCTGCCCCCAGCGTTCTTTCCTTTTCCGCCCGCAACCCTGCTTTCCAGCTCTcccg                                                     >  1:446710/1‑146 (MQ=255)                                                                                              ccAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTg                                                  <  2:432428/151‑1 (MQ=255)                                                                                               cAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCAGATA‑CA‑GCCCTTCCCCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCTCCCGCTgg                                                 >  1:677362/1‑151 (MQ=255)                                                                                                aGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGGCAGGGGAATCCACCGCGACCAGATA‑CA‑GCCCCTCCCCTTTTTCCGTTGCTGCCCCCCGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATTCCGCTCTCCCGCCgg                                                 >  1:905333/1‑150 (MQ=255)                                                                                                                      gTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg                          <  2:614739/151‑1 (MQ=255)                                                                                                                          tGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                       <  1:336415/150‑1 (MQ=255)                                                                                                                          tGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                       <  2:592529/150‑1 (MQ=255)                                                                                                                            gAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGcc                    >  2:535029/1‑151 (MQ=255)                                                                                                                            gAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGcc                    >  1:974463/1‑151 (MQ=255)                                                                                                                                aTCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                <  1:294997/151‑1 (MQ=255)                                                                                                                                  cATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat              <  2:865400/151‑1 (MQ=255)                                                                                                                                   aTA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAata             <  2:417012/151‑1 (MQ=255)                                                                                                                                    tA‑CCGGTCAGGGGAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAatat            >  2:448052/1‑151 (MQ=255)                                                                                                                                       ccGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc          <  1:963941/151‑1 (MQ=255)                                                                                                                                       ccGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc          <  1:524092/151‑1 (MQ=255)                                                                                                                                        cGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTcc         <  2:626849/151‑1 (MQ=255)                                                                                                                                         ggTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACCATATTCCt        >  2:782053/1‑151 (MQ=255)                                                                                                                                          gTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg       <  2:314016/151‑1 (MQ=255)                                                                                                                                             aGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    <  2:413955/151‑1 (MQ=255)                                                                                                                                             aGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    >  1:358881/1‑151 (MQ=255)                                                                                                                                             aGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc    <  1:788654/151‑1 (MQ=255)                                                                                                                                               ggTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat  <  1:709875/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                       CATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATA‑CA‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAT  >  NC_000913/2434739‑2435030

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