BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,223,632	A→G	L15P (CTG→CCG)	yohC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,223,632	0	A	G	99.6%	733.3 / ‑5.0	236	L15P (CTG→CCG)	yohC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (119/116); total (120/116)
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.91e-01

Only 100 of 240 total aligned reads displayed.

GCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTT  >  NC_000913/2223542‑2223721

gccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGt                                                                                         >  1:7928312/1‑93 (MQ=255)

gccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGt                                                                                         <  2:2627304/93‑1 (MQ=255)

gccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGt                                                                                         >  1:6552002/1‑93 (MQ=255)

gccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGt                                                                                         <  1:1802796/93‑1 (MQ=255)

gccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGt                                                                                         >  1:6340527/1‑93 (MQ=255)

 ccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTc                                                                                        >  1:3122702/1‑93 (MQ=255)

 ccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTc                                                                                        <  2:4445296/93‑1 (MQ=255)

  cATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTcc                                                                                       >  1:5249534/1‑93 (MQ=255)

  cATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTcc                                                                                       <  1:272684/93‑1 (MQ=255)

   aTCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTccc                                                                                      >  1:3096128/1‑93 (MQ=255)

    tCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTcccc                                                                                     >  1:7249380/1‑93 (MQ=255)

    tCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTcccc                                                                                     >  1:7660306/1‑93 (MQ=255)

    tCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGCcccc                                                                                     >  1:5733998/1‑93 (MQ=255)

     cAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCa                                                                                    <  2:6592322/93‑1 (MQ=255)

     cAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCa                                                                                    <  1:1650863/93‑1 (MQ=255)

         aaagCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACa                                                                                <  1:4410095/93‑1 (MQ=255)

         aaaaCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACa                                                                                >  2:3373447/1‑93 (MQ=255)

             cGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGc                                                                            <  2:5579020/93‑1 (MQ=255)

              gtggtgGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGcc                                                                           >  1:2281447/1‑92 (MQ=255)

               tggtggGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTc                                                                          >  1:401498/1‑93 (MQ=255)

                ggtgggtgTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCa                                                                         <  2:7384063/93‑1 (MQ=255)

                 gtgggtgTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCAt                                                                        <  2:7628792/93‑1 (MQ=255)

                  tgggtgTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATa                                                                       <  1:7825905/93‑1 (MQ=255)

                  tgggtgTAATGGTGAGAAATCGTTTAGTTTTCGTGATTAATCACTTTCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATa                                                                       <  2:6802599/93‑1 (MQ=255)

                   gggtgTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATAc                                                                      >  2:3363000/1‑93 (MQ=255)

                     gtgtAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCAtcacg                                                                    >  1:2398047/1‑89 (MQ=255)

                      tgtAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGc                                                                   <  2:7004993/93‑1 (MQ=255)

                       gtAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGcc                                                                  >  2:831112/1‑93 (MQ=255)

                        tAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGccc                                                                 >  1:3473292/1‑93 (MQ=255)

                        tAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGccc                                                                 >  1:7437735/1‑93 (MQ=255)

                           tGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATAACGCCCCCt                                                              >  1:4761989/1‑93 (MQ=255)

                             gTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCttg                                                            >  1:5943524/1‑93 (MQ=255)

                               gagaAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCttgtt                                                          >  2:5975615/1‑93 (MQ=255)

                                 gaAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTtgttgt                                                        <  2:489234/93‑1 (MQ=255)

                                 gaAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTtgttgt                                                        >  2:2438694/1‑93 (MQ=255)

                                   aaTCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGa                                                      <  1:2012547/93‑1 (MQ=255)

                                   aaTCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGa                                                      >  1:1850596/1‑93 (MQ=255)

                                    aTCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAg                                                     <  2:6633535/93‑1 (MQ=255)

                                     tcgtttcgttTTCGCGATTAATCACCTGCATTTCACGATCGGCATGGAAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGt                                                    <  2:3503408/93‑1 (MQ=255)

                                       gtttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTc                                                  <  1:1225483/93‑1 (MQ=255)

                                        tttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCa                                                 >  2:6089711/1‑93 (MQ=255)

                                        tttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCa                                                 >  2:2767803/1‑93 (MQ=255)

                                         ttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCAt                                                >  2:6740484/1‑93 (MQ=255)

                                          tcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATg                                               <  1:3573431/93‑1 (MQ=255)

                                           cgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGa                                              <  2:1870544/93‑1 (MQ=255)

                                           cgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGa                                              <  1:4246135/93‑1 (MQ=255)

                                            gtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGaa                                             <  2:1558443/93‑1 (MQ=255)

                                            gtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGaa                                             <  2:3772636/93‑1 (MQ=255)

                                             ttttCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAAc                                            >  2:500882/1‑93 (MQ=255)

                                             ttttCGAGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAAc                                            <  1:5589207/93‑1 (MQ=255)

                                                  gcgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGAc                                       >  2:2579139/1‑93 (MQ=255)

                                                   cgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACa                                      <  2:367426/93‑1 (MQ=255)

                                                    gATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAg                                     <  1:4968804/93‑1 (MQ=255)

                                                    gATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGGCCCCAGACATGGCTCATACCGGCCCCTTGTTGTGAGTTCATGAACCTGACAg                                     >  2:4252733/1‑93 (MQ=255)

                                                     aTTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGt                                    >  1:7160851/1‑93 (MQ=255)

                                                        aaTCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCa                                 >  1:5269392/1‑93 (MQ=255)

                                                        aaTCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCa                                 >  1:2849798/1‑93 (MQ=255)

                                                        aaTCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCa                                 <  1:6656072/93‑1 (MQ=255)

                                                          tCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAg                               >  2:4821291/1‑93 (MQ=255)

                                                          tCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAg                               >  1:5104129/1‑93 (MQ=255)

                                                          tCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAg                               >  2:5747851/1‑93 (MQ=255)

                                                            cccTGCATTTCACGATCGGTATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGta                             <  1:2190234/92‑1 (MQ=255)

                                                            aCCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGta                             >  1:4301125/1‑93 (MQ=255)

                                                            aCCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGta                             >  1:26627/1‑93 (MQ=255)

                                                              cTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGtata                           >  2:7867863/1‑93 (MQ=255)

                                                              cTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGtata                           <  2:2850428/93‑1 (MQ=255)

                                                              cTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGtata                           >  1:2919472/1‑93 (MQ=255)

                                                               tGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATaa                          <  1:7149118/93‑1 (MQ=255)

                                                                gCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATCCCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAg                         <  2:3075414/93‑1 (MQ=255)

                                                                gCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAg                         <  1:5447560/93‑1 (MQ=255)

                                                                gCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAg                         >  2:1665306/1‑93 (MQ=255)

                                                                gCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAg                         >  1:661221/1‑93 (MQ=255)

                                                                  aTTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTc                       <  2:4710708/93‑1 (MQ=255)

                                                                   tttCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCa                      >  2:630207/1‑93 (MQ=255)

                                                                    ttCACGATCTGGATGGGAAAACGGTCCCCAGACCTGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAg                     <  1:1369750/93‑1 (MQ=255)

                                                                     tCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGc                    >  1:1028686/1‑93 (MQ=255)

                                                                      cACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCt                   >  1:6079808/1‑93 (MQ=255)

                                                                      cACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCt                   >  2:5314766/1‑93 (MQ=255)

                                                                       aCGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCtt                  >  2:1399711/1‑93 (MQ=255)

                                                                          aTCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtgtg               >  1:5480540/1‑93 (MQ=255)

                                                                           tCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGa              >  1:4574639/1‑93 (MQ=255)

                                                                            cGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGat             <  2:2189459/93‑1 (MQ=255)

                                                                            cGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGat             <  2:7310652/93‑1 (MQ=255)

                                                                            cGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGat             <  2:6918618/93‑1 (MQ=255)

                                                                             gggATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGatt            >  2:272590/1‑93 (MQ=255)

                                                                               gATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGattat          >  2:6504342/1‑93 (MQ=255)

                                                                               gATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGattat          >  2:7546941/1‑93 (MQ=255)

                                                                               gATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTCCAAGTATAAGTCAGCTTGTGattat          >  1:880928/1‑93 (MQ=255)

                                                                                aTGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGattatt         >  2:1185169/1‑93 (MQ=255)

                                                                                 tGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAttt        <  1:4663713/93‑1 (MQ=255)

                                                                                 tGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAttt        <  1:1435925/93‑1 (MQ=255)

                                                                                    gAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAtttttt     <  2:7526537/93‑1 (MQ=255)

                                                                                    gAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAtttttt     <  2:3984466/93‑1 (MQ=255)

                                                                                     aaaaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTg    <  1:6898549/93‑1 (MQ=255)

                                                                                     aaaaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTg    <  1:6851478/93‑1 (MQ=255)

                                                                                     aaaaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTg    <  1:3793632/93‑1 (MQ=255)

                                                                                      aaaCGGGCCCCAGGCATGGCCCATTCCGGCCCCTTGGTTTGAGTTTATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGt   <  2:67047/93‑1 (MQ=255)

                                                                                       aaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGtt  >  2:131084/1‑93 (MQ=255)

                                                                                       aaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGtt  >  1:6444541/1‑93 (MQ=255)

                                                                                       aaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGtt  >  1:3040572/1‑93 (MQ=255)

GCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTT  >  NC_000913/2223542‑2223721

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