BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,406,298	A→G	Q689Q (CAA→CAG)	b4179 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,406,298	0	A	G	99.7%	1013.9 / ‑5.1	321	Q689Q (CAA→CAG)	b4179	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (159/161); total (159/162)
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.51e-01

Only 100 of 334 total aligned reads displayed.

TTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAA  >  NC_000913/4406209‑4406387

ttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAgcgc                                                                                        >  1:7176945/1‑93 (MQ=255)

ttcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAgcgc                                                                                        <  2:507540/93‑1 (MQ=255)

 tcGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGcc                                                                                       <  1:16087/93‑1 (MQ=255)

  cgtccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTCGGGAAGCGCCt                                                                                      >  2:7531800/1‑93 (MQ=255)

  cgtccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCt                                                                                      <  2:4843347/93‑1 (MQ=255)

  cgtccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCt                                                                                      >  2:6792594/1‑93 (MQ=255)

  cgtccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCt                                                                                      <  2:247613/93‑1 (MQ=255)

   gtccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTg                                                                                     >  1:1603350/1‑93 (MQ=255)

   gtccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTg                                                                                     >  1:857029/1‑93 (MQ=255)

    tccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGa                                                                                    <  1:1408378/93‑1 (MQ=255)

    tccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGa                                                                                    <  2:404807/93‑1 (MQ=255)

    tccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGa                                                                                    <  1:2589026/93‑1 (MQ=255)

    tccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGa                                                                                    <  2:6478594/93‑1 (MQ=255)

    tccgtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGa                                                                                    >  2:6142766/1‑93 (MQ=255)

       gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 <  1:3849907/93‑1 (MQ=255)

       gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 <  1:4712262/93‑1 (MQ=255)

       gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 <  2:367035/93‑1 (MQ=255)

       gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 <  2:578941/93‑1 (MQ=255)

       gtcTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATgg                                                                                 <  1:6658186/93‑1 (MQ=255)

          tGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggggg                                                                              <  2:1940480/93‑1 (MQ=255)

          tGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATggggg                                                                              <  2:2029218/93‑1 (MQ=255)

          tGGACGACTTGTTCATTGATGGGCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGTGCAGCGCCTGATggggg                                                                              >  2:4156009/1‑93 (MQ=255)

           ggACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGa                                                                             >  2:7289636/1‑93 (MQ=255)

                  ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGc                                                                      >  1:7585125/1‑93 (MQ=255)

                  ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGc                                                                      <  1:7925208/93‑1 (MQ=255)

                  ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGc                                                                      >  2:3952590/1‑93 (MQ=255)

                  ttGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCCGATGGGGGAATCCAGc                                                                      <  1:4971914/93‑1 (MQ=255)

                  ctGTTGATTGATGGTCGGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGc                                                                      <  1:6707890/92‑1 (MQ=255)

                   tGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCg                                                                     >  1:4309079/1‑93 (MQ=255)

                   tGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCg                                                                     <  1:6873579/93‑1 (MQ=255)

                    gTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCgg                                                                    >  2:2929444/1‑93 (MQ=255)

                        aTTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAg                                                                >  1:987959/1‑93 (MQ=255)

                         ttGATGGTCTGGTCCCTGCCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGa                                                               <  1:2382205/93‑1 (MQ=255)

                         ttGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGa                                                               >  2:7675443/1‑93 (MQ=255)

                          tGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGAc                                                              <  1:4460141/93‑1 (MQ=255)

                           gATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACt                                                             >  2:7426512/1‑93 (MQ=255)

                             tggtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTa                                                           <  1:6917043/93‑1 (MQ=255)

                              ggtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTAt                                                          >  1:1121961/1‑93 (MQ=255)

                               gtctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATc                                                         <  2:1058486/93‑1 (MQ=255)

                                tctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCg                                                        >  2:3735619/1‑93 (MQ=255)

                                tctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCg                                                        >  2:3881063/1‑93 (MQ=255)

                                 ctggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGc                                                       <  2:6424258/93‑1 (MQ=255)

                                  tggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGcc                                                      <  1:7265904/93‑1 (MQ=255)

                                   ggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCt                                                     <  2:5681687/93‑1 (MQ=255)

                                   ggtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCt                                                     >  2:1774556/1‑93 (MQ=255)

                                    gtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTg                                                    <  2:6213381/93‑1 (MQ=255)

                                    gtcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTg                                                    <  1:1202335/93‑1 (MQ=255)

                                     tcCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTgg                                                   <  2:4571478/93‑1 (MQ=255)

                                        aTGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCg                                                >  1:5502219/1‑93 (MQ=255)

                                         tGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGa                                               >  2:1613279/1‑93 (MQ=255)

                                          gTCTCTTCGCTGGACAATGACTATTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                              <  2:4334047/93‑1 (MQ=255)

                                          gTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGAt                                              >  1:7489091/1‑93 (MQ=255)

                                           tctcTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATc                                             >  2:4907387/1‑93 (MQ=255)

                                           tctcTTCGCTGGACAATGACTACTATCGCTTTGACCAGGGAGGGCAGCGCCTGATGGGGGAATCCAGCGGCAAGAGTTATCGCCTGGGCGacc                                             >  2:1878041/1‑91 (MQ=255)

                                            ctctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcg                                            >  2:1355538/1‑93 (MQ=255)

                                             tctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcgg                                           >  2:2619704/1‑92 (MQ=255)

                                             tctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcgc                                           >  2:7615177/1‑93 (MQ=255)

                                              ctTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATcgcg                                          >  2:816961/1‑93 (MQ=255)

                                                tCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTg                                        <  1:3968454/93‑1 (MQ=255)

                                                tCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGAGGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTg                                        <  2:5484317/93‑1 (MQ=255)

                                                   cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGaa                                     <  1:4864709/93‑1 (MQ=255)

                                                   cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGaa                                     <  1:6057363/93‑1 (MQ=255)

                                                    tGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTGTCGCCTGGGCGATCGCGTGGAAg                                    >  2:4295851/1‑93 (MQ=255)

                                                    tGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAg                                    >  2:4361712/1‑93 (MQ=255)

                                                       aCAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTc                                 <  1:5083307/93‑1 (MQ=255)

                                                       aCAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTc                                 >  1:7049091/1‑93 (MQ=255)

                                                        cAATGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcg                                >  1:1672523/1‑93 (MQ=255)

                                                         aaTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGttggg                               >  2:6099213/1‑90 (MQ=255)

                                                         aaTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgc                               <  2:1558172/93‑1 (MQ=255)

                                                         aaTGACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTcgc                               >  2:3434637/1‑93 (MQ=255)

                                                            gACTACTATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGGTCGCGTc                            <  2:7916323/93‑1 (MQ=255)

                                                             actactATCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCg                           >  1:2952127/1‑93 (MQ=255)

                                                               tactaTCGCTTTGACCAGGCAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATTGCCTGGGCGATCGCGTGGAAGTTCGCGTCgca                         >  2:6761628/1‑91 (MQ=255)

                                                                actaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAg                        >  1:5129132/1‑93 (MQ=255)

                                                                actaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAg                        <  1:569258/93‑1 (MQ=255)

                                                                 ctaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGc                       >  2:5265278/1‑93 (MQ=255)

                                                                 ctaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGc                       >  2:5877258/1‑93 (MQ=255)

                                                                 ctaTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGc                       >  2:7476329/1‑93 (MQ=255)

                                                                  taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCg                      <  1:3058650/93‑1 (MQ=255)

                                                                  taTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCg                      >  1:7727810/1‑93 (MQ=255)

                                                                   aTCGCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCgg                     >  2:234106/1‑93 (MQ=255)

                                                                      gCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTa                  <  1:4070259/93‑1 (MQ=255)

                                                                      gCTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTa                  >  1:364375/1‑93 (MQ=255)

                                                                       cTTTGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTaa                 <  2:4760880/93‑1 (MQ=255)

                                                                        tttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAat                >  1:4191665/1‑93 (MQ=255)

                                                                         ttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAata               >  2:5435500/1‑93 (MQ=255)

                                                                         ttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAata               <  1:5384296/93‑1 (MQ=255)

                                                                         ttGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAata               <  2:7812581/93‑1 (MQ=255)

                                                                          tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat              <  2:3589283/93‑1 (MQ=255)

                                                                          tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat              >  1:321702/1‑93 (MQ=255)

                                                                          tGACCAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAatat              <  2:2136273/93‑1 (MQ=255)

                                                                             ccAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGa           <  1:5245591/93‑1 (MQ=255)

                                                                             ccAGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGa           <  1:6223059/93‑1 (MQ=255)

                                                                               aGGTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACg         <  1:228331/93‑1 (MQ=255)

                                                                                ggTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa        >  2:6555370/1‑93 (MQ=255)

                                                                                ggTAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGa        >  1:3319423/1‑93 (MQ=255)

                                                                                  tAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgc      >  2:3301925/1‑93 (MQ=255)

                                                                                  tAGGGCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgc      <  2:4518709/93‑1 (MQ=255)

                                                                                    gggCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAgcgc    >  1:2111286/1‑93 (MQ=255)

                                                                                      gCAGCGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCaa  >  1:6424021/1‑93 (MQ=255)

TTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCAGACTTATCGCCTGGGCGATCGCGTGGAAGTTCGCGTCGAAGCGGTTAATATGGACGAGCGCAA  >  NC_000913/4406209‑4406387

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