BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,655,245	T→C	L91L (TTA→TTG)	ynfC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,655,245	0	T	C	98.9%	596.8 / ‑3.4	184	L91L (TTA→TTG)	ynfC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (98/84); total (99/85)
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.39e-01

Only 100 of 188 total aligned reads displayed.

TAAACGTACTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGTAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAATGTCTGAGT  >  NC_000913/1655155‑1655334

tAAACGTACTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCaa                                                                                         <  2:3008986/93‑1 (MQ=255)

  aaCGTACTCTTTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAAtt                                                                                       >  2:3055376/1‑93 (MQ=255)

   aCGTACTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTc                                                                                      <  2:2345739/93‑1 (MQ=255)

   aCGTACTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTc                                                                                      <  2:772475/93‑1 (MQ=255)

    cGTACTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCg                                                                                     <  1:4561567/93‑1 (MQ=255)

      tACTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAg                                                                                   >  1:113813/1‑93 (MQ=255)

       aCTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGt                                                                                  >  1:4435954/1‑93 (MQ=255)

       aCTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGt                                                                                  <  2:2143060/93‑1 (MQ=255)

       aCTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTAGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGt                                                                                  >  2:4296091/1‑93 (MQ=255)

          ctctACTCCAGCATCAAAGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGAAACAGCAATTCGAATGaa                                                                               <  2:4219757/93‑1 (MQ=255)

            cttctCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATc                                                                             >  1:911069/1‑93 (MQ=255)

                tcCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCaaaa                                                                         >  1:2545261/1‑93 (MQ=255)

                tcCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCaaaa                                                                         >  1:2719955/1‑93 (MQ=255)

                 ccAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAAc                                                                        >  1:3895009/1‑93 (MQ=255)

                   aGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAg                                                                      <  2:2755316/93‑1 (MQ=255)

                    gCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGc                                                                     <  2:3359243/93‑1 (MQ=255)

                    gCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGc                                                                     >  1:2856025/1‑93 (MQ=255)

                     cGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGcc                                                                    <  2:1332955/93‑1 (MQ=255)

                     cGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGcc                                                                    <  2:3988360/93‑1 (MQ=255)

                      gTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCct                                                                   <  1:2017672/93‑1 (MQ=255)

                      gTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCct                                                                   >  2:641733/1‑93 (MQ=255)

                       tctcGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCctt                                                                  >  2:2593180/1‑93 (MQ=255)

                        ctcGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCcttc                                                                 <  1:3678385/93‑1 (MQ=255)

                        ctcGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCcttc                                                                 <  1:3293793/93‑1 (MQ=255)

                         tcGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCcttct                                                                <  2:2523511/93‑1 (MQ=255)

                         tcGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCcttct                                                                >  2:1131989/1‑93 (MQ=255)

                          cGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCcttctt                                                               <  1:2112486/93‑1 (MQ=255)

                          cGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCcttctt                                                               >  2:2525119/1‑93 (MQ=255)

                           ggCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCttcttc                                                              >  1:4335333/1‑93 (MQ=255)

                           ggCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCttcttc                                                              <  1:509374/93‑1 (MQ=255)

                            gCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTcc                                                             >  1:1749640/1‑93 (MQ=255)

                              aTCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGa                                                           >  1:863245/1‑93 (MQ=255)

                               tCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGAc                                                          >  2:3807568/1‑93 (MQ=255)

                                 aCGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACaa                                                        <  2:3859834/93‑1 (MQ=255)

                                   ggTAATAATTGGCGTCCATTACCAGAGCGACCACGGTATTATTTTCCAGATCCATCAATTCAACTCAATCAAAACAGCCTCCTGCCGACAAAg                                                      >  1:1934746/1‑93 (MQ=255)

                                    gTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGt                                                     <  1:1367393/93‑1 (MQ=255)

                                    gTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGt                                                     <  2:146873/93‑1 (MQ=255)

                                     taataaTTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTc                                                    <  1:208401/93‑1 (MQ=255)

                                        taatTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCa                                                 <  2:3608018/93‑1 (MQ=255)

                                          atTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGa                                               >  2:4117993/1‑93 (MQ=255)

                                          atTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGa                                               >  2:2682338/1‑93 (MQ=255)

                                           ttGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGaa                                              >  2:2069022/1‑93 (MQ=255)

                                           ttGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGaa                                              >  1:3941011/1‑93 (MQ=255)

                                            tGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGaaa                                             <  2:452398/93‑1 (MQ=255)

                                               cGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAacac                                          <  1:4443271/93‑1 (MQ=255)

                                               cGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAacac                                          <  2:1383656/93‑1 (MQ=255)

                                                 tCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGt                                        <  1:1716855/93‑1 (MQ=255)

                                                 tCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGt                                        <  1:3805685/93‑1 (MQ=255)

                                                 gccAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGt                                        <  1:3805684/92‑1 (MQ=255)

                                                  ccAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGtt                                       <  2:1174087/93‑1 (MQ=255)

                                                   cAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGttt                                      >  2:1539505/1‑93 (MQ=255)

                                                   cAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGttt                                      <  1:804624/93‑1 (MQ=255)

                                                   cAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGttt                                      <  1:1801567/93‑1 (MQ=255)

                                                   cAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGttt                                      >  2:693619/1‑93 (MQ=255)

                                                   cAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGttt                                      >  1:3467820/1‑93 (MQ=255)

                                                    aGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTc                                     >  1:4242839/1‑93 (MQ=255)

                                                    aGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTc                                     <  1:4141747/93‑1 (MQ=255)

                                                    aGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTc                                     <  2:3573394/93‑1 (MQ=255)

                                                    aGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTc                                     >  2:3228650/1‑93 (MQ=255)

                                                     gTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCg                                    <  1:3100332/93‑1 (MQ=255)

                                                     gTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCg                                    >  1:2482119/1‑93 (MQ=255)

                                                      tACCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGt                                   <  1:2138356/93‑1 (MQ=255)

                                                       aCCAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTc                                  >  2:1058201/1‑93 (MQ=255)

                                                         cAGAGCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCAc                                <  1:3735576/93‑1 (MQ=255)

                                                          agagCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACt                               >  2:2852341/1‑93 (MQ=255)

                                                          agagCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACt                               <  2:4629886/93‑1 (MQ=255)

                                                          agagCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACt                               >  1:1984304/1‑93 (MQ=255)

                                                           gagCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACtt                              >  1:1229266/1‑93 (MQ=255)

                                                            agCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTc                             <  2:4079038/93‑1 (MQ=255)

                                                            agCGACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTc                             <  1:3557033/93‑1 (MQ=255)

                                                               gACCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCAcc                          <  1:2986285/93‑1 (MQ=255)

                                                                aCCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCt                         <  2:4016098/93‑1 (MQ=255)

                                                                aCCACGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCt                         <  2:1505523/93‑1 (MQ=255)

                                                                 ccACGGTATTATTTTCCAGATCCAGCAATTCGAGCGAATCAAAAAAGCCTTCTTCCGACAAAGTCCCAGAAACACGATTCGTCACTTCACCtt                        >  2:1861085/1‑93 (MQ=255)

                                                                  cACGGTATTATTTGCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTg                       >  1:3136552/1‑93 (MQ=255)

                                                                   ccGGTATTATTTTCCAGATCCAGTAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGc                      <  1:4029258/92‑1 (MQ=255)

                                                                   aCGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGc                      >  2:4649373/1‑93 (MQ=255)

                                                                   aCGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGc                      >  1:1009113/1‑93 (MQ=255)

                                                                   aCGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGc                      >  2:1849678/1‑93 (MQ=255)

                                                                    cGGTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCt                     <  1:2284948/93‑1 (MQ=255)

                                                                     ggTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTc                    <  2:2169876/93‑1 (MQ=255)

                                                                      gTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCa                   >  2:4620480/1‑93 (MQ=255)

                                                                      gTATTATTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCa                   >  2:967104/1‑93 (MQ=255)

                                                                       tattatTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCAt                  >  1:2488201/1‑93 (MQ=255)

                                                                       tattatTTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCAt                  >  2:1989419/1‑93 (MQ=255)

                                                                        attattTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATc                 >  1:3121570/1‑93 (MQ=255)

                                                                         ttattTTCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCCTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATcc                <  2:2404219/93‑1 (MQ=255)

                                                                             tttCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAACCACGTTTCGTCACTTCACCTTGCTCATCCATTa            <  2:4249595/93‑1 (MQ=255)

                                                                              ttCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTaa           >  1:2771157/1‑93 (MQ=255)

                                                                              ttCCAGAACCAGCAATTCGAGTGAATCAAATCAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACGTCACCTTGCTCATCCATTaa           <  2:4746865/93‑1 (MQ=255)

                                                                               tCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAAt          <  1:1376455/93‑1 (MQ=255)

                                                                               tCCAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAAt          >  1:3171811/1‑93 (MQ=255)

                                                                               tCCAGATCCAGCAATTCGACTGAATCAAATCAGCCTTCTACCGACAAAGTCCCAGAAACTCGTTTCGTCACTTCACCTTGCTCATCCATTAAt          <  2:1888560/93‑1 (MQ=255)

                                                                                 cAGATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAATGt        >  2:1152357/1‑93 (MQ=255)

                                                                                   gATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAATGTCt      >  2:55487/1‑93 (MQ=255)

                                                                                   gATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAATGTCt      <  1:2651053/93‑1 (MQ=255)

                                                                                   gATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAATGTCt      <  2:4090671/93‑1 (MQ=255)

                                                                                   gATCCAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAATGTCt      >  2:3380322/1‑93 (MQ=255)

                                                                                      ccAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAATGTCTGAg   <  1:1899437/93‑1 (MQ=255)

                                                                                       cAGCAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCATTTCACCTTGCTCATCCATTAATGTCTGAgg  >  2:724202/1‑92 (MQ=255)

TAAACGTACTCTCTTCTCCAGCGTCTCGGCATCACGGTAATAATTGGCGTCCAGTACCAGAGCGACCACGGTATTATTTTCCAGATCCAGTAATTCGAGTGAATCAAAACAGCCTTCTTCCGACAAAGTCCCAGAAACACGTTTCGTCACTTCACCTTGCTCATCCATTAATGTCTGAGT  >  NC_000913/1655155‑1655334

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