BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,207,023	A→G	D346G (GAC→GGC)	yehP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,207,023	0	A	G	99.1%	706.9 / ‑4.7	229	D346G (GAC→GGC)	yehP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (129/97); total (131/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 242 total aligned reads displayed.

TACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATC  >  NC_000913/2206933‑2207112

    gACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGa                                                                                     <  1:3563979/93‑1 (MQ=255)

    gACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGa                                                                                     >  2:947427/1‑93 (MQ=255)

       gCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATAc                                                                                  <  1:7437242/93‑1 (MQ=255)

       gCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATAc                                                                                  >  2:3556526/1‑93 (MQ=255)

       gCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATAc                                                                                  <  2:6349297/93‑1 (MQ=255)

          tCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGc                                                                               >  2:5195819/1‑93 (MQ=255)

          tCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGc                                                                               >  2:3474290/1‑93 (MQ=255)

          tCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGc                                                                               >  2:3751538/1‑93 (MQ=255)

          tCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGc                                                                               >  1:5944546/1‑93 (MQ=255)

           cAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGcc                                                                              >  2:2873464/1‑93 (MQ=255)

            aGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGccc                                                                             >  1:4940882/1‑93 (MQ=255)

            aGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGccc                                                                             <  2:7431850/93‑1 (MQ=255)

            aGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCACAACACCTTGCTATGGCCGCGATACGGccc                                                                             >  2:6973586/1‑93 (MQ=255)

             ggTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCa                                                                            <  2:2518954/93‑1 (MQ=255)

              gTGAAAAAGTGTGTCCAGAGCGTCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAg                                                                           <  2:3864610/93‑1 (MQ=255)

              gTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAg                                                                           <  1:4032592/93‑1 (MQ=255)

               tGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAgg                                                                          >  1:6059762/1‑93 (MQ=255)

               tGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAgc                                                                          >  1:2071994/1‑92 (MQ=255)

                gAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGgc                                                                         >  2:4458609/1‑93 (MQ=255)

                 aaaaaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGgcg                                                                        >  1:1070794/1‑93 (MQ=255)

                  aaaaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGgcgc                                                                       >  1:1694840/1‑93 (MQ=255)

                    aaGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTg                                                                     <  2:6875380/93‑1 (MQ=255)

                     aGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGTCCGCGATACGGCCCAGGCGCTgg                                                                    >  2:4333822/1‑93 (MQ=255)

                     aGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTgg                                                                    <  1:3197799/93‑1 (MQ=255)

                       tgtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGtt                                                                  >  1:5883849/1‑93 (MQ=255)

                        gtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTa                                                                 >  2:1669649/1‑93 (MQ=255)

                        gtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTa                                                                 >  1:550020/1‑93 (MQ=255)

                         tgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTaa                                                                <  2:4606701/93‑1 (MQ=255)

                         tgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTaa                                                                <  2:2294662/93‑1 (MQ=255)

                         tgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTaa                                                                <  2:2645011/93‑1 (MQ=255)

                           tCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATg                                                              <  2:4864135/93‑1 (MQ=255)

                            ccAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGt                                                             <  1:7865109/93‑1 (MQ=255)

                            ccAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGt                                                             >  2:756550/1‑93 (MQ=255)

                            ccAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGt                                                             >  2:710017/1‑93 (MQ=255)

                             cAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTc                                                            >  2:7375328/1‑93 (MQ=255)

                             cAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTc                                                            <  2:16963/93‑1 (MQ=255)

                             cAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTc                                                            >  1:1360006/1‑93 (MQ=255)

                              agagCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCg                                                           >  1:2019720/1‑93 (MQ=255)

                              agagCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCg                                                           >  1:4686108/1‑93 (MQ=255)

                               gagGGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTCATGGCgg                                                          >  2:57137/1‑93 (MQ=255)

                                 gCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGgcg                                                        >  1:1778695/1‑93 (MQ=255)

                                    gCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGcac                                                     >  2:1130202/1‑93 (MQ=255)

                                     cATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGcaca                                                    >  2:2810564/1‑93 (MQ=255)

                                     cATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGcaca                                                    <  1:825392/93‑1 (MQ=255)

                                       tCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATTTCGGCGCACaaa                                                  >  2:6478581/1‑93 (MQ=255)

                                       tCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaaa                                                  >  1:995592/1‑93 (MQ=255)

                                       tCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaaa                                                  <  1:1104999/93‑1 (MQ=255)

                                       tCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaaa                                                  <  1:1392109/93‑1 (MQ=255)

                                       tCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaaa                                                  >  1:3317308/1‑93 (MQ=255)

                                       tCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaaa                                                  >  1:7160476/1‑93 (MQ=255)

                                       tCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACaaa                                                  >  1:4582985/1‑93 (MQ=255)

                                        cAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGTCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAAt                                                 >  2:6440895/1‑93 (MQ=255)

                                             tGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAgccg                                            >  2:5523568/1‑93 (MQ=255)

                                                tGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAGATAGCCGCCa                                         >  2:71920/1‑93 (MQ=255)

                                                 gggACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCAt                                        <  1:523372/93‑1 (MQ=255)

                                                   gACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGa                                      >  1:5637323/1‑93 (MQ=255)

                                                   gACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGa                                      <  1:3768386/93‑1 (MQ=255)

                                                    aCTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGAc                                     >  1:5422845/1‑93 (MQ=255)

                                                     cTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACg                                    >  1:6627547/1‑93 (MQ=255)

                                                      tGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGc                                   >  2:2062912/1‑93 (MQ=255)

                                                      tGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGc                                   <  2:5341234/93‑1 (MQ=255)

                                                           gcgcTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGc                              >  1:1846623/1‑93 (MQ=255)

                                                            cgcTCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCg                             <  1:3369059/93‑1 (MQ=255)

                                                               tCGATAGCGCCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCACACAAATAGCCGCCATGACGCCGGGCGAGc                          >  2:1577866/1‑93 (MQ=255)

                                                               tCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                          <  2:1964573/93‑1 (MQ=255)

                                                               tCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                          <  1:3177321/93‑1 (MQ=255)

                                                               tCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                          >  2:7704534/1‑93 (MQ=255)

                                                               tCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                          <  1:6233538/93‑1 (MQ=255)

                                                               tCGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                          >  1:1482407/1‑93 (MQ=255)

                                                                cGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCt                         >  2:2694168/1‑93 (MQ=255)

                                                                cGATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCt                         <  1:5135007/93‑1 (MQ=255)

                                                                 gATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTg                        >  2:6300826/1‑93 (MQ=255)

                                                                 gATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTg                        <  2:4243125/93‑1 (MQ=255)

                                                                 gATAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTg                        >  2:1411667/1‑93 (MQ=255)

                                                                  aTAGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTgg                       >  2:7502870/1‑93 (MQ=255)

                                                                    aGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGca                     >  1:17085/1‑93 (MQ=255)

                                                                    aGCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGca                     >  2:5399530/1‑93 (MQ=255)

                                                                     gCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcat                    >  1:674323/1‑93 (MQ=255)

                                                                     gCACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcat                    <  1:216816/93‑1 (MQ=255)

                                                                      cACCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatc                   <  2:2288589/93‑1 (MQ=255)

                                                                       aCCGCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatca                  >  2:3716039/1‑93 (MQ=255)

                                                                          gCAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATgc               >  1:3874985/1‑92 (MQ=255)

                                                                           cAACACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGc              <  2:312998/93‑1 (MQ=255)

                                                                           cAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGc              <  1:2273868/93‑1 (MQ=255)

                                                                            aaCACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCt             >  1:43872/1‑93 (MQ=255)

                                                                            aaCACCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCt             <  1:2944596/93‑1 (MQ=255)

                                                                               acCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCCTCATGGCTTGc          >  2:4084471/1‑93 (MQ=255)

                                                                               acCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCtttt          >  2:6604102/1‑91 (MQ=255)

                                                                               acCTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCtttt          >  2:360792/1‑91 (MQ=255)

                                                                                ccTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTTCg         >  1:1612357/1‑93 (MQ=255)

                                                                                 cTTGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCttttgg        >  2:3308580/1‑89 (MQ=255)

                                                                                  ttGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCggg       >  1:313547/1‑92 (MQ=255)

                                                                                  ttGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGga       <  1:2375537/93‑1 (MQ=255)

                                                                                   tGCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTTCgggg      >  2:1354714/1‑91 (MQ=255)

                                                                                    gCTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGgaga     <  2:792480/93‑1 (MQ=255)

                                                                                     cTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTTCggggga    >  1:3946861/1‑89 (MQ=255)

                                                                                     cTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTTCGGGGaa    >  1:6912901/1‑93 (MQ=255)

                                                                                     cTATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGaa    <  2:2089623/93‑1 (MQ=255)

                                                                                      tATGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGTGGGGAAt   >  1:6327256/1‑93 (MQ=255)

                                                                                       aTGGCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATc  <  1:2640073/93‑1 (MQ=255)

TACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATC  >  NC_000913/2206933‑2207112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`