BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,476,619	A→G	V91A (GTG→GCG)	b4256 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,476,619	0	A	G	97.8%	549.8 / 1.2	180	V91A (GTG→GCG)	b4256	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/2); new base G (80/96); total (82/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

Only 100 of 183 total aligned reads displayed.

ATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAA  >  NC_000913/4476529‑4476708

aTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCt                                                                                         >  1:632146/1‑93 (MQ=255)

 tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCtt                                                                                        >  2:963503/1‑93 (MQ=255)

  ccGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTc                                                                                       <  2:1110327/93‑1 (MQ=255)

   cGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTca                                                                                      <  2:784981/93‑1 (MQ=255)

    ggCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTcac                                                                                     >  2:3638634/1‑93 (MQ=255)

    ggCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTcac                                                                                     <  2:4439062/93‑1 (MQ=255)

      cATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACaa                                                                                   <  2:3975629/93‑1 (MQ=255)

       aTAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAAt                                                                                  <  1:4052315/93‑1 (MQ=255)

       aTAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAAt                                                                                  >  1:3467568/1‑93 (MQ=255)

          aaCAAAACTTACGTTGTTAATGTTGGAAGGCCCGCGACTTGAACACGCGGCCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGa                                                                               <  2:3605385/93‑1 (MQ=255)

          aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGa                                                                               <  2:259394/93‑1 (MQ=255)

          aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGa                                                                               <  2:3869761/93‑1 (MQ=255)

           aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGAc                                                                              >  2:815296/1‑93 (MQ=255)

            cAAAACGTACGTTGTTAATGTTCGAATGCTCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCgcct                                                                             >  2:2497732/1‑90 (MQ=255)

            cAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACa                                                                             >  1:852911/1‑93 (MQ=255)

             aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACAt                                                                            >  2:2909305/1‑93 (MQ=255)

              aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATg                                                                           <  1:2171643/93‑1 (MQ=255)

               aaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATgg                                                                          >  2:2318760/1‑93 (MQ=255)

               aaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATgc                                                                          >  1:2284090/1‑92 (MQ=255)

                acgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGc                                                                         <  2:2819324/93‑1 (MQ=255)

                 cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGcc                                                                        >  1:1675286/1‑93 (MQ=255)

                 cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGcc                                                                        <  2:149655/93‑1 (MQ=255)

                  gtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGccc                                                                       <  2:400836/93‑1 (MQ=255)

                   tacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGTATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCg                                                                      >  2:4195625/1‑93 (MQ=255)

                   tacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCg                                                                      <  1:4628722/93‑1 (MQ=255)

                    tcgtTGTTAATGTTGGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGt                                                                     <  1:1729584/92‑1 (MQ=255)

                     cgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCTACATCCGCGCTTCACAATCGACATGGCCCGTg                                                                    >  2:2189217/1‑93 (MQ=255)

                     cgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTg                                                                    >  2:2360705/1‑93 (MQ=255)

                      gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGc                                                                   <  2:1902391/93‑1 (MQ=255)

                      gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGc                                                                   >  2:1943813/1‑93 (MQ=255)

                      gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGc                                                                   >  1:2493064/1‑93 (MQ=255)

                      gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGc                                                                   <  1:4358594/93‑1 (MQ=255)

                       ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCa                                                                  >  1:2979611/1‑93 (MQ=255)

                        tgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAg                                                                 >  2:1318881/1‑93 (MQ=255)

                         gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGc                                                                <  1:2810415/93‑1 (MQ=255)

                         gtgaaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGc                                                                <  1:4701352/90‑1 (MQ=255)

                          ttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGcc                                                               >  2:3790210/1‑93 (MQ=255)

                            aaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCaa                                                             <  1:3161088/93‑1 (MQ=255)

                            aaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCaa                                                             <  1:3371934/93‑1 (MQ=255)

                            aaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCaa                                                             <  2:2658664/93‑1 (MQ=255)

                             aTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAg                                                            >  1:4187128/1‑93 (MQ=255)

                              tGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCTTGCAGCCAAgc                                                           >  1:3331151/1‑93 (MQ=255)

                              tGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAgc                                                           <  1:2239031/93‑1 (MQ=255)

                                ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAgcgc                                                         <  1:244317/93‑1 (MQ=255)

                                ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAgcgc                                                         >  2:2987419/1‑93 (MQ=255)

                                  cGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                                                       <  1:3712925/93‑1 (MQ=255)

                                    aaTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAg                                                     <  2:503206/93‑1 (MQ=255)

                                    aaTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAg                                                     <  2:2071594/93‑1 (MQ=255)

                                     aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCTCATAGt                                                    >  2:2519416/1‑93 (MQ=255)

                                     aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                                                    <  2:4606373/93‑1 (MQ=255)

                                     aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                                                    <  1:3163971/93‑1 (MQ=255)

                                      tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                                                   <  1:4002617/93‑1 (MQ=255)

                                      tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCACCCAAGCGCCTAGTc                                                   >  1:2615951/1‑93 (MQ=255)

                                        cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                 <  2:2815721/93‑1 (MQ=255)

                                           gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                              <  2:1856961/93‑1 (MQ=255)

                                           gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                              >  1:1067325/1‑93 (MQ=255)

                                           gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                              >  2:472057/1‑93 (MQ=255)

                                           gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                              >  1:4430411/1‑93 (MQ=255)

                                             gAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                            <  1:2612656/93‑1 (MQ=255)

                                              aGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCt                                           >  1:4396441/1‑93 (MQ=255)

                                               gTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                                          <  1:2069145/93‑1 (MQ=255)

                                                  gAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaaa                                       >  1:3519151/1‑93 (MQ=255)

                                                  gAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaaa                                       >  1:3624057/1‑93 (MQ=255)

                                                   aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                      >  1:4314636/1‑93 (MQ=255)

                                                      actcgGGCATTTTTATAGAGTTCGTGCAACATCCTCGTTTTTTAATCGACATTGTCCGTGTAGTCTAGCGCATAGTCGATTTGTTCTAAGccc                                   >  1:4203801/4‑93 (MQ=255)

                                                       cgcgGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCa                                  >  2:2139978/1‑93 (MQ=255)

                                                          gggCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAAt                               >  2:1616939/1‑93 (MQ=255)

                                                          gggCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAAt                               <  1:4550722/93‑1 (MQ=255)

                                                           ggCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATg                              >  2:4707250/1‑93 (MQ=255)

                                                             cATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                            <  2:1288689/93‑1 (MQ=255)

                                                             cATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                            <  1:193853/93‑1 (MQ=255)

                                                             cATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                            >  1:4445365/1‑93 (MQ=255)

                                                              aTTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTc                           <  1:898736/93‑1 (MQ=255)

                                                                   tttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                      <  1:2819573/93‑1 (MQ=255)

                                                                   tttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                      <  1:1139889/93‑1 (MQ=255)

                                                                   tttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                      <  2:711187/93‑1 (MQ=255)

                                                                   tttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                      <  1:1795142/93‑1 (MQ=255)

                                                                   tttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                      <  2:1978430/93‑1 (MQ=255)

                                                                     tAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAgc                    <  1:248492/93‑1 (MQ=255)

                                                                      agagTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAgcg                   >  2:4380659/1‑93 (MQ=255)

                                                                        agTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCa                 >  2:2752439/1‑93 (MQ=255)

                                                                         gtgcgcgCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCaa                <  2:1519411/90‑1 (MQ=255)

                                                                         gTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGAGCaa                >  2:3356668/1‑93 (MQ=255)

                                                                          ttCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAAt               <  1:4271902/93‑1 (MQ=255)

                                                                          ttCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAAt               >  2:2884776/1‑93 (MQ=255)

                                                                           tCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATg              <  1:1212209/93‑1 (MQ=255)

                                                                            cgcgCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATgg             <  1:3058553/93‑1 (MQ=255)

                                                                               gcAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCtt          <  2:4041380/93‑1 (MQ=255)

                                                                               gcAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCtt          >  2:174041/1‑93 (MQ=255)

                                                                                cAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTc         <  2:578803/93‑1 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt       >  2:3159650/1‑93 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt       >  1:1271817/1‑93 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt       >  2:4586403/1‑93 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt       >  1:2399369/1‑93 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt       <  1:1468774/93‑1 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt       >  2:624502/1‑93 (MQ=255)

                                                                                    aTCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCttt     <  2:3313081/93‑1 (MQ=255)

                                                                                    aTCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCttt     <  1:289276/93‑1 (MQ=255)

                                                                                    aTCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCttt     >  1:324560/1‑93 (MQ=255)

                                                                                     tCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTa    <  2:4586374/93‑1 (MQ=255)

ATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAA  >  NC_000913/4476529‑4476708

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