BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,481	1	.	A	29.2%	193.8 / 140.0	89	intergenic (+18/‑16)	metZ/metW	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand): ref base . (6/57); new base A (0/26); total (6/83)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-03

Only 100 of 109 total aligned reads displayed.

GGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGC  >  NC_000913/2947331‑2947596

ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑Gt                                                                                                                      <  1:1421117/151‑1 (MQ=255)

 gctaataattaaGGTTGCATGAATTCTCTTTCGTTGTAGAGTTGGCATCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aa                                                                                                                   <  1:1591409/141‑1 (MQ=11)

  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aa                                                                                                                   <  2:1560291/151‑1 (MQ=255)

  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aa                                                                                                                   <  7:1778705/151‑1 (MQ=255)

  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aa                                                                                                                   <  6:886830/151‑1 (MQ=255)

  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aa                                                                                                                   <  5:1533/151‑1 (MQ=255)

   gTAAAAAACAGGGTGCATGAAAACCCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aaa                                                                                                                  <  5:776105/151‑1 (MQ=255)

   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aaa                                                                                                                  <  2:1282032/151‑1 (MQ=255)

   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑aaa                                                                                                                  <  6:1040508/151‑1 (MQ=255)

   gTAAAAAAAAGGGTCCATCAAAAAGACATCGTAGTAGAGAGCGCATCCGCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaa                                                                                                                 <  3:1902577/151‑2 (MQ=11)

    taacacacacggtGCATCAAAACGCTAGGGGTGTATAGTGCGCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaa                                                                                                                 <  3:1854923/138‑2 (MQ=11)

    aaaaaaaaaGCGTGCATGAAATAGATAGCGGAGTATATTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaag                                                                                                                <  1:1784112/150‑3 (MQ=14)

     ataactaaGGATGCAGGAAAAAGCGTTCGGTGTAGAGTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaag                                                                                                                <  8:1755479/144‑3 (MQ=11)

     aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGc                                                                                                                <  1:6547/151‑1 (MQ=255)

     aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGc                                                                                                                <  6:792254/151‑1 (MQ=255)

     aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TGAAATTTGAAGAA‑GT‑aaaagc                                                                                                               >  3:1514192/1‑148 (MQ=255)

     aaaaaaaaCGGTGATTGAAAACGCGAGCGCTGAATAGTCCCCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagc                                                                                                               <  6:1743499/151‑4 (MQ=11)

      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCa                                                                                                               <  4:1261898/151‑1 (MQ=255)

      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCa                                                                                                               <  7:35215/151‑1 (MQ=255)

      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCa                                                                                                               <  1:1562696/151‑1 (MQ=255)

      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCa                                                                                                               <  6:193149/151‑1 (MQ=255)

      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCa                                                                                                               <  3:30693/151‑1 (MQ=255)

      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCa                                                                                                               <  3:38357/151‑1 (MQ=255)

      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTTATGAA‑GT‑aaaagc                                                                                                               >  1:1417958/1‑147 (MQ=255)

      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagca                                                                                                              <  6:1425278/151‑5 (MQ=255)

      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGAATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCa                                                                                                               <  1:1247342/151‑1 (MQ=255)

       aaacaaGGCTGCATGAAATCCCGAGCGGAGTATAGTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagcag                                                                                                             <  7:691415/151‑6 (MQ=14)

       aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAg                                                                                                              <  6:473738/151‑1 (MQ=255)

       aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAg                                                                                                              <  6:973238/151‑1 (MQ=255)

       aaaaaaGGTTGAATGAAAAACATATCGGTGTATATTAGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GT‑aaaagcag                                                                                                             <  2:382529/151‑6 (MQ=11)

        accaagggtcaaTGAATAAGCGAGCGTAGTATACTGCCAATCCGCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGt                                                                                                             <  6:570188/140‑1 (MQ=1)

        aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGt                                                                                                             <  4:1850971/151‑1 (MQ=255)

        aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGt                                                                                                             <  8:1463258/151‑1 (MQ=255)

        aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGt                                                                                                             <  3:238575/151‑1 (MQ=255)

          cccgccgccATGAAATCGCTTGCGCAGTATAGAGCCCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTAc                                                                                                           <  1:1010816/143‑1 (MQ=11)

          aacggTTGCATGAAATCTCTTGCGGAGTATTTTTCGATTCCAGTGACGCGGGGTGTAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTAc                                                                                                           <  5:1343500/148‑1 (MQ=11)

          aacggTTGCATGAAAAAGAGAGAGGAGTATAGTCCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTAc                                                                                                           <  6:1318923/148‑1 (MQ=14)

          aaaGGTTGCATGAAAACTCGAGCGGAGTATAGACCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTAc                                                                                                           <  5:1448738/151‑1 (MQ=14)

           aaGGTTGCATGAATACTCGAGCGTAGTATTGTCGGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACg                                                                                                          <  5:1184268/151‑1 (MQ=11)

           aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACg                                                                                                          >  3:113423/1‑151 (MQ=255)

           aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTGAATAAA‑TA‑AAAGCAGTACg                                                                                                          >  1:1884958/1‑151 (MQ=18)

           aaGGTTGAATGAAAACGAGATCGGAGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACg                                                                                                          <  6:1418475/151‑1 (MQ=14)

            aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACgg                                                                                                         <  3:1223896/151‑1 (MQ=255)

             gggtGCATGAAAACGCTAGGGGAGTATATTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGt                                                                                                        <  1:1071335/148‑1 (MQ=14)

             ggTTGCATGAAATCGCGAGCGGAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGt                                                                                                        <  6:481621/151‑1 (MQ=14)

              gtggCATGAAAACGCCTGCGCTGTATTCTCCGCTGCCCCGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTg                                                                                                       <  3:843318/148‑1 (MQ=11)

              gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTg                                                                                                       <  6:760684/151‑1 (MQ=255)

              gTTGAATGAAATCTCGTGCGTAGTATAGTGCGCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTg                                                                                                       <  3:802056/151‑1 (MQ=11)

               gtgaatcaaaaCTAGATCGCAGTATAGTCGGAATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGa                                                                                                      <  6:495590/144‑1 (MQ=11)

                tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGa                                                                                                      <  3:590918/150‑1 (MQ=255)

                tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGa                                                                                                      <  3:1605865/150‑1 (MQ=255)

                tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGa                                                                                                      <  8:751745/150‑1 (MQ=255)

                 gCATGAAACCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGAcg                                                                                                    <  6:1472716/151‑1 (MQ=255)

                 gCATGAAAACGCGATCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGATCAAATCCGGCGCCCGCAACCAATTAAAATTTGATGAA‑GT‑ACAGCAGTACGGTGAcg                                                                                                    <  3:1748495/151‑1 (MQ=255)

                 gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGAcg                                                                                                    <  5:1163363/151‑1 (MQ=255)

                 ccaggAAAACGCGATGGGAGTATTGTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGAc                                                                                                     <  6:384058/146‑1 (MQ=11)

                  aatcaaaaCGAGAGGGTTGTATTGTGCGCCTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGAcgc                                                                                                   <  6:238/147‑1 (MQ=11)

                  aaggAAAACGCGAGGGGTGTATTGTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGAcgc                                                                                                   <  8:1515199/148‑1 (MQ=11)

                  aaTGAAATCTCTAGCGTAGTATATTGCGCATCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGAcgc                                                                                                   <  6:1885740/150‑1 (MQ=1)

                   aTGAAAACGCGAGCGGAGTATAGTGCGGATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGAcgcg                                                                                                  <  3:1347565/151‑1 (MQ=255)

                   aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGAcgcg                                                                                                  <  3:1338398/151‑1 (MQ=255)

                   aTGAAAACGAGCTGGGTGTATAGTCCGCTTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGAcgcg                                                                                                  <  3:235205/151‑1 (MQ=11)

                   aTGAAAACGAGAGCGGAGTAGAGTCGGAATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGAcgcg                                                                                                  <  8:1902060/151‑1 (MQ=11)

                    ggAAAACGCTAGCGGAGAATAGTGCGCTTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGACGCgg                                                                                                 <  6:896133/150‑1 (MQ=11)

                     gaaaactctattggAGTATATTTGTCATACTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATGCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACTGTGACGCgg                                                                                                 <  3:1116599/138‑1 (MQ=11)

                     gAAAACTCGATCGTAGTATATTTCGAATCCAGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTa                                                                                                            <  1:1537737/139‑1 (MQ=11)

                     gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAGGAG‑GT‑AAAGCAGTACGGTGAcgcgcg                                                                                                >  2:479657/1‑149 (MQ=255)

                     gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGGA‑AAAGCAGTaagatcggtcgcg                                                                                                >  6:1538256/1‑139 (MQ=14)

                      aattagagttggtaGTATTGTGGGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGACGCgggg                                                                                               <  8:768622/138‑1 (MQ=11)

                      aaatagcgtttgttgTATAGTTGGCATCCCGGGACGCGGGGTGGAGGAGCGTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑GTAAAAGCAGTACGGTGACGCgggg                                                                                               <  3:800989/137‑1 (MQ=11)

                       acaCGCGCGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGt                                                                                              <  1:785212/149‑1 (MQ=255)

                       aatcgagatcGGAGTATAGTGCGCATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAA‑Gt                                                                                                                      <  6:1487916/118‑1 (MQ=14)

                       aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGt                                                                                              <  8:603271/151‑1 (MQ=255)

                        aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTg                                                                                             <  3:920468/151‑1 (MQ=255)

                        aaCGCGAGCGGAGTATAGTGCGAATCCACGGACGCGTGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAATGTCGTCGGTGCAAATCCTGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTg                                                                                             <  5:40018/151‑1 (MQ=255)

                           gcgAGCGGAGTATAGTGCGCATCCACGGACGCTGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGag                                                                                          <  7:753468/151‑1 (MQ=255)

                           gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGag                                                                                          <  2:838368/151‑1 (MQ=255)

                           gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGag                                                                                          <  6:1345542/151‑1 (MQ=255)

                             gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGCCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGagca                                                                                        <  4:1697201/151‑1 (MQ=255)

                              cgCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGGCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGagcag                                                                                       <  8:422552/150‑1 (MQ=255)

                                 gtaggatagggatcATCCCCGGAAGCTGGGTGGAGCAGCCTGGTAGGTCGGCGGTCTCATAACCCGAAGGGCGGCGGTTCAAATCAGTCCCCCGCAACCAATTATATTTTGATGAA‑GT‑ATCGGAGTACGGGGACGCGGGGTGGAGCAGCCt                                                                                    <  3:902323/138‑1 (MQ=9)

                                   aGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCTGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTgg                                                                                  <  2:12718/151‑1 (MQ=255)

                                      atacggataaTCCACGGACTCGGGTGGGATCAGCCTTTTAGCTCTTCGGGATCAGAACCCTAAGGGCGTAGGTTCAAATCCGGCCCCCGCAACCAATAAAACTTTGATGAA‑GA‑AACCCAGTACCGTGACGCGGGGTGGAGCAGCCTGGTAg                                                                               <  4:279973/142‑1 (MQ=9)

                                      ataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAGATTTCGGAA‑GA‑AAAGCAGTAGGCTGTCGCGtggtggccccgtctgtgcg                                                                               >  6:624958/1‑138 (MQ=17)

                                      ataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAg                                                                               <  1:955376/151‑1 (MQ=255)

                                        agggagcATCCCCGCACGCGGTGGGGATGAGCCTGGTCGCTCGGCGGGCACACAAGCCGAAGGTCGGCGGTTCAAATCCCCCCCCCGCAACCAATTACACATTGCTGAA‑TT‑CAAGCACTACCCTGACGCGGGGTGGAGCAGCCTGGTAGCt                                                                             <  2:1209899/146‑1 (MQ=9)

                                           gcgcATGGACCGACTCGTGGTGGAGTAGCCTTGTAGCTCTTCGGGATCATATCCCGAATCGCGGAGCTTCAAATACCGCCGGCGCAAGCAAATACAAATAGTAGAA‑GA‑AACGGACTACCGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgt                                                                          <  4:839812/151‑1 (MQ=2)

                                           gcgcATCCACGGAAGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTGATATCCGTCCGCCGCAACGAATTAAAATTTGATGAA‑GT‑TAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgt                                                                          <  1:1598324/151‑1 (MQ=17)

                                            cgcATCCACGGAAGCGGGTGGGAGCAGCCTGGTAGCTCGTCGGGCTCAAAACCCGAAGGTCGTCGGTGCAAATCAGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCACTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                         <  8:1539715/151‑1 (MQ=17)

                                               aTCCTCTGACGAGTGTTGGAGCAGCCTGGTCGCTCGTCGGTCTCATAACCCGAAGGTCGTCGGTTCAAATCATTCCCGCGCAACGAATTAAAATTTGATGAA‑GT‑AACGCAGTACTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCggg                                                                      <  5:1844640/151‑1 (MQ=14)

                                               aTCCTCGGACGATGGGTGGAGAAGGGTGGTATCTCTTCGGTCTCTGAACCCTAAGGTCTTCGGTTCAAATCCTCCCGCCGCAACCAATTAACCTTTGCTGAA‑TT‑CTCGCACTACCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCggg                                                                      <  2:494529/151‑1 (MQ=9)

                                                  cACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTc                                                                   <  8:115112/151‑1 (MQ=255)

                                                              ggggtgcagcCTGGTCGCTCTGCTGTCTCACAAGCCGAAGGTCGGCGCGTCCAAACCCCCCGGCGGAACCAATAATCCATTGCTGAA‑GT‑TACGGATTACCGGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAgg                                                       <  5:230143/146‑1 (MQ=2)

                                                                             cgatctgcgctataaTAAGCCTAATGGCGTCGGGTCCAATCATCCCCCCGCAACGAATAAAACATTGCTGAA‑TA‑CAAGCAGTACCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAAt                                         <  2:643733/136‑1 (MQ=2)

                                                                               gtgtgcgGGCTCATAACCCGAATCTCGGCTGGTCATATCCTGCCGCCTCAAGCAATTAAAATGTGATGTA‑TA‑ATCGCAGTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCg                                      <  4:677758/146‑1 (MQ=2)

                                                                                  gtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtcg                                                                        >  7:447576/1‑114 (MQ=255)

                                                                                           aTAACCCGAAGGTCGTCGGTTCAAATACGTCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAcc                          <  6:634445/151‑1 (MQ=255)

                                                                                               ccggAAGGTCGTAGGTGGATATCCGACCGCTGCAACGAATTATATTTTGATGTA‑TT‑CTAGCAGTACGGTGACGCGGGGTGGAGCAGCGTGGTAGCTGGTCGGGCTCATAAGCCGAAGGTCGTCGGTTCAAATCCGGCCCGCGCAACGAAt                       <  7:320598/147‑1 (MQ=2)

                                                                                                     gggcggcgGTTCAACTCCGGCCCCCGCAAACAATTAAACATGGATGAA‑GT‑AACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                <  8:113695/145‑1 (MQ=14)

                                                                                                                   atccttcccccGCAATCAATTAAAATTTGATGAA‑TT‑AACGCATTACCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGc  <  3:1406510/145‑1 (MQ=18)

GGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAA‑GT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGC  >  NC_000913/2947331‑2947596

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