BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,676,843	T→C	R131R (CGT→CGC)	b1604 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,676,843	0	T	C	99.5%	605.3 / ‑5.0	187	R131R (CGT→CGC)	b1604	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (88/98); total (88/99)
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.73e-01

Only 100 of 195 total aligned reads displayed.

TCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAA  >  NC_000913/1676753‑1676932

  gcgTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagcc                                                                                       <  1:1026608/93‑1 (MQ=255)

  gcgTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagcc                                                                                       <  1:1142534/93‑1 (MQ=255)

    gTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagccag                                                                                     >  2:5548407/1‑93 (MQ=255)

    gTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagccag                                                                                     >  2:3920173/1‑93 (MQ=255)

     tAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagcaagc                                                                                    >  2:4464718/1‑93 (MQ=255)

      aaTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagccagcc                                                                                   <  2:6553938/93‑1 (MQ=255)

      aaTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagccagcc                                                                                   <  1:2425232/93‑1 (MQ=255)

      aaTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagccagcc                                                                                   <  1:1678263/93‑1 (MQ=255)

      aaTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCagccagcc                                                                                   <  2:1628705/93‑1 (MQ=255)

       aTTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTAATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAgccagcca                                                                                  <  1:6244899/93‑1 (MQ=255)

        ttAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGccagccag                                                                                 >  2:3205862/1‑93 (MQ=255)

        ttAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGccagccag                                                                                 >  1:1012870/1‑93 (MQ=255)

         tAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGa                                                                                >  1:2701760/1‑93 (MQ=255)

          aaCGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGaa                                                                               <  1:2233450/93‑1 (MQ=255)

           aCGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAg                                                                              >  1:5772122/1‑93 (MQ=255)

            cGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGt                                                                             <  1:1418117/93‑1 (MQ=255)

            cGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGt                                                                             <  1:5709506/93‑1 (MQ=255)

              gtgtTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCa                                                                           >  2:5859864/1‑93 (MQ=255)

                gttgtCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAat                                                                         >  2:4137703/1‑93 (MQ=255)

                 ttgtCGAACTTCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAatg                                                                        <  1:5507096/93‑1 (MQ=255)

                 ttgtCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAatg                                                                        >  1:5465958/1‑93 (MQ=255)

                  tgtCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAatga                                                                       <  2:138210/93‑1 (MQ=255)

                   gtCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAatgat                                                                      >  2:4820413/1‑93 (MQ=255)

                     cGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGc                                                                    <  2:2815204/93‑1 (MQ=255)

                       aaCTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCa                                                                  <  1:6314422/93‑1 (MQ=255)

                        aCTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCAt                                                                 >  2:718122/1‑93 (MQ=255)

                        aCTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCAt                                                                 >  2:5397216/1‑93 (MQ=255)

                        aCTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCAt                                                                 <  1:6602797/93‑1 (MQ=255)

                         cTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATc                                                                >  2:1966590/1‑93 (MQ=255)

                         cTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATc                                                                >  1:1062056/1‑93 (MQ=255)

                           gCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCAc                                                              <  1:6444245/93‑1 (MQ=255)

                           gCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCAc                                                              <  1:1692114/93‑1 (MQ=255)

                           gCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCAc                                                              <  1:4126581/93‑1 (MQ=255)

                            ccaaaaGATCAGGCTGTTCTGATTGACCCGTTTGACCCGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCAcc                                                             <  1:3454160/90‑1 (MQ=255)

                               aaaGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCaaa                                                          <  2:2127782/93‑1 (MQ=255)

                                 aGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGc                                                        <  1:4483667/93‑1 (MQ=255)

                                 aGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGc                                                        <  2:4751140/93‑1 (MQ=255)

                                  gATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCg                                                       >  1:339081/1‑93 (MQ=255)

                                  gATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCg                                                       >  2:332460/1‑93 (MQ=255)

                                   aTCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCgg                                                      >  2:2354302/1‑93 (MQ=255)

                                    tCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGc                                                     >  2:4202561/1‑93 (MQ=255)

                                    tCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGc                                                     >  2:6103331/1‑93 (MQ=255)

                                     cAGGCTGTTCTGATTGAACCGTTTGACACGGTCCCCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCa                                                    <  2:4581040/93‑1 (MQ=255)

                                      aGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                   >  2:2541460/1‑93 (MQ=255)

                                      aGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                   >  1:4624198/1‑93 (MQ=255)

                                      aGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                   >  1:2078085/1‑93 (MQ=255)

                                      aGGCTGTTCTGATTGAACCGTTTGACACGGGCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaa                                                   >  1:2209054/1‑93 (MQ=255)

                                       ggCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaa                                                  >  2:1469290/1‑93 (MQ=255)

                                        gCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaaa                                                 <  1:1838062/93‑1 (MQ=255)

                                        gCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCaaaa                                                 <  2:6268290/93‑1 (MQ=255)

                                         cTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAg                                                <  1:3091601/93‑1 (MQ=255)

                                         cTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCACCCATCCAGAAGTCAATGATGCCATCACCCAAGCGGCAAAAg                                                >  1:1390782/1‑93 (MQ=255)

                                          tGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGc                                               >  2:6436235/1‑93 (MQ=255)

                                          tGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGc                                               <  1:2678365/93‑1 (MQ=255)

                                           gTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCg                                              <  2:3609972/93‑1 (MQ=255)

                                           gTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCg                                              >  2:2985961/1‑93 (MQ=255)

                                            ttCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGa                                             >  2:4000238/1‑93 (MQ=255)

                                            ttCTGATTGAACCGCTTGACAAGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACTAAAGCGGCAAAAGCGa                                             >  1:2658032/1‑93 (MQ=255)

                                              cTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGaaa                                           <  1:2961257/93‑1 (MQ=255)

                                                gATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAgg                                         <  1:5881805/93‑1 (MQ=255)

                                                 aTTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAggg                                        >  2:5978864/1‑92 (MQ=255)

                                                  ttGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTg                                       >  1:4925568/1‑93 (MQ=255)

                                                    taaCCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGcc                                     <  2:6186005/92‑1 (MQ=255)

                                                    gAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGcc                                     >  2:5039526/1‑93 (MQ=255)

                                                       ccGTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTAc                                  >  2:3125006/1‑93 (MQ=255)

                                                         gTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTActc                                <  2:403199/93‑1 (MQ=255)

                                                         gTTTGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTActc                                <  2:434148/93‑1 (MQ=255)

                                                          tttGACGCGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTActct                               <  2:3045573/93‑1 (MQ=255)

                                                          tttGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTActct                               <  1:3168796/93‑1 (MQ=255)

                                                          tttGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTActct                               >  2:1949862/1‑93 (MQ=255)

                                                           ttGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGAGGCAAAAGCGAAAGGTGACTACTAtt                              >  1:3054048/1‑93 (MQ=255)

                                                            tGACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCttt                             <  1:2868569/93‑1 (MQ=255)

                                                             gACACGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTc                            <  1:4800040/93‑1 (MQ=255)

                                                              acacGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCt                           <  1:1887555/93‑1 (MQ=255)

                                                              acacGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCt                           <  1:2575422/93‑1 (MQ=255)

                                                              acacGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCt                           <  1:5196101/93‑1 (MQ=255)

                                                              acacGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCt                           <  2:5617177/93‑1 (MQ=255)

                                                              acacGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCt                           >  2:5409332/1‑93 (MQ=255)

                                                               cacGGTCACCGTCTAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTa                          >  1:330025/1‑93 (MQ=255)

                                                               cacGGTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTa                          <  1:5974186/93‑1 (MQ=255)

                                                                  ggTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAt                       >  2:4855655/1‑93 (MQ=255)

                                                                  ggTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAt                       <  1:6071911/93‑1 (MQ=255)

                                                                   gTCACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtc                      <  2:3806079/93‑1 (MQ=255)

                                                                     cACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgt                    >  1:3362180/1‑93 (MQ=255)

                                                                     cACCGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgt                    >  1:6278231/1‑93 (MQ=255)

                                                                        cGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgttcg                 <  1:4969811/93‑1 (MQ=255)

                                                                        cGTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgttcg                 <  2:5739752/93‑1 (MQ=255)

                                                                         gTCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACAtcgttcgt                >  2:5709443/1‑93 (MQ=255)

                                                                          tCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTc               <  2:4932168/93‑1 (MQ=255)

                                                                          tCCAGGGCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTc               <  1:1552975/93‑1 (MQ=255)

                                                                                gCTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCg         >  1:1654131/1‑93 (MQ=255)

                                                                                 cTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa        <  1:2477558/93‑1 (MQ=255)

                                                                                 cTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa        <  1:5348470/93‑1 (MQ=255)

                                                                                 cTTCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGa        <  1:929666/93‑1 (MQ=255)

                                                                                  ttCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGAt       >  2:2699600/1‑93 (MQ=255)

                                                                                   tCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATg      >  1:1643107/1‑93 (MQ=255)

                                                                                   tCTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACACCGTTCGTCAAATCGATg      <  2:944540/93‑1 (MQ=255)

                                                                                    cTATCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGc     <  2:257266/93‑1 (MQ=255)

                                                                                       tCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCaa  <  1:1967151/93‑1 (MQ=255)

                                                                                       tCGCAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAACGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCaa  <  1:5439122/93‑1 (MQ=255)

TCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAA  >  NC_000913/1676753‑1676932

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