BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	226,978	A→G	noncoding (1220/2904 nt)	rrlH →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	226,978	0	A	G	77.6%	150.5 / 112.9	165	noncoding (1220/2904 nt)	rrlH	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (19/18); new base G (53/75); total (72/93)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 166 total aligned reads displayed.

GCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTC  >  NC_000913/226878‑227077

gcgcGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTa                                                                                                     >  2:424970/1‑101 (MQ=33)

         gtGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGt                                                                                            <  2:424527/100‑1 (MQ=11)

         aTGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGt                                                                                            >  1:424978/1‑101 (MQ=11)

         aTGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGt                                                                                            >  2:424979/1‑101 (MQ=11)

          tGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTa                                                                                           <  1:424560/101‑1 (MQ=17)

          tGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTa                                                                                           <  2:424881/101‑1 (MQ=17)

            tAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATg                                                                                         >  1:424982/1‑101 (MQ=17)

            tAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATg                                                                                         >  2:424984/1‑101 (MQ=17)

            tAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATg                                                                                         >  2:424983/1‑101 (MQ=40)

             aaCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGc                                                                                        <  1:424495/101‑1 (MQ=17)

              aCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCt                                                                                       >  2:424987/1‑101 (MQ=17)

              aCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCt                                                                                       >  1:424988/1‑101 (MQ=17)

              aCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCt                                                                                       >  1:424989/1‑101 (MQ=17)

               cGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTg                                                                                      <  1:424393/101‑1 (MQ=17)

               cGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTg                                                                                      <  1:424262/101‑1 (MQ=40)

                 gggCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGa                                                                                    <  1:424573/101‑1 (MQ=17)

                 gggCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGa                                                                                    <  1:424665/101‑1 (MQ=17)

                 gggCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGa                                                                                    <  1:424758/101‑1 (MQ=17)

                  ggCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAg                                                                                   <  1:424652/101‑1 (MQ=17)

                   gCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAgg                                                                                  >  2:424992/1‑101 (MQ=17)

                    cTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGt                                                                                 <  2:424485/101‑1 (MQ=14)

                     tAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTa                                                                                <  1:424483/101‑1 (MQ=17)

                        aCCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCa                                                                             >  2:424994/1‑101 (MQ=17)

                        aCCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCa                                                                             <  2:424694/101‑1 (MQ=40)

                        aCCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCa                                                                             <  2:424439/101‑1 (MQ=40)

                         ccATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAg                                                                            >  1:424996/1‑101 (MQ=17)

                         ccATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAg                                                                            >  2:424998/1‑101 (MQ=17)

                          cATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGa                                                                           >  1:424999/1‑101 (MQ=17)

                          cATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGa                                                                           <  1:424402/101‑1 (MQ=17)

                             gCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAg                                                                         >  2:425002/1‑100 (MQ=18)

                             gCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGt                                                                        >  2:425001/1‑101 (MQ=17)

                             gCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGt                                                                        <  1:424587/101‑1 (MQ=18)

                              cACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTg                                                                       >  1:425003/1‑101 (MQ=17)

                               aCCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGc                                                                      <  1:424605/101‑1 (MQ=18)

                                 cGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGa                                                                    <  1:424469/101‑1 (MQ=18)

                                  gAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGaa                                                                   >  1:425009/1‑101 (MQ=17)

                                  gAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGaa                                                                   <  1:424476/101‑1 (MQ=18)

                                  gAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGaa                                                                   <  2:424436/101‑1 (MQ=255)

                                   aaGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAAt                                                                  >  1:425011/1‑101 (MQ=17)

                                   aaGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAAt                                                                  >  2:425012/1‑101 (MQ=17)

                                    aGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATg                                                                 >  1:425013/1‑101 (MQ=17)

                                    aGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGCAGGTATCAGAACTGCGAATg                                                                 <  1:424773/101‑1 (MQ=255)

                                     gCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGc                                                                <  1:424563/101‑1 (MQ=255)

                                      cTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCt                                                               >  1:425017/1‑101 (MQ=17)

                                      cTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCt                                                               >  2:425018/1‑101 (MQ=17)

                                      cTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCt                                                               <  1:424507/101‑1 (MQ=17)

                                       tGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTg                                                              <  1:424857/101‑1 (MQ=18)

                                       tGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTg                                                              <  1:424543/101‑1 (MQ=255)

                                        gCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGa                                                             <  1:424491/101‑1 (MQ=18)

                                          ggCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACa                                                           >  2:425021/1‑101 (MQ=17)

                                          ggCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACa                                                           <  1:424782/101‑1 (MQ=40)

                                          ggCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACa                                                           >  1:425022/1‑101 (MQ=40)

                                          ggCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACa                                                           <  1:424342/101‑1 (MQ=40)

                                            cAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATa                                                         >  1:425023/1‑101 (MQ=17)

                                            cAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATa                                                         <  1:424697/101‑1 (MQ=18)

                                             aGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATaa                                                        <  1:424390/101‑1 (MQ=18)

                                             aGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATaa                                                        <  2:424622/101‑1 (MQ=25)

                                               cGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGt                                                      <  2:424668/101‑1 (MQ=18)

                                               cGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGt                                                      <  2:424636/101‑1 (MQ=18)

                                                gACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTa                                                     >  2:425026/1‑101 (MQ=17)

                                                gACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTa                                                     >  2:425028/1‑101 (MQ=17)

                                                gACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTCGAGGTATCAGAAGTGCGAATGCTGACATAAGTa                                                     <  2:424562/101‑1 (MQ=12)

                                                 aCGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGa                                                                           <  1:424704/78‑1 (MQ=255)

                                                 aCGCTTATGCGGTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGATGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTaa                                                    >  2:425029/1‑101 (MQ=255)

                                                  cGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAAc                                                   <  2:424341/101‑1 (MQ=17)

                                                  cGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAAc                                                   <  2:424743/101‑1 (MQ=18)

                                                   gCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACg                                                  >  1:425034/1‑101 (MQ=18)

                                                   gCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACg                                                  >  1:425031/1‑101 (MQ=18)

                                                   gCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACg                                                  >  1:425033/1‑101 (MQ=40)

                                                     ttATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGAt                                                >  2:425037/1‑101 (MQ=17)

                                                     ttATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGAt                                                <  1:424553/101‑1 (MQ=18)

                                                     ttATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGAt                                                <  1:424588/101‑1 (MQ=18)

                                                     ttATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGAt                                                <  1:424608/101‑1 (MQ=18)

                                                     ttATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGAt                                                >  1:425036/1‑101 (MQ=14)

                                                     ttATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGAt                                                <  2:424570/101‑1 (MQ=18)

                                                     ttATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGAt                                                <  2:424583/101‑1 (MQ=18)

                                                     ttATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGAt                                                <  2:424754/101‑1 (MQ=255)

                                                      tATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTTAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTg                                                              <  2:424252/86‑1 (MQ=255)

                                                      tATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATa                                               <  2:424609/101‑1 (MQ=17)

                                                      tATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATa                                               <  1:424710/101‑1 (MQ=14)

                                                       aTGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATaa                                              >  2:425040/1‑101 (MQ=17)

                                                        tGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATaaa                                             <  2:424519/101‑1 (MQ=17)

                                                        tGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATaaa                                             <  1:424591/101‑1 (MQ=17)

                                                        tGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATaaa                                             <  2:424650/101‑1 (MQ=17)

                                                        tGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATaaa                                             <  2:424600/101‑1 (MQ=17)

                                                        tGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATaaa                                             <  1:424623/101‑1 (MQ=40)

                                                         gCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAg                                            <  2:424338/101‑1 (MQ=14)

                                                         gCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAg                                            >  2:425041/1‑101 (MQ=14)

                                                          cGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGc                                           <  1:424725/101‑1 (MQ=11)

                                                           gttgttGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCg                                          <  2:424770/101‑1 (MQ=38)

                                                             tgttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCggg                                        >  1:425044/1‑101 (MQ=35)

                                                              gttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCggg                                        >  2:425045/1‑100 (MQ=37)

                                                               ttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCggg                                        >  1:425046/1‑99 (MQ=37)

                                                                   gTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATg                                                                                         <  2:424381/46‑1 (MQ=255)

                                                                    tAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAg                                 <  1:424538/101‑1 (MQ=38)

                                                                                 tGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTg                                      >  2:425061/1‑83 (MQ=39)

                                                                                 tGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                    >  2:425059/1‑101 (MQ=35)

                                                                                   tAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAg                  >  1:425065/1‑101 (MQ=38)

                                                                                     aGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGAc                >  1:425067/1‑101 (MQ=38)

                                                                                        cTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa             >  2:425068/1‑101 (MQ=35)

GCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTC  >  NC_000913/226878‑227077

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`