BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,345,093	A→G	E113G (GAA→GGA)	yhbJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,345,093	0	A	G	84.8%	681.3 / 116.2	296	E113G (GAA→GGA)	yhbJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (15/30); new base G (126/124); total (141/155)
Fisher's exact test for biased strand distribution p-value = 3.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

Only 100 of 302 total aligned reads displayed.

TGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAATGTCCGTTCACGAGC  >  NC_000913/3345003‑3345185

tGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGgaa                                                                                            >  1:7072062/1‑93 (MQ=255)

  ttCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGt                                                                                          <  2:6020289/93‑1 (MQ=255)

  ttCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGt                                                                                          <  1:4586977/93‑1 (MQ=255)

    ccTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGc                                                                                        <  1:2411159/93‑1 (MQ=255)

    ccTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGc                                                                                        <  2:849295/93‑1 (MQ=255)

     cTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCt                                                                                       <  2:6967620/93‑1 (MQ=255)

     cTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCt                                                                                       <  2:4114248/93‑1 (MQ=255)

     cTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCt                                                                                       <  1:3856104/93‑1 (MQ=255)

      tGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTa                                                                                      <  1:1331740/93‑1 (MQ=255)

          tGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGa                                                                                  >  2:1750299/1‑93 (MQ=255)

          tGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGa                                                                                  <  2:6668306/93‑1 (MQ=255)

           gCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGAc                                                                                 <  1:6641205/93‑1 (MQ=255)

            ccGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACa                                                                                <  1:2227042/93‑1 (MQ=255)

               aCCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACaaag                                                                             <  1:4468545/93‑1 (MQ=255)

                 cGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACaaagaa                                                                           <  2:517115/93‑1 (MQ=255)

                 cGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACaaagaa                                                                           <  1:1347861/93‑1 (MQ=255)

                 cGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACaaagaa                                                                           <  1:5315929/93‑1 (MQ=255)

                  gTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACaaagaaa                                                                          >  1:1446630/1‑93 (MQ=255)

                  gTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACaaagaaa                                                                          >  2:230594/1‑93 (MQ=255)

                   tAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACaaagaaag                                                                         >  2:3462886/1‑93 (MQ=255)

                   tAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACaaagaaag                                                                         <  2:3614926/93‑1 (MQ=255)

                    aaTACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGc                                                                        <  2:1398639/93‑1 (MQ=255)

                    aaTACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGc                                                                        <  1:891928/93‑1 (MQ=255)

                     aTACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCg                                                                       >  1:3230647/1‑93 (MQ=255)

                      tACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGa                                                                      >  2:6714207/1‑93 (MQ=255)

                      tACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGa                                                                      >  2:6511257/1‑93 (MQ=255)

                      tACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGa                                                                      >  1:1806870/1‑93 (MQ=255)

                       aCCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGAt                                                                     <  2:3683469/93‑1 (MQ=255)

                        ccTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGAtt                                                                    <  1:4074154/93‑1 (MQ=255)

                        ccTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGAtt                                                                    <  2:6167810/93‑1 (MQ=255)

                         cTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGAttt                                                                   <  1:6024020/93‑1 (MQ=255)

                          ttAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTg                                                                  >  1:4851765/1‑93 (MQ=255)

                          ttAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTg                                                                  >  2:3882838/1‑93 (MQ=255)

                            aaTTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGTAAGTGCTATCGACAAAGAAAGCGATTTGCt                                                                <  1:7360593/93‑1 (MQ=255)

                            aaTTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCt                                                                >  2:1148083/1‑93 (MQ=255)

                             aTTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTg                                                               >  2:1338371/1‑93 (MQ=255)

                             aTTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTg                                                               <  1:4662537/93‑1 (MQ=255)

                              ttCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTgg                                                              >  2:5096280/1‑93 (MQ=255)

                              ttCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTgg                                                              >  1:7248341/1‑93 (MQ=255)

                               tcgtcgTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGa                                                             <  1:2965662/93‑1 (MQ=255)

                                 gtagtTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGc                                                           <  1:5448955/90‑1 (MQ=255)

                                  tcgtTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGcc                                                          <  2:2009985/93‑1 (MQ=255)

                                  tcgtTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGACGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGCAGcc                                                          >  2:4736328/1‑93 (MQ=255)

                                   cgtTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCct                                                         <  1:5165124/93‑1 (MQ=255)

                                    gtTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCctc                                                        <  2:3274106/93‑1 (MQ=255)

                                    gtTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCctc                                                        <  2:7086111/93‑1 (MQ=255)

                                    gtTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCctc                                                        <  2:303887/93‑1 (MQ=255)

                                     ttACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCctct                                                       >  2:1621861/1‑93 (MQ=255)

                                         aGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGt                                                   >  2:7209527/1‑93 (MQ=255)

                                         aGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGt                                                   <  2:2670337/93‑1 (MQ=255)

                                          gTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGtt                                                  <  1:3390986/93‑1 (MQ=255)

                                          gTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGtt                                                  >  1:163330/1‑93 (MQ=255)

                                            gACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTcg                                                <  2:3126513/93‑1 (MQ=255)

                                            gACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTcg                                                >  1:6818867/1‑93 (MQ=255)

                                             acacGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTcgc                                               <  2:4579545/93‑1 (MQ=255)

                                               nngcgcCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGa                                             <  1:140705/91‑1 (MQ=255)

                                               acGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGa                                             >  1:5521222/1‑93 (MQ=255)

                                               acGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGa                                             <  2:6223362/93‑1 (MQ=255)

                                                  cgcCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCg                                          <  1:5742489/93‑1 (MQ=255)

                                                   gcCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCgg                                         <  2:5720449/93‑1 (MQ=255)

                                                    ccGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGa                                        <  2:6836925/93‑1 (MQ=255)

                                                      gACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATc                                      >  1:262624/1‑93 (MQ=255)

                                                      gACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATc                                      >  2:5427982/1‑93 (MQ=255)

                                                       aCTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCt                                     <  2:1916912/93‑1 (MQ=255)

                                                         tGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGa                                   >  2:2454853/1‑93 (MQ=255)

                                                          gCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGAt                                  >  1:1685241/1‑93 (MQ=255)

                                                          gCATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGAt                                  <  1:7477105/93‑1 (MQ=255)

                                                          gCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGAt                                  >  1:6321203/1‑93 (MQ=255)

                                                           cATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGAtt                                 <  1:4382187/93‑1 (MQ=255)

                                                           cATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGAtt                                 <  1:4237059/93‑1 (MQ=255)

                                                           cATCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGAtt                                 <  1:6582024/93‑1 (MQ=255)

                                                           cATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGAtt                                 <  2:6253753/93‑1 (MQ=255)

                                                            aTCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTg                                <  1:442113/93‑1 (MQ=255)

                                                            aTCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTg                                <  1:4741594/93‑1 (MQ=255)

                                                            aTCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTg                                <  2:4445809/93‑1 (MQ=255)

                                                             tCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGt                               >  1:43441/1‑93 (MQ=255)

                                                             tCCGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGt                               >  1:3720074/1‑93 (MQ=255)

                                                               cGCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCg                             >  1:4151176/1‑93 (MQ=255)

                                                                gCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGa                            <  1:6902685/93‑1 (MQ=255)

                                                                gCTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGa                            <  2:4921832/93‑1 (MQ=255)

                                                                gCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGa                            <  1:2740892/93‑1 (MQ=255)

                                                                 cTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGAt                           <  1:3539804/93‑1 (MQ=255)

                                                                 cTTTCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGAt                           >  1:5032678/1‑93 (MQ=255)

                                                                  tttCCAGCAAAACCCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATa                          <  1:2336668/93‑1 (MQ=255)

                                                                  tttCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATa                          >  1:7188438/1‑93 (MQ=255)

                                                                   ttCCAGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATAc                         >  2:3567951/1‑93 (MQ=255)

                                                                       aGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCa                     >  2:4532240/1‑93 (MQ=255)

                                                                       aGCAAAAACCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCa                     <  1:1919621/93‑1 (MQ=255)

                                                                          aaaaaCCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGaa                  >  1:6482554/1‑93 (MQ=255)

                                                                           aaaaCCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGaaa                 >  2:3991622/1‑93 (MQ=255)

                                                                           aaaaCCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGaaa                 <  1:2480020/93‑1 (MQ=255)

                                                                           aaaaCCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGaaa                 <  2:7402795/93‑1 (MQ=255)

                                                                            aaaCCTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAAt                >  2:3751003/1‑93 (MQ=255)

                                                                               ccTGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAATGTc             <  2:1490010/93‑1 (MQ=255)

                                                                                 tGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAATGTCCg           >  1:848676/1‑93 (MQ=255)

                                                                                 tGTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAATGTCCg           >  1:5822958/1‑93 (MQ=255)

                                                                                  gTCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAATGTCCGt          <  2:3961577/93‑1 (MQ=255)

                                                                                   tCGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAATGTCCGtt         <  2:2589797/93‑1 (MQ=255)

                                                                                    cGCTGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAATGTCCGTTc        <  2:5318182/93‑1 (MQ=255)

                                                                                       tGGGAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAATGTCCGTTCACg     >  1:1058569/1‑93 (MQ=255)

TGTTCCTGGATGCCGACCGTAATACCTTAATTCGTCGTTACAGTGACACGCGCCGACTGCATCCGCTTTCCAGCAAAAACCTGTCGCTGGAAAGTGCTATCGACAAAGAAAGCGATTTGCTGGAGCCTCTGCGTTCGCGAGCGGATCTGATTGTCGATACCTCAGAAATGTCCGTTCACGAGC  >  NC_000913/3345003‑3345185

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