BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,923,720	T→C	V86A (GTT→GCT)	b1843 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,923,720	0	T	C	99.6%	756.4 / ‑4.9	234	V86A (GTT→GCT)	b1843	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (119/114); total (120/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

Only 100 of 237 total aligned reads displayed.

CCTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCAT  >  NC_000913/1923631‑1923809

ccTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGa                                                                                        <  2:4261428/93‑1 (MQ=255)

 cTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGaa                                                                                       <  2:619649/93‑1 (MQ=255)

  ttCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAat                                                                                      >  1:9146353/1‑93 (MQ=255)

  ttCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAat                                                                                      >  1:819097/1‑93 (MQ=255)

   tCCTGCCCCACCCGATCTCTCCCTTTTAGACACGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCTTGCCTTCCCGCTGCata                                                                                     >  2:7261024/1‑93 (MQ=255)

   tCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAata                                                                                     >  2:560090/1‑93 (MQ=255)

    ccTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAatat                                                                                    >  2:6011656/1‑93 (MQ=255)

     cTGCCCCACCCGATCTCTCACTGTTAGCCCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAAtata                                                                                   <  1:4815245/93‑1 (MQ=255)

     cTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAAtata                                                                                   <  2:3097919/93‑1 (MQ=255)

     cTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAAtata                                                                                   <  2:659032/93‑1 (MQ=255)

      tGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATaa                                                                                  <  2:8965311/93‑1 (MQ=255)

      tGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATaa                                                                                  <  2:2450501/93‑1 (MQ=255)

      tGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATaa                                                                                  <  2:6430738/93‑1 (MQ=255)

      cgCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATaa                                                                                  <  2:1893848/92‑1 (MQ=255)

       ggcccgcccgCTCTCTCTCTGGTAGACCCTCTTTGTTATGCTGCAACGACCTGGCGAATGACCATCATTGCCGGCCTCCCCGCTGAATAGAAc                                                                                 <  2:2675819/91‑1 (MQ=255)

         cccACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGa                                                                               <  2:3881837/93‑1 (MQ=255)

          ccACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGAt                                                                              <  1:4161713/93‑1 (MQ=255)

          ccACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGAt                                                                              <  2:5676344/93‑1 (MQ=255)

           cACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATc                                                                             <  2:8639242/93‑1 (MQ=255)

           cACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTCCCGGCCTCCCCGCTGAATATAACGATc                                                                             <  2:8510114/93‑1 (MQ=255)

             cccGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGc                                                                           <  1:4555022/93‑1 (MQ=255)

             cccGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGc                                                                           >  1:6365631/1‑93 (MQ=255)

              ccGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCt                                                                          <  2:6688192/93‑1 (MQ=255)

              ccGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCt                                                                          >  1:5805573/1‑93 (MQ=255)

                gATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCttt                                                                        >  1:8032353/1‑93 (MQ=255)

                 aTCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCtttt                                                                       >  2:5093999/1‑92 (MQ=255)

                 aTCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTa                                                                       >  1:9750548/1‑93 (MQ=255)

                 aTCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTa                                                                       >  2:2175223/1‑93 (MQ=255)

                 aTCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTa                                                                       >  1:10026473/1‑93 (MQ=255)

                 aTCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTa                                                                       <  1:7645136/93‑1 (MQ=255)

                   ctctcACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTAtt                                                                     <  2:8308964/93‑1 (MQ=255)

                   ctctcACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTAtt                                                                     >  2:6769409/1‑93 (MQ=255)

                    tctcACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTc                                                                    >  2:1102987/1‑93 (MQ=255)

                    tctcACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTc                                                                    <  1:5923213/93‑1 (MQ=255)

                      tcACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGt                                                                  <  1:2442596/93‑1 (MQ=255)

                      tcACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGt                                                                  >  1:2127398/1‑93 (MQ=255)

                      tcACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGt                                                                  >  1:3724133/1‑93 (MQ=255)

                       cACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTg                                                                 >  1:7242500/1‑93 (MQ=255)

                        aCTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATCTAACGATCGCTTTATTCGTgg                                                                >  2:5306524/1‑93 (MQ=255)

                         cTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGa                                                               >  2:4690035/1‑93 (MQ=255)

                         cTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGa                                                               >  1:9854541/1‑93 (MQ=255)

                          tGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGaa                                                              >  2:2296484/1‑93 (MQ=255)

                           gTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGTAAt                                                             >  1:7440574/1‑93 (MQ=255)

                            ttAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCGTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAAtt                                                            <  2:4915731/93‑1 (MQ=255)

                             tAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTg                                                           <  2:1288239/93‑1 (MQ=255)

                                aCCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCgg                                                        <  2:2882419/93‑1 (MQ=255)

                                 cccGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGgt                                                       >  2:5938356/1‑93 (MQ=255)

                                  ccGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGgtg                                                      <  1:5538834/93‑1 (MQ=255)

                                  ccGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGgtg                                                      >  1:3417033/1‑93 (MQ=255)

                                     cTTTGTTATGCCGCAACGACCTGGCGAATGCCCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTc                                                   <  1:3040432/93‑1 (MQ=255)

                                      tttGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCg                                                  >  1:1747067/1‑93 (MQ=255)

                                        tGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCa                                                >  2:9513722/1‑93 (MQ=255)

                                         gTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCAc                                               >  1:264894/1‑93 (MQ=255)

                                          ttATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCAcc                                              >  1:6993454/1‑93 (MQ=255)

                                            aTGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGt                                            >  1:7104877/1‑93 (MQ=255)

                                            aTGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGt                                            <  1:8284629/93‑1 (MQ=255)

                                              gCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTgg                                          <  1:4381405/93‑1 (MQ=255)

                                               ccGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGc                                         <  2:7141436/93‑1 (MQ=255)

                                                cGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCg                                        >  1:9742003/1‑93 (MQ=255)

                                                 gCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGa                                       <  1:6588904/93‑1 (MQ=255)

                                                     cGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGaaaaa                                   >  1:8516710/1‑93 (MQ=255)

                                                     cGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGaaaaa                                   <  1:7562873/93‑1 (MQ=255)

                                                     cGACCTGGCGAATGACCATCATTCCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGaaaaa                                   <  2:5532304/93‑1 (MQ=255)

                                                       aCCTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACg                                 >  1:9581313/1‑93 (MQ=255)

                                                         cTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATTACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGcc                               <  2:2006607/93‑1 (MQ=255)

                                                         cTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGcc                               >  2:9701469/1‑93 (MQ=255)

                                                         cTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGcc                               >  1:288102/1‑93 (MQ=255)

                                                         cTGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGcc                               >  1:3904325/1‑93 (MQ=255)

                                                          tGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCt                              <  2:1706335/93‑1 (MQ=255)

                                                          tGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCt                              >  2:8254368/1‑93 (MQ=255)

                                                          tGGCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCt                              >  1:4436382/1‑93 (MQ=255)

                                                           ggCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTg                             <  1:2582176/93‑1 (MQ=255)

                                                           ggCGAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTg                             <  2:6946658/93‑1 (MQ=255)

                                                              gAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGa                          >  2:9177828/1‑93 (MQ=255)

                                                              gAATGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGa                          <  2:3604805/93‑1 (MQ=255)

                                                               aaTGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGAt                         <  1:1643061/93‑1 (MQ=255)

                                                                aTGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATc                        <  2:3999042/93‑1 (MQ=255)

                                                                aTGACCATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATc                        <  2:2613241/93‑1 (MQ=255)

                                                                   gccAGCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTAc                     <  2:1576390/92‑1 (MQ=255)

                                                                    ccATCATTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTAcc                    <  1:8328540/93‑1 (MQ=255)

                                                                     catcatTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCa                   >  1:3962883/1‑93 (MQ=255)

                                                                     catcatTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCa                   >  1:1058046/1‑93 (MQ=255)

                                                                       tcatTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATc                 <  1:5140122/93‑1 (MQ=255)

                                                                       tcatTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATc                 <  1:2977129/93‑1 (MQ=255)

                                                                       tcatTGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATc                 <  2:5168995/93‑1 (MQ=255)

                                                                          ttGCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCaga              <  2:420636/93‑1 (MQ=255)

                                                                            gCCGGCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGt            <  1:5366552/93‑1 (MQ=255)

                                                                               ggCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCAt         >  2:7974142/1‑93 (MQ=255)

                                                                               ggCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCAt         >  1:2241841/1‑93 (MQ=255)

                                                                               ggCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGGGTTCGCACCGCGGCGAAAAACGCCTGGGATCTACCATCAGAGTCAt         >  2:3939665/1‑93 (MQ=255)

                                                                                gCCTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATg        >  1:2023980/1‑93 (MQ=255)

                                                                                 ccTTCCCGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATgg       >  1:7975986/1‑93 (MQ=255)

                                                                                 ccTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATgg       <  2:7877933/93‑1 (MQ=255)

                                                                                 ccTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATgg       <  1:6202246/93‑1 (MQ=255)

                                                                                 ccTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATgg       <  1:1365862/93‑1 (MQ=255)

                                                                                  cTTCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGt      <  2:7547054/93‑1 (MQ=255)

                                                                                   ttCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTg     >  2:229953/1‑93 (MQ=255)

                                                                                   ttCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTg     >  1:3537629/1‑93 (MQ=255)

                                                                                    tCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGc    >  1:4639698/1‑93 (MQ=255)

                                                                                    tCCCGCTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGc    >  1:8770166/1‑93 (MQ=255)

CCTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCAT  >  NC_000913/1923631‑1923809

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: `‑`