BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,097,623	T→C	M162V (ATG→GTG)	rhaT ←	RhaT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,097,623	0	T	C	99.7%	994.4 / ‑5.5	319	M162V (ATG→GTG)	rhaT	RhaT
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (197/120); total (198/121)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 331 total aligned reads displayed.

TCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCGCCT  >  NC_000913/4097533‑4097711

  gCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTg                                                                                      >  1:2594310/1‑93 (MQ=255)

    aaaGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTgcg                                                                                    >  2:3700133/1‑93 (MQ=255)

    aaaGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTgcg                                                                                    >  1:7861457/1‑93 (MQ=255)

     aaGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTcgcgt                                                                                   >  1:2295913/1‑92 (MQ=255)

     aaGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTgcgc                                                                                   >  1:4606708/1‑93 (MQ=255)

       ggACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGctc                                                                                 >  1:3941486/1‑93 (MQ=255)

       ggACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGctc                                                                                 <  1:1063764/93‑1 (MQ=255)

          cATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCttt                                                                              <  1:6969984/93‑1 (MQ=255)

           aTCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTc                                                                             <  1:977586/93‑1 (MQ=255)

           aTCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTc                                                                             <  2:2956678/93‑1 (MQ=255)

             cccGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCaa                                                                           <  1:2115246/93‑1 (MQ=255)

             cccGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCaa                                                                           <  1:2013007/93‑1 (MQ=255)

              ccGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAAc                                                                          <  1:5272424/93‑1 (MQ=255)

              ccGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAAc                                                                          <  1:1667313/93‑1 (MQ=255)

               cGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACt                                                                         >  2:7645094/1‑93 (MQ=255)

                 gCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCGCCTTGCGCTCTTTCAACTGc                                                                       <  1:5053874/93‑1 (MQ=255)

                  cAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGcc                                                                      <  2:5020182/93‑1 (MQ=255)

                   agagAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGccc                                                                     <  2:3746569/93‑1 (MQ=255)

                   agagAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGccc                                                                     <  1:262559/93‑1 (MQ=255)

                   agagAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGccc                                                                     >  1:6626738/1‑93 (MQ=255)

                    gagaAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCg                                                                    >  1:3524828/1‑93 (MQ=255)

                      gaAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGgc                                                                  >  2:5834218/1‑93 (MQ=255)

                       aaaaTGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGgcg                                                                 >  2:72417/1‑93 (MQ=255)

                         aaTGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGcgcg                                                               >  2:7825049/1‑93 (MQ=255)

                           tGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAg                                                             <  2:1080430/93‑1 (MQ=255)

                           tGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAg                                                             >  1:1171348/1‑93 (MQ=255)

                                  caTCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAAt                                                      >  1:6003525/1‑93 (MQ=255)

                                   aTCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTATTTCAACTGCCCGGCGCGAGTTACAATc                                                     >  1:4793764/1‑93 (MQ=255)

                                    tCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATcc                                                    >  2:6026316/1‑93 (MQ=255)

                                    tCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATcc                                                    >  1:2359419/1‑93 (MQ=255)

                                      aCCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCt                                                  <  2:3227461/93‑1 (MQ=255)

                                       ccGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTa                                                 <  1:4285873/93‑1 (MQ=255)

                                       ccGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTa                                                 <  2:6095111/93‑1 (MQ=255)

                                        cGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTAc                                                <  2:638165/93‑1 (MQ=255)

                                        cGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTAc                                                <  1:2684901/93‑1 (MQ=255)

                                        cGCCAGCACCAGCCCGTTTGTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGGGACAAACCCTAc                                                >  1:7027596/1‑93 (MQ=255)

                                        cGCCAACACCAGCCCTTTTTTCAGATTGACCTCTTCGGCTTTAATGCCCCCCTTGCGCTCTTTCAACTTCCCGGCGCGAGTTACAATCCCTAc                                                <  1:1801908/93‑1 (MQ=255)

                                         gCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACg                                               <  1:4107016/93‑1 (MQ=255)

                                           cAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGcc                                             <  2:2320020/93‑1 (MQ=255)

                                            aGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCa                                            >  2:551673/1‑93 (MQ=255)

                                             gCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCaa                                           <  1:1102531/93‑1 (MQ=255)

                                              cACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCact                                          >  1:7535396/1‑91 (MQ=255)

                                               aCCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATc                                         >  2:2199451/1‑93 (MQ=255)

                                               aCCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATc                                         >  1:4494536/1‑93 (MQ=255)

                                                ccAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCa                                        >  2:1323880/1‑93 (MQ=255)

                                                ccAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCa                                        >  2:5190106/1‑93 (MQ=255)

                                                 cAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAg                                       >  1:6785866/1‑93 (MQ=255)

                                                  aGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAgc                                      >  2:4097604/1‑93 (MQ=255)

                                                   gCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAgcg                                     >  2:4572014/1‑93 (MQ=255)

                                                   gCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCAACTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAgcg                                     >  2:7692473/1‑93 (MQ=255)

                                                     ccTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGcc                                   >  1:3628265/1‑93 (MQ=255)

                                                     ccTTTTTTCAAATTGATCTCTTCGGCTTTAATCCCCCCCGTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGcc                                   <  1:1092987/93‑1 (MQ=255)

                                                      cTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGcca                                  >  2:4279026/1‑93 (MQ=255)

                                                       ttttttCAGATTGAACTCTTCGGCTTTGATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccac                                 >  2:3996293/1‑93 (MQ=255)

                                                       ttttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccac                                 >  2:2548539/1‑93 (MQ=255)

                                                       ttttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccac                                 >  2:7903193/1‑93 (MQ=255)

                                                       ttttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccac                                 >  1:4350817/1‑93 (MQ=255)

                                                       ttttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccac                                 >  1:7361876/1‑93 (MQ=255)

                                                          tttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAg                              <  2:1682524/93‑1 (MQ=255)

                                                           ttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGa                             <  1:5828045/93‑1 (MQ=255)

                                                            tCAGATTGACCTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGaa                            <  1:6917034/93‑1 (MQ=255)

                                                            tCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGaa                            <  1:3409553/93‑1 (MQ=255)

                                                               gATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGc                         >  2:5923384/1‑93 (MQ=255)

                                                               gATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGc                         >  2:6487783/1‑93 (MQ=255)

                                                               gATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGc                         >  2:4644664/1‑93 (MQ=255)

                                                               gATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGc                         >  1:6504315/1‑93 (MQ=255)

                                                                atggAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGcc                        <  1:622855/90‑1 (MQ=255)

                                                                 ttGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCg                       <  2:713779/93‑1 (MQ=255)

                                                                 ttGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCg                       >  2:7479407/1‑93 (MQ=255)

                                                                   gAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAg                     >  1:4911520/1‑93 (MQ=255)

                                                                     aCTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCa                   >  2:700987/1‑93 (MQ=255)

                                                                      ctctTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCaa                  <  1:6488118/93‑1 (MQ=255)

                                                                      ctctTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCaa                  <  1:3141768/93‑1 (MQ=255)

                                                                      ctctTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCaa                  <  1:7279629/93‑1 (MQ=255)

                                                                       tctTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGATACAATCCCTACGCCAATCAACGCCACCAGAACGCCGAGCAAc                 >  2:5581069/1‑93 (MQ=255)

                                                                        ctTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACg                >  2:2744114/1‑93 (MQ=255)

                                                                         ttCGGCTTTAATGCCCACCTTGCGCTCTTTCACCTGCCCGGCGCGAGTTACAATCCCTACGCCCATAAGCGCCACCAGAAAGCCGAGCAACGt               >  2:5148815/1‑93 (MQ=255)

                                                                         ttCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGt               >  1:880282/1‑93 (MQ=255)

                                                                          tCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTc              <  2:4868111/93‑1 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa             <  2:676272/93‑1 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa             >  1:538887/1‑93 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa             <  2:2024399/93‑1 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa             <  2:3515957/93‑1 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa             <  2:1502407/93‑1 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa             <  1:406624/93‑1 (MQ=255)

                                                                             gCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATg           >  1:1337904/1‑93 (MQ=255)

                                                                             gCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATg           <  1:7174927/93‑1 (MQ=255)

                                                                             gCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATg           >  2:1726914/1‑93 (MQ=255)

                                                                              cTTTGATGCCCACCTTGCGCTCTTTCAAGTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGACAGGc          >  2:7420023/1‑93 (MQ=255)

                                                                              cTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGc          >  1:515431/1‑93 (MQ=255)

                                                                              cTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGc          >  2:3353720/1‑93 (MQ=255)

                                                                               tttAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCg         >  1:2376359/1‑93 (MQ=255)

                                                                                 tAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTc       >  1:3765777/1‑93 (MQ=255)

                                                                                 tAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTc       >  1:2909268/1‑93 (MQ=255)

                                                                                  aaTGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTcc      <  1:1537757/93‑1 (MQ=255)

                                                                                   aTGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCg     >  2:4428762/1‑93 (MQ=255)

                                                                                   aTGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCg     >  2:2813807/1‑93 (MQ=255)

                                                                                    tGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCGc    >  1:360289/1‑93 (MQ=255)

                                                                                     gCCCACCTTGCGCTCTTACAACTGCCCGGCGCGAGTTACAATCCCTCCGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGCACGcc   <  2:4905929/93‑1 (MQ=255)

                                                                                      cccACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCGCCt  <  1:7729758/93‑1 (MQ=255)

TCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCGCCT  >  NC_000913/4097533‑4097711

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