BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,387,335	T→C	V68A (GTA→GCA)	b2275 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,387,335	0	T	C	98.4%	582.7 / ‑1.1	184	V68A (GTA→GCA)	b2275	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (81/100); total (82/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

Only 100 of 187 total aligned reads displayed.

ACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATC  >  NC_000913/2387245‑2387424

 cGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCAcc                                                                                        >  2:1857701/1‑93 (MQ=255)

 cGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCAcc                                                                                        <  1:1506404/93‑1 (MQ=255)

  ggAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCt                                                                                       <  1:1167159/93‑1 (MQ=255)

  ggAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCt                                                                                       >  2:3269532/1‑93 (MQ=255)

   gagaCTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTg                                                                                      >  2:4456173/1‑93 (MQ=255)

   gagaCTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAAAGGCGAATCGCACCTg                                                                                      >  2:4219426/1‑93 (MQ=255)

     gaCTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGcc                                                                                    <  2:4637905/93‑1 (MQ=255)

        cgTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCttt                                                                                 <  1:3552549/92‑1 (MQ=255)

         gTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCtttt                                                                                >  2:4009741/1‑93 (MQ=255)

         gTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCtttt                                                                                >  2:3090478/1‑93 (MQ=255)

          tttCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTc                                                                               <  2:4045122/93‑1 (MQ=255)

            tCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCa                                                                             >  2:3486557/1‑93 (MQ=255)

            tCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTACAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCa                                                                             >  2:4623454/1‑93 (MQ=255)

              ggTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGc                                                                           >  2:3448067/1‑93 (MQ=255)

               gTTTGTACTGGCCTATTGGTAACGACTATAAACTTCCTATAGTGTGCGAAACGTCCCAGGGCGAATGGCCACTCCCACCTGCCTTTTCCAGCa                                                                          <  2:4157449/93‑1 (MQ=255)

                  tGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTa                                                                       >  1:4573550/1‑93 (MQ=255)

                  tGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTa                                                                       <  2:1595522/93‑1 (MQ=255)

                    tACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTaag                                                                     <  1:1716517/93‑1 (MQ=255)

                    tACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTaag                                                                     <  1:227538/93‑1 (MQ=255)

                     aCTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTaaga                                                                    <  2:3837920/93‑1 (MQ=255)

                      cTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTaagaa                                                                   >  1:310903/1‑93 (MQ=255)

                       gggCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTaagaag                                                                  <  2:3700246/92‑1 (MQ=255)

                        ggCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGTCGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGt                                                                 >  1:3550905/1‑93 (MQ=255)

                        ggCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGt                                                                 <  2:2884816/93‑1 (MQ=255)

                        ggCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGt                                                                 <  2:3296606/93‑1 (MQ=255)

                          ccTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTc                                                               >  2:1113276/1‑93 (MQ=255)

                          ccTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTc                                                               <  1:1394464/93‑1 (MQ=255)

                           cTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCg                                                              >  2:3950752/1‑93 (MQ=255)

                            tATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCga                                                             >  1:3540207/1‑93 (MQ=255)

                             aTTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCgaa                                                            <  1:2914185/93‑1 (MQ=255)

                             aTTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCgaa                                                            <  1:4493887/93‑1 (MQ=255)

                                 gTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAAt                                                        <  2:2080114/93‑1 (MQ=255)

                                  tCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATg                                                       >  1:4373423/1‑93 (MQ=255)

                                   cAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATgg                                                      >  2:4248707/1‑93 (MQ=255)

                                    aGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGc                                                     <  1:3010916/93‑1 (MQ=255)

                                    aGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGc                                                     <  1:3494624/93‑1 (MQ=255)

                                     ggAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCt                                                    >  1:363465/1‑93 (MQ=255)

                                      gAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTg                                                   >  2:2967772/1‑93 (MQ=255)

                                      gAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTg                                                   >  1:3863651/1‑93 (MQ=255)

                                        gCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGa                                                 <  2:1529556/93‑1 (MQ=255)

                                         cATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGaa                                                <  1:308227/93‑1 (MQ=255)

                                          aTAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAg                                               <  1:2555353/93‑1 (MQ=255)

                                           tAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGt                                              >  2:2519250/1‑93 (MQ=255)

                                             aaCTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGa                                            >  1:1505333/1‑93 (MQ=255)

                                             aaCTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGa                                            >  2:1743852/1‑93 (MQ=255)

                                             aaCTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGa                                            >  1:2033556/1‑93 (MQ=255)

                                              aCTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGaa                                           >  2:3163668/1‑93 (MQ=255)

                                              aCTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGaa                                           >  2:2295988/1‑93 (MQ=255)

                                              aCTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGaa                                           <  1:2339219/93‑1 (MQ=255)

                                              aCTACCGATAGTGTGCGAAACGTACCAAGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGaa                                           <  1:653784/93‑1 (MQ=255)

                                                     ataGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAatgatgct                                    >  1:4697081/1‑91 (MQ=255)

                                                      taGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAatgatgatg                                   >  2:4294187/1‑93 (MQ=255)

                                                      taGTGTGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAatgatgatg                                   >  2:1803956/1‑93 (MQ=255)

                                                        gtgtGCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATgatgatgat                                 >  2:2687304/1‑93 (MQ=255)

                                                         tgtgCGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATc                                >  2:626738/1‑93 (MQ=255)

                                                           tacGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATccc                              <  1:3161104/91‑1 (MQ=255)

                                                             cGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCa                            <  1:2284873/93‑1 (MQ=255)

                                                             cGAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCa                            <  2:1775254/93‑1 (MQ=255)

                                                              gAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCAt                           <  2:3457126/93‑1 (MQ=255)

                                                              gAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCAt                           >  1:1858502/1‑93 (MQ=255)

                                                              gAAACGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCAt                           >  2:2823113/1‑93 (MQ=255)

                                                                aaCGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGa                         >  1:2248179/1‑93 (MQ=255)

                                                                aaCGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGa                         >  1:3750176/1‑93 (MQ=255)

                                                                aaCGTACCATAGCGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGa                         <  2:3618836/93‑1 (MQ=255)

                                                                  cGTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGaaa                       <  1:1804804/93‑1 (MQ=255)

                                                                   gTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGaaaa                      <  1:2804781/93‑1 (MQ=255)

                                                                   gTACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGaaaa                      >  2:832375/1‑93 (MQ=255)

                                                                    tACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAAt                     <  1:1132297/93‑1 (MQ=255)

                                                                    tACCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAAt                     <  2:1628745/93‑1 (MQ=255)

                                                                     aCCATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATg                    <  1:2490772/93‑1 (MQ=255)

                                                                       cATGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGc                  >  2:1379255/1‑93 (MQ=255)

                                                                        aTGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCa                 <  2:104554/93‑1 (MQ=255)

                                                                        aTGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCa                 >  2:1157268/1‑93 (MQ=255)

                                                                        aTGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCa                 <  2:1392697/93‑1 (MQ=255)

                                                                         tGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCAt                <  2:94887/93‑1 (MQ=255)

                                                                         tGACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCAt                >  2:1637260/1‑93 (MQ=255)

                                                                          ggcGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATc               <  2:2226186/91‑1 (MQ=255)

                                                                          gACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATc               <  2:463154/93‑1 (MQ=255)

                                                                          gACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCGTc               <  2:875401/93‑1 (MQ=255)

                                                                          gACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATc               <  1:981864/93‑1 (MQ=255)

                                                                          gACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATc               <  1:862237/93‑1 (MQ=255)

                                                                          gACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATc               <  1:737882/93‑1 (MQ=255)

                                                                          gACGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATc               <  2:736848/93‑1 (MQ=255)

                                                                           aCGAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCa              <  1:4144868/93‑1 (MQ=255)

                                                                             gAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCACCATGAAAATGGCATCAGc            >  1:1032902/1‑93 (MQ=255)

                                                                             gAATGGCGAATCGCCCCTGCCTTTTCCGGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGc            <  1:1357435/93‑1 (MQ=255)

                                                                             gAATGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGc            <  2:4410798/93‑1 (MQ=255)

                                                                              aaTGGCGAATCGCACCTGCCTTTTCCAGCATTAGGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCa           >  1:1427641/1‑93 (MQ=255)

                                                                              aaTGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCa           >  2:1828985/1‑93 (MQ=255)

                                                                              aaTGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCa           <  1:3314064/93‑1 (MQ=255)

                                                                               aTGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCAt          >  2:3892897/1‑93 (MQ=255)

                                                                               aTGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCAt          <  2:2817753/93‑1 (MQ=255)

                                                                               aTGGCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCAt          >  1:2592931/1‑93 (MQ=255)

                                                                                  gCGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGa       >  1:668407/1‑93 (MQ=255)

                                                                                   cGAATCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGaa      <  2:2451368/93‑1 (MQ=255)

                                                                                     aaTCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGa    <  1:4255857/93‑1 (MQ=255)

                                                                                      aTCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCTCATGAAAATGGCATCAGCATCGAAGAt   >  2:2584989/1‑93 (MQ=255)

                                                                                      aTCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGAt   >  1:4736639/1‑93 (MQ=255)

                                                                                       tCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATc  <  1:1215112/93‑1 (MQ=255)

                                                                                       tCGCACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAATGAATAATGATGATCCCCATGAAAATGGCATCAGCATCGAAGAta  >  1:1606524/1‑92 (MQ=255)

ACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACGTACCATGACGAATGGCGAATCGTACCTGCCTTTTCCAGCATTAAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATC  >  NC_000913/2387245‑2387424

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