BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	731,579	C→T	P925L (CCG→CTG)	rhsC →	RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	731,579	0	C	T	99.1%	311.0 / ‑5.5	115	P925L (CCG→CTG)	rhsC	RhsC
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (46/68); total (46/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 121 total aligned reads displayed.

TGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGAC  >  NC_000913/731489‑731655

tGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTgg                                                                            >  2:3775528/1‑93 (MQ=255)

tGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTgg                                                                            <  2:2908522/93‑1 (MQ=37)

  gATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGAt                                                                          >  1:186120/1‑93 (MQ=37)

  gATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGAt                                                                          <  1:39644/93‑1 (MQ=255)

  gATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGAt                                                                          <  2:572876/93‑1 (MQ=37)

   atatCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATa                                                                         <  2:3303361/93‑1 (MQ=37)

   atatCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATa                                                                         <  2:287187/93‑1 (MQ=255)

    tttCCGCCTCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCCACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATaa                                                                        <  1:536511/91‑1 (MQ=11)

      tCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAAcc                                                                      <  2:4546927/93‑1 (MQ=37)

       ccGCATCCCGTATCCCACAGCCCGGGCAGGTAACCGCC‑GGCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGt                                                                    <  2:3096978/93‑1 (MQ=11)

          cATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTAt                                                                  >  2:2900892/1‑93 (MQ=255)

          cATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTAt                                                                  <  1:707723/93‑1 (MQ=37)

           aTCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATc                                                                 >  2:2845828/1‑93 (MQ=255)

           aTCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATc                                                                 >  1:3394708/1‑93 (MQ=37)

           aTCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATc                                                                 <  1:2907752/93‑1 (MQ=255)

            tCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCg                                                                >  1:631731/1‑93 (MQ=255)

             cccGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGc                                                               <  1:3989678/93‑1 (MQ=37)

              ccGTATCCCACAGGCCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGcc                                                              <  2:3454454/93‑1 (MQ=255)

              ccGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGcc                                                              >  1:944134/1‑93 (MQ=255)

              ccGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGcc                                                              <  1:129664/93‑1 (MQ=37)

              ccGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGcc                                                              >  2:3003889/1‑93 (MQ=255)

                gTATACCACAGACCCGGCGGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCg                                                            >  1:2270680/1‑93 (MQ=255)

                gTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCg                                                            <  2:857729/93‑1 (MQ=255)

                gTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCg                                                            >  2:1907786/1‑93 (MQ=37)

                 tataCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGt                                                           >  1:1647103/1‑93 (MQ=37)

                 tataCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGt                                                           <  2:1899873/93‑1 (MQ=37)

                 tataCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGt                                                           <  2:1420590/93‑1 (MQ=37)

                 tataCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGt                                                           >  1:4602599/1‑93 (MQ=37)

                  ataCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTg                                                          >  2:3594719/1‑93 (MQ=37)

                   taCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGa                                                         <  1:1732926/93‑1 (MQ=255)

                   taCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGa                                                         <  1:3730562/93‑1 (MQ=255)

                   taCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGGGa                                                         >  1:3635171/1‑93 (MQ=255)

                    aCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGAc                                                        >  1:1928127/1‑93 (MQ=21)

                    aCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGAc                                                        >  2:3136651/1‑93 (MQ=21)

                      cacaGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGAcgc                                                      >  2:3184106/1‑93 (MQ=18)

                      cacaGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGAcgc                                                      >  2:4189079/1‑93 (MQ=18)

                          gACCCTGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCAc                                                  >  1:2455436/1‑93 (MQ=255)

                          gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCAc                                                  <  1:2847665/93‑1 (MQ=18)

                          gACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCAc                                                  >  1:2522782/1‑93 (MQ=18)

                            cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTa                                                <  2:1662284/93‑1 (MQ=18)

                             ccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTAt                                               >  1:612298/1‑93 (MQ=18)

                                gCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCtt                                            <  2:3232478/93‑1 (MQ=18)

                                 cAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCttt                                           <  2:2472763/93‑1 (MQ=18)

                                 cAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCttt                                           <  1:91478/93‑1 (MQ=18)

                                 cAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGAGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCttt                                           >  2:3763643/1‑93 (MQ=16)

                                   ggTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAc                                         <  1:1681975/93‑1 (MQ=18)

                                   ggTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAc                                         <  1:1714711/93‑1 (MQ=18)

                                    gTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAcc                                        <  2:3819029/93‑1 (MQ=18)

                                    gTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAcc                                        <  2:782963/93‑1 (MQ=18)

                                    gTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAcc                                        <  2:3130986/93‑1 (MQ=18)

                                     tAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCg                                       <  2:3486475/93‑1 (MQ=18)

                                     tAACCGCCTGCCCGACCCGGAGCGGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCg                                       <  2:2211579/93‑1 (MQ=16)

                                      agccGCCTGCCAGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTAGATACCCGTATCGCCCGTGACGCGCACTATCTTTACCgg                                      <  2:2648858/91‑1 (MQ=11)

                                       aCCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGt                                     >  2:2870784/1‑93 (MQ=18)

                                          gCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATg                                  >  1:3780183/1‑93 (MQ=18)

                                          gCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATg                                  >  2:2140402/1‑93 (MQ=18)

                                           ccTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGa                                 <  2:429532/93‑1 (MQ=18)

                                           ccTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGa                                 <  2:821262/93‑1 (MQ=18)

                                            cTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGAc                                <  2:4043174/93‑1 (MQ=18)

                                            cTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGAc                                <  1:2603086/93‑1 (MQ=18)

                                            cTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGAc                                <  2:2128853/93‑1 (MQ=18)

                                            cTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGAc                                <  2:2128854/93‑1 (MQ=18)

                                             tGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGAcc                               >  2:1854565/1‑93 (MQ=18)

                                             tGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGAcc                               >  1:2585246/1‑93 (MQ=18)

                                             tGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGAcc                               <  1:2930255/93‑1 (MQ=18)

                                               cccGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGt                             <  2:606002/93‑1 (MQ=18)

                                               cccGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGt                             >  2:3610048/1‑93 (MQ=18)

                                                ccGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTc                            >  1:3198948/1‑93 (MQ=21)

                                                 cGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCa                           <  2:37899/93‑1 (MQ=18)

                                                 cGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCa                           <  1:3110435/93‑1 (MQ=18)

                                                 cGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCa                           <  2:3307391/93‑1 (MQ=18)

                                                  gACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCAc                          <  2:4197091/93‑1 (MQ=18)

                                                   aCCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACg                         >  2:2272531/1‑93 (MQ=21)

                                                   aCCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACg                         >  1:2992845/1‑93 (MQ=21)

                                                       ggAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAg                     <  1:1515388/93‑1 (MQ=18)

                                                       ggAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAg                     <  1:3972738/93‑1 (MQ=18)

                                                         aGCTGCCCCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGc                   <  2:2259248/93‑1 (MQ=16)

                                                         aGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGc                   <  1:2581508/93‑1 (MQ=18)

                                                         aGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGc                   <  1:1782741/93‑1 (MQ=18)

                                                         aGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGc                   >  2:3030887/1‑93 (MQ=18)

                                                         aGCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGc                   <  1:1316706/93‑1 (MQ=18)

                                                          gCTGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt                  <  1:4329211/93‑1 (MQ=18)

                                                            tGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGa                <  1:299914/93‑1 (MQ=18)

                                                            tGCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGa                >  2:1741905/1‑93 (MQ=18)

                                                            tGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGa                <  1:2883523/93‑1 (MQ=33)

                                                             gCACCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGAc               <  1:177197/93‑1 (MQ=18)

                                                               aCCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACag             >  2:468952/1‑93 (MQ=18)

                                                               aCCCGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACag             >  1:1249477/1‑93 (MQ=18)

                                                                cccGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCTCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACaga            <  1:4395759/93‑1 (MQ=16)

                                                                cccGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACaga            <  1:3200660/93‑1 (MQ=18)

                                                                cccGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACaga            <  1:1321704/93‑1 (MQ=18)

                                                                 ccGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACagag           >  2:3235264/1‑93 (MQ=21)

                                                                 ccGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACagag           <  2:593552/93‑1 (MQ=18)

                                                                  cGGACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAgaga          <  2:3705022/93‑1 (MQ=18)

                                                                   ggACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaa         <  1:3924881/93‑1 (MQ=18)

                                                                   ggACAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaa         <  1:903499/93‑1 (MQ=18)

                                                                      cAGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAAc      <  1:1178535/93‑1 (MQ=18)

                                                                       aGCGCCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAAcc     >  2:731583/1‑93 (MQ=18)

                                                                         cgcCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGa   <  2:2133755/93‑1 (MQ=17)

                                                                          gcCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGAc  <  2:4677761/93‑1 (MQ=14)

TGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGAC  >  NC_000913/731489‑731655

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