BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,632,141	T→C	S166P (TCA→CCA)	b1545 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,632,141	0	T	C	99.6%	691.0 / ‑4.9	238	S166P (TCA→CCA)	b1545	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (103/134); total (103/135)
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.39e-01

Only 100 of 241 total aligned reads displayed.

ATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTC  >  NC_000913/1632051‑1632219

aTAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCa                                                                              <  1:1119395/93‑1 (MQ=255)

aTAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCa                                                                              <  1:4753782/93‑1 (MQ=255)

aTAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCa                                                                              <  2:48048/93‑1 (MQ=255)

 tAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAg                                                                             >  2:8309162/1‑93 (MQ=37)

 tAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAg                                                                             >  2:6315191/1‑93 (MQ=37)

 tAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAg                                                                             <  2:6261406/93‑1 (MQ=37)

 tAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAg                                                                             >  1:4803119/1‑93 (MQ=37)

  aGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAgg                                                                            <  2:7893153/93‑1 (MQ=255)

   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                           <  2:3498858/93‑1 (MQ=255)

   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                           <  1:5418111/93‑1 (MQ=255)

     aaGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtgt                                                                         <  1:764485/93‑1 (MQ=255)

      agagCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTa                                                                        >  2:7840870/1‑93 (MQ=37)

       gagCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTaa                                                                       <  2:8910118/93‑1 (MQ=255)

         gCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGt                                                                     <  2:6465420/93‑1 (MQ=255)

         gCCCGCGGCGCAGGGAAACGTTTCGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGt                                                                     <  2:5526942/93‑1 (MQ=255)

            cgcgGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGtata                                                                  >  2:7190377/1‑93 (MQ=255)

            cgcgGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGtata                                                                  <  1:5978592/93‑1 (MQ=255)

            cgcgGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGtata                                                                  <  1:4905261/93‑1 (MQ=255)

            cgcgGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGtata                                                                  <  1:927812/93‑1 (MQ=255)

             gcggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATaa                                                                 <  2:139253/93‑1 (MQ=255)

             gcggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATaa                                                                 <  1:4531923/93‑1 (MQ=255)

             gcggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATaa                                                                 <  1:6416355/93‑1 (MQ=255)

             gcggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATaa                                                                 <  2:850028/93‑1 (MQ=255)

              cggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAg                                                                >  1:2323670/1‑93 (MQ=255)

                  gcAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTCTTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGcc                                                            >  1:6498271/1‑93 (MQ=255)

                   cAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCa                                                           <  1:1624747/93‑1 (MQ=255)

                    aGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAt                                                          >  2:8590424/1‑93 (MQ=255)

                    aGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAt                                                          <  1:8030584/93‑1 (MQ=255)

                    aGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAt                                                          <  1:8713268/93‑1 (MQ=255)

                    aGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAt                                                          <  1:3556380/93‑1 (MQ=255)

                    aGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAt                                                          >  1:6919387/1‑93 (MQ=255)

                    aGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAt                                                          <  2:203529/93‑1 (MQ=255)

                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                         <  1:805078/93‑1 (MQ=255)

                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                         >  2:3918910/1‑93 (MQ=255)

                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                         >  1:697872/1‑93 (MQ=255)

                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                         >  1:6814887/1‑93 (MQ=255)

                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                         >  2:4963845/1‑93 (MQ=255)

                      ggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTg                                                        >  1:8186068/1‑93 (MQ=255)

                        aaaCGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGcc                                                      <  2:2291351/93‑1 (MQ=255)

                        aaaCGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGcc                                                      >  2:5105453/1‑93 (MQ=255)

                        aaaCGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGcc                                                      >  2:2826892/1‑93 (MQ=255)

                           cGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCgg                                                   >  2:4390375/1‑93 (MQ=255)

                             ttttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGa                                                 >  1:7476428/1‑93 (MQ=255)

                             ttttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGa                                                 >  2:3110251/1‑93 (MQ=255)

                             ttttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGa                                                 >  2:892523/1‑93 (MQ=255)

                              tttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGaa                                                <  1:4641263/93‑1 (MQ=37)

                               ttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAt                                               <  1:5634651/93‑1 (MQ=37)

                               ttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAt                                               <  2:2926936/93‑1 (MQ=37)

                                 ggTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTc                                             <  2:8372307/93‑1 (MQ=37)

                                 ggTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTc                                             >  2:1497112/1‑93 (MQ=37)

                                  gTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCa                                            >  2:943190/1‑93 (MQ=255)

                                    cGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaa                                          <  1:3171186/93‑1 (MQ=255)

                                     gACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaaa                                         <  1:3856966/93‑1 (MQ=255)

                                     gACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaaa                                         <  2:4707287/93‑1 (MQ=255)

                                      accaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                        <  1:4788085/93‑1 (MQ=255)

                                       ccaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAcc                                       >  1:605479/1‑93 (MQ=255)

                                        caccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCt                                      <  1:1898420/93‑1 (MQ=25)

                                        caccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCt                                      >  1:2921113/1‑93 (MQ=37)

                                        caccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCt                                      <  2:5228517/93‑1 (MQ=25)

                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTc                                     >  1:9157167/1‑93 (MQ=21)

                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAACTc                                     >  2:5459063/1‑93 (MQ=255)

                                          ccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcg                                    >  1:6203245/1‑93 (MQ=21)

                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgc                                   >  1:2487892/1‑93 (MQ=21)

                                            tgtgTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgcg                                  >  2:1708802/1‑93 (MQ=18)

                                             gtgtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCgg                                 <  1:739620/93‑1 (MQ=18)

                                             gtgtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCgg                                 >  1:4334221/1‑93 (MQ=21)

                                             gtgtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCgg                                 >  2:2760304/1‑93 (MQ=18)

                                               gtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCa                               <  1:1753469/93‑1 (MQ=18)

                                               gtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACATGGCggaa                               >  1:6077694/1‑91 (MQ=11)

                                                ttAAATGAAGAACAGAAACAGGCGGTCTTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGACTGCACAACCTCGCGGCaa                              >  1:4136518/1‑93 (MQ=11)

                                                ttAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCaa                              >  1:220312/1‑93 (MQ=21)

                                                 tAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCaaa                             >  2:1055994/1‑93 (MQ=21)

                                                  aaaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                            >  2:9487050/1‑93 (MQ=18)

                                                  aaaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                            >  2:8159584/1‑93 (MQ=18)

                                                   aaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAcc                           >  2:7591264/1‑93 (MQ=21)

                                                   aaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAcc                           >  2:6052460/1‑93 (MQ=18)

                                                    aTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCa                          <  1:2653395/93‑1 (MQ=18)

                                                       aagaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAttt                       <  2:6880120/93‑1 (MQ=18)

                                                          aaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTaa                    <  1:1288439/93‑1 (MQ=18)

                                                          aaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTaa                    <  1:7548451/93‑1 (MQ=18)

                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag                   >  1:5168429/1‑93 (MQ=21)

                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag                   <  2:1835588/93‑1 (MQ=18)

                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag                   <  2:4082983/93‑1 (MQ=18)

                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag                   <  2:445303/93‑1 (MQ=18)

                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag                   <  2:5962117/93‑1 (MQ=18)

                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag                   <  1:3917485/93‑1 (MQ=18)

                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga                  <  1:5699205/93‑1 (MQ=18)

                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga                  <  2:2623284/93‑1 (MQ=18)

                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga                  <  1:5488379/93‑1 (MQ=18)

                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga                  <  2:5698096/93‑1 (MQ=18)

                                                              gAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGc                >  2:7269085/1‑93 (MQ=21)

                                                               aaaCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGcc               <  2:6800296/93‑1 (MQ=18)

                                                               aaaCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGcc               >  1:6602049/1‑93 (MQ=21)

                                                                aaCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCa              >  1:4857417/1‑93 (MQ=18)

                                                                 aCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaa             <  1:9808780/93‑1 (MQ=18)

                                                                  cAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaaa            <  1:9649721/93‑1 (MQ=18)

                                                                  cAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaaa            <  2:2360164/93‑1 (MQ=18)

                                                                   aGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAg           <  1:5066324/93‑1 (MQ=18)

                                                                    ggCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGc          <  1:1032522/93‑1 (MQ=18)

                                                                      cGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCaa        <  1:8424940/93‑1 (MQ=14)

ATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTC  >  NC_000913/1632051‑1632219

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