BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,876,794	G→A	G318G (GGC→GGT)	gyrB ←	DNA gyrase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,876,794	0	G	A	99.6%	878.6 / NA	251	G318G (GGC→GGT)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (145/105); minor base C (0/1); total (145/106)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

CGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCA  >  NC_000913/3876704‑3876882

 gccgATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCt                                                                                       >  1:1789505/1‑93 (MQ=255)

 gccgATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCt                                                                                       <  1:337157/93‑1 (MQ=255)

 gccgATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCt                                                                                       >  2:4022795/1‑93 (MQ=255)

  ccgATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCtt                                                                                      <  2:2827747/93‑1 (MQ=255)

    gATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCa                                                                                    >  1:1035148/1‑93 (MQ=255)

     aTTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGGATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCAc                                                                                   >  2:3147122/1‑93 (MQ=255)

     aTTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCAc                                                                                   <  1:1371712/93‑1 (MQ=255)

     aTTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCAc                                                                                   >  2:3800168/1‑93 (MQ=255)

     aTTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCAc                                                                                   <  1:910020/93‑1 (MQ=255)

      tttCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCAcg                                                                                  >  1:1516560/1‑93 (MQ=255)

      tttCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCAcg                                                                                  >  1:1240109/1‑93 (MQ=255)

       ttCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCAcgc                                                                                 >  1:1560613/1‑93 (MQ=255)

        tCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTATCACGGAAACGACCGCAATCAGACCTTCAcgcg                                                                                >  1:3768929/1‑93 (MQ=255)

           ccTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCAt                                                                             >  1:2836016/1‑93 (MQ=255)

            cTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGCAACGACCGCAATCAGACCTTCACGCGCATc                                                                            >  1:2134940/1‑93 (MQ=255)

            cTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATc                                                                            >  2:1631437/1‑93 (MQ=255)

             tCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCg                                                                           >  2:1754447/1‑93 (MQ=255)

             tCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCg                                                                           >  2:4626599/1‑93 (MQ=255)

              cAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGt                                                                          <  2:4428819/93‑1 (MQ=255)

              cAGAAGAAACAAGTTTGTCTTTGGTCTGGGAGGAGAGTTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGt                                                                          <  1:997923/93‑1 (MQ=255)

                gaagaaACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCa                                                                        <  2:1151661/93‑1 (MQ=255)

                 aagaaACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCAc                                                                       <  1:2704480/93‑1 (MQ=255)

                  agaaACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCAcc                                                                      >  2:3506243/1‑93 (MQ=255)

                   gaaACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCg                                                                     >  1:410260/1‑93 (MQ=255)

                   gaaACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCg                                                                     <  1:449008/93‑1 (MQ=255)

                    aaaCCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCgg                                                                    <  1:4673001/93‑1 (MQ=255)

                    aaaCCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCgg                                                                    >  2:1397141/1‑93 (MQ=255)

                     aaCCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCggg                                                                   >  1:3873464/1‑92 (MQ=255)

                       ccAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTgg                                                                 >  1:3086677/1‑93 (MQ=255)

                        cAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGgc                                                                >  2:2478283/1‑93 (MQ=255)

                         aGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGgcg                                                               >  1:3534180/1‑93 (MQ=255)

                         aGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGgcg                                                               >  2:2907334/1‑93 (MQ=255)

                          gTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGgcgc                                                              >  2:3544941/1‑93 (MQ=255)

                           tttGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCt                                                             >  2:2676434/1‑93 (MQ=255)

                           tttGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCt                                                             >  2:2789658/1‑93 (MQ=255)

                           tttGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCt                                                             >  2:4370978/1‑93 (MQ=255)

                            ttGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTg                                                            >  2:2107350/1‑93 (MQ=255)

                             tGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGa                                                           <  2:2176708/93‑1 (MQ=255)

                              gTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGAc                                                          <  2:717588/93‑1 (MQ=255)

                              gTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGAc                                                          >  1:3255548/1‑93 (MQ=255)

                               tCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACt                                                         <  2:256611/93‑1 (MQ=255)

                               tCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACt                                                         <  1:4300828/93‑1 (MQ=255)

                               tCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACt                                                         <  1:1008791/93‑1 (MQ=255)

                                cTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACtt                                                        <  2:830105/93‑1 (MQ=255)

                                 tttGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACttt                                                       >  1:2672711/1‑93 (MQ=255)

                                  ttGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTg                                                      <  1:4082501/93‑1 (MQ=255)

                                  ttGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTg                                                      <  2:4486205/93‑1 (MQ=255)

                                     gTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGGCCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCt                                                   >  1:2514363/1‑93 (MQ=255)

                                     gTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCt                                                   <  2:4107723/93‑1 (MQ=255)

                                      tCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCtt                                                  >  1:32845/1‑93 (MQ=255)

                                      tCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCtt                                                  <  1:3971658/93‑1 (MQ=255)

                                        tGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCtttt                                                >  1:2314872/1‑93 (MQ=255)

                                        tGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCtttt                                                >  2:2011368/1‑93 (MQ=255)

                                        tGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCtttt                                                >  2:3211084/1‑93 (MQ=255)

                                         gggAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCttttt                                               <  1:4626291/93‑1 (MQ=255)

                                          ggaggaGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCtttttt                                              >  2:1206949/1‑93 (MQ=255)

                                            aggagAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGCCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGc                                            <  1:895790/93‑1 (MQ=255)

                                             ggagAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCt                                           >  2:3992177/1‑93 (MQ=255)

                                              gagaATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTg                                          >  2:4448678/1‑93 (MQ=255)

                                                gaATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTa                                        >  2:3346720/1‑93 (MQ=255)

                                                 aaTTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                       >  2:4274545/1‑93 (MQ=255)

                                                  aTTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGc                                      >  2:199710/1‑93 (MQ=255)

                                                  aTTTCGGGTCCGGCACTTTCACGGAAACGACCGCAAGCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGc                                      <  2:614847/93‑1 (MQ=255)

                                                    ttCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCct                                    <  2:424005/93‑1 (MQ=255)

                                                      cGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttc                                  <  1:938224/93‑1 (MQ=255)

                                                       gggTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttct                                 <  1:1106959/93‑1 (MQ=255)

                                                       cggTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttct                                 <  2:4674433/92‑1 (MQ=255)

                                                        ggTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttctt                                <  2:1454667/93‑1 (MQ=255)

                                                        ggTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttctt                                <  2:4049224/93‑1 (MQ=255)

                                                         gTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCttt                               <  1:1156895/93‑1 (MQ=255)

                                                         gTCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCttt                               <  2:3667390/93‑1 (MQ=255)

                                                          tCCGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTg                              <  2:2344558/93‑1 (MQ=255)

                                                           ccGGCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGt                             >  1:4219670/1‑93 (MQ=255)

                                                             ggCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTcc                           <  2:668407/93‑1 (MQ=255)

                                                              gCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCa                          <  1:3492970/93‑1 (MQ=255)

                                                              gCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCa                          <  2:2611722/93‑1 (MQ=255)

                                                              gCACTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCa                          >  1:320533/1‑93 (MQ=255)

                                                                aCTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATg                        >  1:2909588/1‑93 (MQ=255)

                                                                 cTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGt                       >  1:347437/1‑93 (MQ=255)

                                                                 cTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGt                       >  1:630789/1‑93 (MQ=255)

                                                                 cTTTCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGt                       >  2:2721944/1‑93 (MQ=255)

                                                                  cgtCACGGAAACGACCGCAATCGGCGCTTCACGCGCATCTTCACAGGTGGACATGAGTTTGACTTTTGTGCTGTAGCCTTCTTTGTCCATGTa                      <  2:3676954/91‑1 (MQ=255)

                                                                   ttCACGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAg                     <  1:778958/93‑1 (MQ=255)

                                                                       cGGAAACGACCGCAATCAGACCTTCACGCGCATCGTCTCCGGTGGCGCTGACTTTGGCATTTTTGCTGTAGCCTTGTTTGTCCATGTAGGCGt                 >  2:3822558/1‑93 (MQ=255)

                                                                        ggAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                >  2:967968/1‑93 (MQ=255)

                                                                         gAAACGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTc               >  2:1313926/1‑93 (MQ=255)

                                                                          aaaCGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCa              >  1:1073063/1‑93 (MQ=255)

                                                                           aaCGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAg             >  1:3780229/1‑93 (MQ=255)

                                                                           aaCGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAg             >  1:4561484/1‑93 (MQ=255)

                                                                           aaCGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAg             >  1:3935464/1‑93 (MQ=255)

                                                                           aaCGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAg             >  1:3640003/1‑93 (MQ=255)

                                                                            aCGACCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAgg            <  1:3819290/93‑1 (MQ=255)

                                                                               aCCGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGtc         >  1:1998359/1‑92 (MQ=255)

                                                                                ccGCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTAc        >  2:482704/1‑93 (MQ=255)

                                                                                ccGCAATCAGACCTTCACGCGCATCGTCACCGGTGGAGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTAc        <  1:3006/93‑1 (MQ=255)

                                                                                  gCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACgg      >  2:3452006/1‑93 (MQ=255)

                                                                                  gCAATCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACgg      >  2:3611210/1‑93 (MQ=255)

                                                                                    aaTCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGt    >  1:2464830/1‑93 (MQ=255)

                                                                                     aTCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTc   >  1:1665444/1‑93 (MQ=255)

                                                                                      tCAGACCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCa  <  1:4378366/93‑1 (MQ=255)

CGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCA  >  NC_000913/3876704‑3876882

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