BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,457,112	A→G	V23A (GTA→GCA)	b2341 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,457,112	0	A	G	83.4%	526.0 / 102.8	241	V23A (GTA→GCA)	b2341	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (15/25); new base G (79/122); total (94/147)
Fisher's exact test for biased strand distribution p-value = 8.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 242 total aligned reads displayed.

CTCGCAACTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCC  >  NC_000913/2457020‑2457203

   gCAACTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGt                                                                                          <  2:831292/93‑1 (MQ=255)

    cAACTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTc                                                                                         <  2:351519/93‑1 (MQ=255)

      aCTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGa                                                                                       <  2:4676127/93‑1 (MQ=255)

      aCTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACACCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGa                                                                                       <  1:6223443/93‑1 (MQ=255)

       ctctTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGAt                                                                                      >  2:5746326/1‑93 (MQ=255)

       ctctTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGAt                                                                                      <  1:2627219/93‑1 (MQ=255)

        tctTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATg                                                                                     <  2:218513/93‑1 (MQ=255)

        tctTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATg                                                                                     >  1:2255897/1‑93 (MQ=255)

          tttGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGt                                                                                   >  1:6572945/1‑93 (MQ=255)

          tttGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGt                                                                                   >  2:1049796/1‑93 (MQ=255)

            tGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGa                                                                                 <  1:796501/93‑1 (MQ=255)

             gTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGAt                                                                                <  1:6475523/93‑1 (MQ=255)

               tttCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATaa                                                                              <  1:6635436/93‑1 (MQ=255)

                ttCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAAc                                                                             <  1:1480604/93‑1 (MQ=255)

                 tCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACg                                                                            >  1:6059593/1‑93 (MQ=255)

                  cACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACgg                                                                           <  2:3923747/93‑1 (MQ=255)

                   aCGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGc                                                                          >  2:7375212/1‑93 (MQ=255)

                     ggAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCaa                                                                        >  2:3474824/1‑93 (MQ=255)

                     ggAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCaa                                                                        <  1:5740283/93‑1 (MQ=255)

                      gAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAAt                                                                       <  1:7082590/93‑1 (MQ=255)

                      gAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAAt                                                                       <  2:3267364/93‑1 (MQ=255)

                        gTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAAtgt                                                                     >  2:5795013/1‑93 (MQ=255)

                        gTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAAtgt                                                                     <  2:5723230/93‑1 (MQ=255)

                        gTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAAtgt                                                                     <  2:3670222/93‑1 (MQ=255)

                         ttGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAAtgtt                                                                    >  1:6717644/1‑93 (MQ=255)

                          tGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAAtgttg                                                                   <  2:2525011/93‑1 (MQ=255)

                           gCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAAtgttgt                                                                  >  2:6293985/1‑93 (MQ=255)

                           gCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAAtgttgt                                                                  >  2:3329376/1‑93 (MQ=255)

                           gCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTATCACCCGGTGCGTCGATGGTGATAACGGCAAtgttgt                                                                  <  1:40567/93‑1 (MQ=255)

                            cTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTc                                                                 >  1:2985275/1‑93 (MQ=255)

                            cTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTc                                                                 <  1:5861221/93‑1 (MQ=255)

                            cTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTc                                                                 >  2:1311622/1‑93 (MQ=255)

                             ttAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTcc                                                                >  2:1372558/1‑93 (MQ=255)

                              taataaTGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCa                                                               <  1:2652890/93‑1 (MQ=255)

                               aataatGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAg                                                              >  2:7247514/1‑93 (MQ=255)

                                ataatGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGa                                                             <  2:1557708/93‑1 (MQ=255)

                                ataatGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGa                                                             <  2:4025244/93‑1 (MQ=255)

                                 taatGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGAc                                                            <  2:6355959/93‑1 (MQ=255)

                                 taatGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGAc                                                            <  2:4507698/93‑1 (MQ=255)

                                  aatGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGGTGGTGATAACGGCAATGTTGTCCAGACg                                                           >  2:3446080/1‑93 (MQ=255)

                                    tGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGaa                                                         <  2:6260569/93‑1 (MQ=255)

                                     ggCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAAc                                                        <  1:4487208/93‑1 (MQ=255)

                                      gcgcgcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACa                                                       >  1:1044268/1‑93 (MQ=255)

                                       cgcgcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACAt                                                      >  1:4176509/1‑93 (MQ=255)

                                          gcACCTGCGAGGCAATCTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTaa                                                   <  1:3460059/93‑1 (MQ=255)

                                          gcACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTaa                                                   <  1:171364/93‑1 (MQ=255)

                                           cACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGCCCAGACGAACATTAAg                                                  <  2:2708317/93‑1 (MQ=255)

                                            aCCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAgg                                                 <  2:335724/93‑1 (MQ=255)

                                             ccTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAggg                                                >  1:2278324/1‑93 (MQ=255)

                                               tGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTa                                              <  1:2485427/93‑1 (MQ=255)

                                                 cGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTaaa                                            <  2:2983395/93‑1 (MQ=255)

                                                  gAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAAc                                           <  1:3085515/93‑1 (MQ=255)

                                                   aGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACg                                          >  2:6818345/1‑93 (MQ=255)

                                                   aGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACg                                          >  1:4606986/1‑93 (MQ=255)

                                                   aGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACg                                          >  2:6707404/1‑93 (MQ=255)

                                                   aGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACg                                          >  1:5391232/1‑93 (MQ=255)

                                                    ggCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGc                                         <  1:1910555/93‑1 (MQ=255)

                                                       aaaCTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGa                                      <  1:2580616/93‑1 (MQ=255)

                                                       aaaCTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGa                                      <  1:7234369/93‑1 (MQ=255)

                                                        aaCTCCGCCTTCAGGGTATTCATTTTCTCCCCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGAt                                     <  1:5623433/93‑1 (MQ=255)

                                                        aaCTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGAt                                     >  1:4962068/1‑93 (MQ=255)

                                                          cTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGt                                   >  2:5945859/1‑93 (MQ=255)

                                                           tCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTc                                  >  2:161655/1‑93 (MQ=255)

                                                           tCCGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTc                                  >  1:2849404/1‑93 (MQ=255)

                                                            ccGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCa                                 <  2:3131603/93‑1 (MQ=255)

                                                            ccGCCTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCa                                 <  1:5092580/93‑1 (MQ=255)

                                                            ccGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCa                                 <  2:4243413/93‑1 (MQ=255)

                                                               ccTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCAttt                              >  1:5438962/1‑93 (MQ=255)

                                                                cTTCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTc                             <  1:1950491/93‑1 (MQ=255)

                                                                 ttCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTcc                            <  1:3228916/93‑1 (MQ=255)

                                                                 ttCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTcc                            >  1:83614/1‑93 (MQ=255)

                                                                  tCAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCa                           <  2:3622712/93‑1 (MQ=255)

                                                                   cAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCat                          <  1:3797788/93‑1 (MQ=255)

                                                                   cAGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCat                          >  1:738569/1‑93 (MQ=255)

                                                                   cAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCat                          <  2:2900268/93‑1 (MQ=255)

                                                                    aGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCatt                         >  1:2749563/1‑93 (MQ=255)

                                                                    aGGGTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCatt                         >  1:2144908/1‑93 (MQ=255)

                                                                     gggTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCatta                        <  1:481148/93‑1 (MQ=255)

                                                                     gggTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCatta                        >  1:3135677/1‑93 (MQ=255)

                                                                     gggTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCatta                        >  1:4552300/1‑93 (MQ=255)

                                                                      ggTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCattat                       >  2:5017232/1‑93 (MQ=255)

                                                                       gTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCattatt                      >  2:7007771/1‑93 (MQ=255)

                                                                       gTATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCattatt                      >  1:431716/1‑93 (MQ=255)

                                                                        tATTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTATACGCTGATGTCATTTCCATTATTc                     <  2:7363224/93‑1 (MQ=255)

                                                                         aTTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTcc                    <  2:6890381/93‑1 (MQ=255)

                                                                         aTTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTcc                    >  1:128986/1‑93 (MQ=255)

                                                                         aTTCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTcc                    <  1:81085/93‑1 (MQ=255)

                                                                          ttCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCg                   <  2:5982470/93‑1 (MQ=255)

                                                                           tCATTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGc                  <  1:3376634/93‑1 (MQ=255)

                                                                             aTTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCt                >  1:5106864/1‑93 (MQ=255)

                                                                             aTTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCt                >  2:90514/1‑93 (MQ=255)

                                                                             aTTTTCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACACTGATGTCATTTCCATTATTCCGCCt                >  2:2271112/1‑93 (MQ=255)

                                                                               tttCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTcc              <  1:1813167/93‑1 (MQ=255)

                                                                                ttCTCACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCa             >  1:3343665/1‑93 (MQ=255)

                                                                                 tctcACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAg            >  2:5576904/1‑93 (MQ=255)

                                                                                 tctcACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAg            >  2:3735714/1‑93 (MQ=255)

                                                                                   tcACCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGaa          >  1:4924876/1‑93 (MQ=255)

                                                                                     aCCCGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAAcc        <  1:3989546/93‑1 (MQ=255)

                                                                                       ccGGTGCGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCAt      <  2:1798244/93‑1 (MQ=255)

                                                                                           tACGTCGATGGTGAGAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGcc  >  2:6437372/1‑93 (MQ=255)

CTCGCAACTCTTTGTTTTCACGGAGTTGCTTAATAATGGCGCGCACCTGCGAGGCAAACTCCGCCTTCAGGGTATTCATTTTCTCACCCGGTACGTCGATGGTGATAACGGCAATGTTGTCCAGACGAACATTAAGGGTAAACGCTGATGTCATTTCCATTATTCCGCCTCCAGAACCATTGCC  >  NC_000913/2457020‑2457203

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