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.3%	378.2 / ‑4.7	135	S166P (TCA→CCA)	b1545	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (68/66); total (68/67)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 137 total aligned reads displayed.

ATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACACC  >  NC_000913/1632051‑1632226

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

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

   ctAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTCTTAAATGAAGAACAGAAACAGGCGGTATTCCAACGAATTAAGCCAGgt                                                                                  <  2:393496/92‑1 (MQ=11)

    tAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAggcg                                                                                 >  1:2887680/1‑91 (MQ=255)

    tAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtg                                                                                 >  1:5922553/1‑93 (MQ=37)

    tAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtg                                                                                 >  2:7905287/1‑93 (MQ=37)

    tAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtg                                                                                 <  1:2199222/93‑1 (MQ=255)

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

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

      agagCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGGGTa                                                                               >  2:806781/1‑93 (MQ=255)

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

        agCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAg                                                                             <  2:6290474/93‑1 (MQ=255)

        agCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAg                                                                             >  1:4965415/1‑93 (MQ=37)

          cccGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGta                                                                           <  1:916347/93‑1 (MQ=255)

          cccGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGta                                                                           <  1:5955360/93‑1 (MQ=255)

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

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

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

               ggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGt                                                                      <  1:4026645/93‑1 (MQ=255)

               ggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGt                                                                      <  1:764060/93‑1 (MQ=255)

                 cgcAGGGAAACGTTTTGGTCGACCACCTGTGTTAAGTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGc                                                                    <  1:3513739/93‑1 (MQ=255)

                 cgcAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGc                                                                    <  1:1867363/93‑1 (MQ=255)

                  gcAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGcc                                                                   >  1:1857344/1‑93 (MQ=255)

                  gcAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGcc                                                                   <  2:6339900/93‑1 (MQ=255)

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

                   cAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCa                                                                  >  1:7365650/1‑93 (MQ=255)

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

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

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

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

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

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

                         aaCGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGccc                                                            >  1:58462/1‑93 (MQ=255)

                            gTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCggg                                                         >  1:5165117/1‑93 (MQ=255)

                            gTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCggg                                                         >  1:6782783/1‑93 (MQ=255)

                            gTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCggg                                                         >  1:3647591/1‑93 (MQ=255)

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

                             ttttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGa                                                        <  1:3045186/93‑1 (MQ=255)

                             ttttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCAATTGCCCGGGa                                                        >  1:2741986/1‑93 (MQ=255)

                             ttttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGGATTCGAACGAATTAAGCCAGGTGTAAGTATAAGGGCCATTGCCCGGGa                                                        >  2:5308424/1‑93 (MQ=255)

                               ttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAt                                                      >  2:1677174/1‑93 (MQ=255)

                                tGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAtt                                                     <  1:3934723/93‑1 (MQ=255)

                                tGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAtt                                                     <  1:5412440/93‑1 (MQ=37)

                                   tCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaa                                                  <  2:6840268/93‑1 (MQ=37)

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

                                     gACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTACCCGGGAATTCaaaa                                                >  1:4371866/1‑93 (MQ=255)

                                      accaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                               <  2:5379959/93‑1 (MQ=37)

                                      accaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                               >  2:6666177/1‑93 (MQ=37)

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

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

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

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

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

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

                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTc                                            <  2:3620807/93‑1 (MQ=21)

                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTc                                            <  2:6814151/93‑1 (MQ=21)

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

                                          ccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcg                                           <  1:7403880/93‑1 (MQ=18)

                                            tgtgTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgcg                                         <  1:2241952/93‑1 (MQ=18)

                                            tgtgTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgcg                                         >  2:2912596/1‑93 (MQ=21)

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

                                               gtTAAATGAAGAACAGAAACAGGCAGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCa                                      <  1:7034018/93‑1 (MQ=16)

                                                ttAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCAGGGAATTCAAAACCTCGCGGCaa                                     >  2:5388792/1‑93 (MQ=255)

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

                                                  aaaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                                   <  2:6517700/93‑1 (MQ=18)

                                                  aaaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                                   <  1:2211106/93‑1 (MQ=18)

                                                  aaaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                                   >  1:4506850/1‑93 (MQ=21)

                                                  aaaTGAAGAAAAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                                   <  1:5829557/93‑1 (MQ=16)

                                                    aTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCa                                 >  2:3000799/1‑93 (MQ=21)

                                                     tGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCTCGGGAATTCAAAACCTCGCGGCAAACCAt                                <  2:5878242/93‑1 (MQ=255)

                                                     tGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAt                                >  1:1869259/1‑93 (MQ=21)

                                                     tGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAt                                >  2:1348049/1‑93 (MQ=21)

                                                     tGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAt                                <  2:5878243/93‑1 (MQ=18)

                                                      gaagaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAtt                               >  1:169486/1‑93 (MQ=21)

                                                      gaagaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAtt                               >  1:2566301/1‑93 (MQ=21)

                                                      gaagaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAtt                               <  1:7259095/93‑1 (MQ=18)

                                                        agaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGACTTCAAAACCTCGCGGCAAACCAtttt                             >  1:2982624/1‑93 (MQ=255)

                                                        agaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAtttt                             <  1:7733805/93‑1 (MQ=18)

                                                        agaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAtttt                             >  2:427192/1‑93 (MQ=21)

                                                         gaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTa                            >  1:5046744/1‑93 (MQ=21)

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

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

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

                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCACAACCTCGCGGCAAACCATTTTAaga                         >  2:4896089/1‑93 (MQ=16)

                                                            cAGAAACAGACGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga                         <  1:4419025/93‑1 (MQ=16)

                                                             aGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAagag                        >  1:3896265/1‑93 (MQ=21)

                                                             aGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAagag                        >  1:3384909/1‑93 (MQ=21)

                                                             aGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAagag                        >  1:3896264/1‑93 (MQ=21)

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

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

                                                                 aCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaa                    >  1:6244885/1‑93 (MQ=18)

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

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

                                                                  cAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaaa                   >  1:5132211/1‑93 (MQ=18)

                                                                   aGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAg                  >  2:6576013/1‑93 (MQ=18)

                                                                    ggCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGc                 >  2:2417132/1‑93 (MQ=18)

                                                                     gCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCa                >  1:3683165/1‑93 (MQ=17)

                                                                      cGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCaa               >  2:2080531/1‑93 (MQ=14)

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

                                                                                   aactAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACAcc  <  1:3615164/90‑1 (MQ=1)

ATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACACC  >  NC_000913/1632051‑1632226

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