BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,455,409	G→A	R591C (CGT→TGT)	b2341 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,455,409	0	G	A	99.5%	695.7 / NA	204	R591C (CGT→TGT)	b2341	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (89/114); minor base T (0/1); total (89/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

GTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCAGGCG  >  NC_000913/2455319‑2455498

gTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACAcc                                                                                         >  1:2629154/1‑93 (MQ=255)

 ttCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACAccc                                                                                        <  1:5116325/93‑1 (MQ=255)

 ttCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACAccc                                                                                        >  2:505979/1‑93 (MQ=255)

 ttCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACAccc                                                                                        >  2:457595/1‑93 (MQ=255)

 ttCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACAccc                                                                                        <  1:2121956/93‑1 (MQ=255)

 tgcagcagCCGGCGGTGCGGAGATGCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACAccc                                                                                        <  2:4648393/91‑1 (MQ=255)

     gCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCtttt                                                                                    >  2:5632537/1‑92 (MQ=255)

     gCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTc                                                                                    <  2:689069/93‑1 (MQ=255)

     gCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTc                                                                                    <  1:5162986/93‑1 (MQ=255)

     gCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTc                                                                                    <  2:164613/93‑1 (MQ=255)

      cAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCt                                                                                   <  2:128915/93‑1 (MQ=255)

      cAACCTGCCGTGCGGAGATTCGCCCCTGGCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCt                                                                                   <  1:1213755/93‑1 (MQ=255)

       aaCCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTg                                                                                  >  2:5700433/1‑93 (MQ=255)

       aaCCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTg                                                                                  >  2:1638289/1‑93 (MQ=255)

        aCCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGa                                                                                 >  1:3717282/1‑93 (MQ=255)

        aCCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGa                                                                                 >  2:5916593/1‑93 (MQ=255)

            gcggtgcggAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCCTTCTGACCAt                                                                             >  1:2198417/1‑93 (MQ=255)

               gtgcggAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATaaa                                                                          <  2:1776223/93‑1 (MQ=255)

                ggcggAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAg                                                                         <  1:877934/92‑1 (MQ=255)

                 gcggAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGGTCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGa                                                                        <  1:4440176/93‑1 (MQ=255)

                  cggAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGAt                                                                       <  1:1382733/93‑1 (MQ=255)

                     agaTTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGa                                                                    <  2:4613245/93‑1 (MQ=255)

                      gaTTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACATTAAAGATAGaa                                                                   >  2:2484130/1‑93 (MQ=255)

                       aTTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGaaa                                                                  <  1:6118618/93‑1 (MQ=255)

                       aTTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGaaa                                                                  <  1:3877639/93‑1 (MQ=255)

                        ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATACAGATAGAAAc                                                                 >  1:4010786/1‑93 (MQ=255)

                         tCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAAcc                                                                <  2:2579465/93‑1 (MQ=255)

                         tCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAAcc                                                                <  1:1670149/93‑1 (MQ=255)

                           gCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAAcccc                                                              <  2:3364864/93‑1 (MQ=255)

                           gCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAAcccc                                                              <  2:1037032/93‑1 (MQ=255)

                           gCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAAcccc                                                              <  1:1840039/93‑1 (MQ=255)

                            ccccTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCg                                                             <  1:2775544/93‑1 (MQ=255)

                             cccTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCgg                                                            >  2:4132934/1‑93 (MQ=255)

                                tGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCa                                                         >  1:1093122/1‑93 (MQ=255)

                                tGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCa                                                         >  1:3250310/1‑93 (MQ=255)

                                  cccTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCAtt                                                       <  2:3799067/93‑1 (MQ=255)

                                  cccTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCAtt                                                       <  2:203907/93‑1 (MQ=255)

                                   ccTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCAttt                                                      >  2:3109148/1‑93 (MQ=255)

                                    cTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCAtttt                                                     <  2:2439358/93‑1 (MQ=255)

                                    cTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCAtttt                                                     >  2:151559/1‑93 (MQ=255)

                                     ttGTGTGCCAATCAGCGGGTAAATGGCGGGATCGAACTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCAttttt                                                    >  1:2002833/1‑93 (MQ=255)

                                       gtgtgCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTc                                                  <  1:1886759/93‑1 (MQ=255)

                                       gtgtgCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTc                                                  <  1:1467410/93‑1 (MQ=255)

                                        tgtgCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCt                                                 >  2:3580942/1‑93 (MQ=255)

                                          tgCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGc                                               >  1:964463/1‑93 (MQ=255)

                                             cAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCtt                                            <  2:5364591/93‑1 (MQ=255)

                                              aaTCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCttt                                           <  2:5549099/93‑1 (MQ=255)

                                               aTCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTg                                          >  1:2383510/1‑93 (MQ=255)

                                                tCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGc                                         <  1:46286/93‑1 (MQ=255)

                                                 cAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCg                                        <  2:1654299/93‑1 (MQ=255)

                                                   gCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGAt                                      <  2:109184/93‑1 (MQ=255)

                                                   gCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGAt                                      <  2:1218008/93‑1 (MQ=255)

                                                      gctAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGAtcgt                                   <  1:1630563/91‑1 (MQ=255)

                                                        tAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGAtcgtcg                                 <  2:523758/93‑1 (MQ=255)

                                                        tAAATGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGAtcgtcg                                 >  1:1301046/1‑93 (MQ=255)

                                                         aaaTGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATcgtcgt                                <  1:2319891/93‑1 (MQ=255)

                                                         aaaTGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTACTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATcgtcgt                                <  1:2500110/93‑1 (MQ=255)

                                                           aTGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTc                              <  1:1269664/93‑1 (MQ=255)

                                                            tGGCGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCa                             >  1:4908100/1‑93 (MQ=255)

                                                               cGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCaaaa                          <  2:2138347/93‑1 (MQ=255)

                                                               cGGGATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCaaaa                          <  2:5040763/93‑1 (MQ=255)

                                                                gggATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAAt                         >  1:3361639/1‑93 (MQ=255)

                                                                gggATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAAt                         <  1:350382/93‑1 (MQ=255)

                                                                  gATCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTg                       >  2:6242880/1‑93 (MQ=255)

                                                                    tCGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTgaa                     <  1:2855717/93‑1 (MQ=255)

                                                                     cGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTgaag                    <  2:5990464/93‑1 (MQ=255)

                                                                     cGACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTgaag                    <  1:422986/93‑1 (MQ=255)

                                                                      gACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTgaaga                   <  2:4420151/93‑1 (MQ=255)

                                                                      gACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTgaaga                   >  2:3193024/1‑93 (MQ=255)

                                                                      gACCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTgaaga                   >  2:886612/1‑93 (MQ=255)

                                                                       aCCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTgaagaa                  <  2:4916024/93‑1 (MQ=255)

                                                                       aCCTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTgaagaa                  >  2:3286320/1‑93 (MQ=255)

                                                                         cTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAaac                >  2:765263/1‑93 (MQ=255)

                                                                         cTGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAaac                <  1:5528824/93‑1 (MQ=255)

                                                                         cTGTTTGTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAaac                <  1:3945861/93‑1 (MQ=255)

                                                                          tGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAaaca               >  2:3549356/1‑93 (MQ=255)

                                                                          tGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAaaca               >  1:4266741/1‑93 (MQ=255)

                                                                          tGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAaaca               <  1:5403437/93‑1 (MQ=255)

                                                                          tGTTTTTTGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAaaca               >  1:6013051/1‑93 (MQ=255)

                                                                            ttttttGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAaacaac             >  2:4983309/1‑93 (MQ=255)

                                                                            ttttttGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAaacaac             <  2:1441614/93‑1 (MQ=255)

                                                                             tttttGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAacaaca            <  1:2491987/93‑1 (MQ=255)

                                                                              ttttGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACAt           <  2:2812716/93‑1 (MQ=255)

                                                                              ttttGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACAt           <  2:3866032/93‑1 (MQ=255)

                                                                              ttttGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACAt           >  1:5807207/1‑93 (MQ=255)

                                                                               tttGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACAtt          >  1:4970593/1‑93 (MQ=255)

                                                                               tttGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACAtt          <  1:5843718/93‑1 (MQ=255)

                                                                                ttGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACAttt         <  1:3415042/93‑1 (MQ=255)

                                                                                 tGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTg        <  1:5305161/93‑1 (MQ=255)

                                                                                 tGCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTg        <  1:2994185/93‑1 (MQ=255)

                                                                                  gCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGc       >  2:3371334/1‑93 (MQ=255)

                                                                                  gCTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATAGTCGTTCAAAATTGAAGAAACAACATTTGc       >  1:6308348/1‑93 (MQ=255)

                                                                                   cTTTTACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCa      <  2:3653684/93‑1 (MQ=255)

                                                                                    ttttACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCAg     <  2:4774924/93‑1 (MQ=255)

                                                                                    ttttACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCAg     <  1:1019848/93‑1 (MQ=255)

                                                                                    ttttACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCAg     <  2:450781/93‑1 (MQ=255)

                                                                                     tttACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCAgg    <  1:2057904/93‑1 (MQ=255)

                                                                                      ttACACCCTTTCTGACCATAAAGATAGAACCCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCAGgc   <  1:4327006/93‑1 (MQ=255)

                                                                                      ttACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCAGgc   <  1:3789948/93‑1 (MQ=255)

                                                                                       tACACCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCAGgcg  <  1:1718879/93‑1 (MQ=255)

GTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTTGCGATCGTCGTTCAAAATTGAAGAAACAACATTTGCAGGCG  >  NC_000913/2455319‑2455498

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