BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,407,232	G→A	A176T (GCG→ACG)	prmA →	PrmA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,407,232	0	G	A	98.5%	691.0 / ‑2.8	198	A176T (GCG→ACG)	prmA	PrmA
Reads supporting (aligned to +/- strand): ref base G (0/3); new base A (86/109); total (86/112)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 214 total aligned reads displayed.

TCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGC  >  NC_000913/3407142‑3407321

tctctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCgcgacg                                                                                         >  1:1147064/1‑93 (MQ=255)

tctctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCgcgacg                                                                                         >  1:3537024/1‑93 (MQ=255)

 ctctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGacgc                                                                                        >  2:352304/1‑93 (MQ=255)

 ctctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGacgc                                                                                        <  2:2894309/93‑1 (MQ=255)

  tctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCt                                                                                       >  1:4321065/1‑93 (MQ=255)

  tctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCt                                                                                       <  2:2685948/93‑1 (MQ=255)

    tgtgCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGa                                                                                     >  1:4421973/1‑93 (MQ=255)

     gtgCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGag                                                                                    <  2:1231271/93‑2 (MQ=255)

     gtgCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCCCGACGCTgca                                                                                    >  1:2817155/1‑91 (MQ=255)

      tgCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGaaa                                                                                   <  2:1903204/93‑1 (MQ=255)

       gCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAAc                                                                                  <  2:4675609/93‑1 (MQ=255)

         cTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTg                                                                                <  2:614400/93‑1 (MQ=255)

         cTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTg                                                                                >  2:1840405/1‑93 (MQ=255)

         cTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTg                                                                                <  1:958662/93‑1 (MQ=255)

         cTGCAATGGCTCGACAGCCTCGATTTAACCGGGAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATAGCGACGCTGAAACTg                                                                                >  1:4054900/1‑93 (MQ=255)

          tGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGGTCCGGCATTCTGGCGATCGCGACGCTGAAACTgg                                                                               >  1:2793584/1‑93 (MQ=255)

           gCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTggg                                                                              >  1:1173540/1‑93 (MQ=255)

             aaTGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTg                                                                            >  1:4267154/1‑93 (MQ=255)

             aaTGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTg                                                                            >  2:3671226/1‑93 (MQ=255)

             aaTGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTg                                                                            >  2:3613440/1‑93 (MQ=255)

              aTGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTgc                                                                           >  2:746334/1‑93 (MQ=255)

               tGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTgca                                                                          <  1:1237361/93‑1 (MQ=255)

                   tCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCaa                                                                      >  1:3651971/1‑93 (MQ=255)

                     gACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCaaaa                                                                    >  2:2376342/1‑93 (MQ=255)

                      aCAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAg                                                                   <  2:3454856/93‑1 (MQ=255)

                      aCAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAg                                                                   >  1:2100179/1‑93 (MQ=255)

                        aGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGcc                                                                 >  2:3975162/1‑93 (MQ=255)

                        aGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGcc                                                                 <  1:2592285/93‑1 (MQ=255)

                           cTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCAtt                                                              >  2:4502755/1‑93 (MQ=255)

                           cTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCAtt                                                              >  2:328290/1‑93 (MQ=255)

                            tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTg                                                             >  2:43506/1‑93 (MQ=255)

                            tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTg                                                             >  2:3471515/1‑93 (MQ=255)

                            tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTg                                                             >  1:4672725/1‑93 (MQ=255)

                             cGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTgg                                                            <  2:2584815/93‑1 (MQ=255)

                              gATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGt                                                           <  2:3059172/93‑1 (MQ=255)

                                tttAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTAt                                                         >  1:859150/1‑93 (MQ=255)

                                   aaCCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGa                                                      <  1:3866621/93‑1 (MQ=255)

                                      cGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGatat                                                   >  2:3341678/1‑93 (MQ=255)

                                      cGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGatat                                                   <  2:2152255/93‑1 (MQ=255)

                                       ggTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATc                                                  <  1:4017789/93‑1 (MQ=255)

                                       ggTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATc                                                  <  2:799288/93‑1 (MQ=255)

                                        gTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCg                                                 <  1:1973275/93‑1 (MQ=255)

                                         tAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGa                                                <  1:3395849/93‑1 (MQ=255)

                                         tAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGa                                                <  1:3438790/93‑1 (MQ=255)

                                         gaaaaCAGTCATCGACTTTGGCTGTTGTTCCGGCATGCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGa                                                <  2:2454643/92‑1 (MQ=255)

                                          aaaaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGAt                                               <  2:4405347/93‑1 (MQ=255)

                                           aaaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATc                                              <  1:3795469/93‑1 (MQ=255)

                                            aaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATcc                                             <  2:4461395/93‑1 (MQ=255)

                                            aaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATcc                                             <  1:1399987/93‑1 (MQ=255)

                                              cAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGc                                           >  1:101861/1‑93 (MQ=255)

                                              cAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGc                                           >  1:1435283/1‑93 (MQ=255)

                                                gTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAg                                         <  1:1732872/93‑1 (MQ=255)

                                                gTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAg                                         <  1:3053933/93‑1 (MQ=255)

                                                    tCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGa                                     <  2:2740479/93‑1 (MQ=255)

                                                    tCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGa                                     <  1:1525716/93‑1 (MQ=255)

                                                    tCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGa                                     >  1:2536519/1‑93 (MQ=255)

                                                    tCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGa                                     >  1:2765253/1‑93 (MQ=255)

                                                    tCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGa                                     <  2:3323970/93‑1 (MQ=255)

                                                     cGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGAt                                    <  2:4009711/93‑1 (MQ=255)

                                                     cGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGAt                                    <  1:3879866/93‑1 (MQ=255)

                                                     cGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGAt                                    <  2:4200582/93‑1 (MQ=255)

                                                       aCTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTc                                  >  1:1908928/1‑93 (MQ=255)

                                                       aCTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTc                                  >  1:3402638/1‑93 (MQ=255)

                                                        cTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCa                                 >  2:1671426/1‑93 (MQ=255)

                                                         tttGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAg                                <  2:1433228/93‑1 (MQ=255)

                                                          ttGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAgg                               >  1:3580119/1‑93 (MQ=255)

                                                           tGGCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGc                              >  1:3604350/1‑93 (MQ=255)

                                                            ggttgtgGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGcc                             <  1:4149776/90‑1 (MQ=255)

                                                            ggCTGTGGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGcc                             <  2:105741/93‑1 (MQ=255)

                                                            cgttgtgGTTCGGGCATTGTGGCGATCGAGTCGTGGAATCTGGGTTCAGCAAAAGCAATTGGTATTGATATCGATCCGCAGGCGATTCAGGcc                             <  2:1082973/90‑1 (MQ=255)

                                                               tgtgGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGc                          >  2:4139177/1‑93 (MQ=255)

                                                                gtgGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGcc                         >  2:3153905/1‑93 (MQ=255)

                                                                 tgGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCcg                        <  1:2981219/93‑1 (MQ=255)

                                                                 tgGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCcg                        <  1:3703081/93‑1 (MQ=255)

                                                                 tgGTTCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGCATTGATATCGATCCGCAGGCGATTCAGGCCAGcgg                        >  1:468468/1‑91 (MQ=255)

                                                                    ttCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGa                     <  1:663252/93‑1 (MQ=255)

                                                                    ttCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGa                     <  1:164285/93‑1 (MQ=255)

                                                                     tCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                    >  1:3335465/1‑93 (MQ=255)

                                                                     tCCGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                    <  2:4256873/93‑1 (MQ=255)

                                                                       cGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGTCAGCCGCGATaa                  >  1:2961065/1‑93 (MQ=255)

                                                                       cGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATaa                  <  1:1541050/93‑1 (MQ=255)

                                                                       cGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATaa                  <  2:1984022/93‑1 (MQ=255)

                                                                       cGGCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATaa                  <  1:4222217/93‑1 (MQ=255)

                                                                         gCATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGACATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACg                >  1:3301201/1‑93 (MQ=255)

                                                                          cATTCTGGCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGc               <  1:454384/93‑1 (MQ=255)

                                                                                ggCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACg         <  2:821975/93‑1 (MQ=255)

                                                                                ggCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACg         >  1:3185393/1‑93 (MQ=255)

                                                                                 gCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt        >  1:3628074/1‑93 (MQ=255)

                                                                                 gCGATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt        <  2:1472845/93‑1 (MQ=255)

                                                                                   gATCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTaa      >  1:2135194/1‑93 (MQ=255)

                                                                                    aTCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTaag     >  2:2983288/1‑92 (MQ=255)

                                                                                    aTCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAAt     <  1:2994811/93‑1 (MQ=255)

                                                                                    aTCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAAt     >  2:4639226/1‑93 (MQ=255)

                                                                                     tCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATg    <  1:4302164/93‑1 (MQ=255)

                                                                                     tCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATg    <  1:4501044/93‑1 (MQ=255)

                                                                                     tCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATg    >  2:3867154/1‑93 (MQ=255)

                                                                                     tCGCGACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATg    >  1:901608/1‑93 (MQ=255)

                                                                                      cgcgACGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATgg   <  2:3389636/93‑1 (MQ=255)

                                                                                      cgcgACGCTGAAACTGGGTGCAGCAAAAGCAATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATgg   <  2:835262/93‑1 (MQ=255)

                                                                                       gcgacgCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGc  <  2:1783861/93‑1 (MQ=255)

TCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGC  >  NC_000913/3407142‑3407321

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