BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,019,057	A→G	V415V (GTA→GTG)	yigR →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,019,057	0	A	G	99.7%	1031.9 / ‑5.3	333	V415V (GTA→GTG)	yigR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (161/171); total (161/172)
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.71e-01

Only 100 of 351 total aligned reads displayed.

AAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGC  >  NC_000913/4018967‑4019146

aaaTCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTggg                                                                                         >  1:5242637/1‑93 (MQ=255)

aaaTCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTggg                                                                                         >  1:46598/1‑93 (MQ=255)

 aaTCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGa                                                                                        <  2:2668500/93‑1 (MQ=255)

  aTCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGAc                                                                                       >  2:1408525/1‑93 (MQ=255)

   tCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACg                                                                                      >  1:6353892/1‑93 (MQ=255)

     tGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGcc                                                                                    >  1:4410500/1‑93 (MQ=255)

        ttAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGc                                                                                 <  1:3653226/93‑1 (MQ=255)

        ttAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGc                                                                                 >  2:1253996/1‑93 (MQ=255)

          aaTGCGGCGCGTCGCTTCAATATGCAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCtt                                                                               <  2:4475714/93‑1 (MQ=255)

          aaTACGGCGCGTCGCTTCAATATGGAAGTCCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCtt                                                                               <  1:5547732/93‑1 (MQ=255)

           aTACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCttt                                                                              >  1:4234480/1‑93 (MQ=255)

            tgcGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTa                                                                             <  1:530624/91‑1 (MQ=255)

             aCGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTAt                                                                            >  1:2634453/1‑93 (MQ=255)

             aCGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTAt                                                                            <  1:2671609/93‑1 (MQ=255)

               ggCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATcc                                                                          <  1:3766721/93‑1 (MQ=255)

                  gcgTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCa                                                                       >  1:3363335/1‑93 (MQ=255)

                   cgTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCac                                                                      >  1:5862070/1‑92 (MQ=255)

                   cgTCACTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCTCaa                                                                      >  2:2022792/1‑93 (MQ=255)

                    gTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCGGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAAc                                                                     >  2:7451435/1‑93 (MQ=255)

                       gCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCg                                                                  >  2:5228544/1‑93 (MQ=255)

                       gCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCg                                                                  <  2:2644658/93‑1 (MQ=255)

                        cTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGa                                                                 <  2:5608252/93‑1 (MQ=255)

                         ttCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAt                                                                <  2:4833281/93‑1 (MQ=255)

                         ttCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAt                                                                >  2:5350723/1‑93 (MQ=255)

                         ttCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAt                                                                <  2:4131970/93‑1 (MQ=255)

                          tCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAtt                                                               >  2:604895/1‑93 (MQ=255)

                          tCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAtt                                                               <  1:1499250/93‑1 (MQ=255)

                           cAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGAttt                                                              <  1:1344544/93‑1 (MQ=255)

                            aaTATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTa                                                             <  2:5721626/93‑1 (MQ=255)

                            aaTATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTa                                                             >  1:4576738/1‑93 (MQ=255)

                            aaTATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTa                                                             >  1:7455423/1‑93 (MQ=255)

                              tatGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATg                                                           >  2:4975162/1‑93 (MQ=255)

                                tGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGa                                                         >  2:4686634/1‑93 (MQ=255)

                                tGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGa                                                         <  2:6271713/93‑1 (MQ=255)

                                 ggAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGaa                                                        <  2:148353/93‑1 (MQ=255)

                                 ggAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGaa                                                        >  2:7593930/1‑93 (MQ=255)

                                  gAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGaaa                                                       <  2:2548098/93‑1 (MQ=255)

                                   aaGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGaaaa                                                      <  2:4643078/93‑1 (MQ=255)

                                    aGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAAc                                                     <  2:5994243/93‑1 (MQ=255)

                                    aGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAAc                                                     <  2:6375406/93‑1 (MQ=255)

                                      tGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACgg                                                   >  2:288112/1‑93 (MQ=255)

                                       gCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGc                                                  <  2:6838623/93‑1 (MQ=255)

                                       gCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGc                                                  >  1:7176122/1‑93 (MQ=255)

                                          gCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGaa                                               >  1:1200069/1‑93 (MQ=255)

                                             gCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGcc                                            >  1:3492535/1‑93 (MQ=255)

                                             gCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGcc                                            >  1:2263585/1‑93 (MQ=255)

                                              cAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCt                                           >  1:4917052/1‑93 (MQ=255)

                                              cAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCt                                           <  2:1792129/93‑1 (MQ=255)

                                                 cTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTc                                        >  1:1781895/1‑93 (MQ=255)

                                                 cTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTc                                        >  2:678802/1‑93 (MQ=255)

                                                  tGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTcc                                       <  1:1841037/93‑1 (MQ=255)

                                                  tGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTcc                                       >  2:4751355/1‑93 (MQ=255)

                                                  tGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTcc                                       >  1:2392947/1‑93 (MQ=255)

                                                    gtgtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTg                                     >  2:659456/1‑93 (MQ=255)

                                                    gtgtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTg                                     <  1:3915334/93‑1 (MQ=255)

                                                    gtgtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTg                                     >  2:4426095/1‑93 (MQ=255)

                                                     tgtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTgg                                    <  2:228983/93‑1 (MQ=255)

                                                      gtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGa                                   <  2:6885527/93‑1 (MQ=255)

                                                      gtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGa                                   >  2:1055796/1‑93 (MQ=255)

                                                      gtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGa                                   <  2:3086722/93‑1 (MQ=255)

                                                        tACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGt                                 <  2:6965097/93‑1 (MQ=255)

                                                        tACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGt                                 >  2:5065350/1‑93 (MQ=255)

                                                          cTCCAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCg                               >  2:3359223/1‑93 (MQ=255)

                                                           cccAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGt                              <  1:2144277/92‑1 (MQ=255)

                                                             cAGAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTgg                            <  1:6458185/93‑1 (MQ=255)

                                                               gAAAACCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGAt                          <  1:1356913/93‑1 (MQ=255)

                                                                 aaaCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTa                        >  1:7804561/1‑93 (MQ=255)

                                                                   aCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGCAGCCTTTCCTGGAGTCGTGGATTaaa                      >  1:4647196/1‑93 (MQ=255)

                                                                   aCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTaaa                      >  2:3879099/1‑93 (MQ=255)

                                                                   aCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTaaa                      >  1:893975/1‑93 (MQ=255)

                                                                   aCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTaaa                      <  1:2180125/93‑1 (MQ=255)

                                                                   aCCCTGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTaaa                      <  1:1623185/93‑1 (MQ=255)

                                                                       tGCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATc                  <  2:7814145/93‑1 (MQ=255)

                                                                        gCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCa                 >  1:1025900/1‑93 (MQ=255)

                                                                        gCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCa                 <  1:2873993/93‑1 (MQ=255)

                                                                        gCTCTACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCa                 <  2:2504487/93‑1 (MQ=255)

                                                                         ctctACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAg                <  2:1741477/93‑1 (MQ=255)

                                                                          tctACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGTAAAATGGCGAAGCCTTTCCTGGAGTCGTGGTTTAAAGATCAgg               >  1:6005002/1‑93 (MQ=255)

                                                                          tctACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAgg               <  2:6670141/93‑1 (MQ=255)

                                                                          tctACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAgg               >  2:1940154/1‑93 (MQ=255)

                                                                          tctACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAgg               <  1:1890935/93‑1 (MQ=255)

                                                                           ctACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGt              >  2:5142711/1‑93 (MQ=255)

                                                                            tACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTc             >  1:3062707/1‑93 (MQ=255)

                                                                            tACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTc             >  1:1397764/1‑93 (MQ=255)

                                                                            tACGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTc             <  2:2406326/93‑1 (MQ=255)

                                                                              cGTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCgg           >  1:731021/1‑93 (MQ=255)

                                                                               gTCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGt          <  1:6595045/93‑1 (MQ=255)

                                                                               gTCGAAGGGGGGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGt          <  1:5649706/93‑1 (MQ=255)

                                                                                tCGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTa         <  1:6238451/93‑1 (MQ=255)

                                                                                 cGAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTAt        <  1:7100322/93‑1 (MQ=255)

                                                                                  gAAGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGCGGATTAAAGATCAGGTCGGTAtt       >  2:772687/1‑93 (MQ=255)

                                                                                   aaGGGGTGGGACGCCAGCTTTATCCGCAACTCGTTTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTc      <  1:7639983/93‑1 (MQ=255)

                                                                                   aaGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTc      >  1:6137699/1‑93 (MQ=255)

                                                                                   aaGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTc      <  2:7231657/93‑1 (MQ=255)

                                                                                   aaGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTAGGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTc      <  2:4048859/93‑1 (MQ=255)

                                                                                    aGGGGTGGGGCGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTcc     <  2:7128184/93‑1 (MQ=255)

                                                                                    aGGGGTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTcc     <  2:5465639/93‑1 (MQ=255)

                                                                                     ggggTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCt    <  2:6525149/93‑1 (MQ=255)

                                                                                     ggggTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCt    <  1:6886979/93‑1 (MQ=255)

                                                                                      gggTGGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTg   <  1:4391525/93‑1 (MQ=255)

AAATCTGTTTAATACGGCGCGTCGCTTCAATATGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGC  >  NC_000913/4018967‑4019146

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