BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,751,091	G→C	A426P (GCT→CCT)	recN →	RecN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,751,091	0	G	C	99.7%	1026.6 / NA	311	A426P (GCT→CCT)	recN	RecN
Reads supporting (aligned to +/- strand): ref base G (1/0); new base C (159/151); total (160/151)
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 320 total aligned reads displayed.

ATGTTAAATTTGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGGCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAACCCCGGC  >  NC_000913/2751002‑2751179

   ttAAATTTGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTAttc                                                                                    >  1:1546955/1‑92 (MQ=255)

    tAAATTTGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGc                                                                                   <  1:4315367/93‑1 (MQ=255)

        tttGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCaaa                                                                               <  2:373914/93‑1 (MQ=255)

        tttGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCaaa                                                                               <  1:5644363/93‑1 (MQ=255)

         ttGACGAGCATCGCCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAGCCCAGGTCAGCCAATGCAGCCTATTGCCAAAg                                                                              <  1:3751246/93‑1 (MQ=255)

         ttGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAg                                                                              <  2:3312267/93‑1 (MQ=255)

          tGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGt                                                                             <  1:9756647/93‑1 (MQ=255)

           gACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTc                                                                            >  2:9384532/1‑93 (MQ=255)

            aCGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCg                                                                           >  2:5210762/1‑93 (MQ=255)

                gCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATc                                                                       <  1:4968088/93‑1 (MQ=255)

                gCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATc                                                                       <  2:4592524/93‑1 (MQ=255)

                   tCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCgg                                                                    >  1:9986203/1‑93 (MQ=255)

                     aCCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtg                                                                  >  2:6810315/1‑93 (MQ=255)

                     aCCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtg                                                                  <  2:4386806/93‑1 (MQ=255)

                     aCCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtg                                                                  >  1:443799/1‑93 (MQ=255)

                      ccTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtgg                                                                 <  2:1549785/93‑1 (MQ=255)

                      ccTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtgg                                                                 >  2:1323773/1‑93 (MQ=255)

                      ccTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtgg                                                                 >  1:7548291/1‑93 (MQ=255)

                      ccTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtgg                                                                 <  2:9123325/93‑1 (MQ=255)

                       cTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtggt                                                                >  1:2756105/1‑93 (MQ=255)

                        tGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGgtggtg                                                               >  1:3919884/1‑93 (MQ=255)

                        tGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGgtggtg                                                               >  1:8609111/1‑93 (MQ=255)

                        tGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAACTCGCATCCGgtggtg                                                               >  2:2722739/1‑93 (MQ=255)

                         gggCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGa                                                              >  2:22358/1‑93 (MQ=255)

                         gggCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGa                                                              <  1:8426722/93‑1 (MQ=255)

                         gggCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGa                                                              >  1:2274465/1‑93 (MQ=255)

                          ggCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGaa                                                             <  1:6762027/93‑1 (MQ=255)

                           gcgcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAAt                                                            <  2:6995368/93‑1 (MQ=255)

                            cgcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAAtt                                                           <  2:2025419/93‑1 (MQ=255)

                            cgcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGAAGCCTATTGCCAAAGTCGCATCCGGTGGTGagtg                                                           >  1:6993554/1‑90 (MQ=255)

                             gcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTg                                                          <  2:5504075/93‑1 (MQ=255)

                             gcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTg                                                          >  2:1080511/1‑93 (MQ=255)

                             gcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTg                                                          >  1:6287027/1‑93 (MQ=255)

                                gACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTcc                                                       >  2:7432448/1‑93 (MQ=255)

                                gACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTcc                                                       <  1:7729343/93‑1 (MQ=255)

                                 aCGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTccc                                                      <  1:454155/93‑1 (MQ=255)

                                    gTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCa                                                   >  1:8112787/1‑93 (MQ=255)

                                    gTGCCGATCGTATTGAGTTTCGGGGAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGTCAAAGTCGCATCCGGTGGCGAATTGTCCCGCa                                                   >  2:9652011/1‑93 (MQ=255)

                                     tGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCAt                                                  <  1:195529/93‑1 (MQ=255)

                                     tGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCAt                                                  <  1:8700620/93‑1 (MQ=255)

                                      gCCGATCGTATTGAGTTTGGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATc                                                 <  2:6277151/93‑1 (MQ=255)

                                      gCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATc                                                 <  2:8005379/93‑1 (MQ=255)

                                         gATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCa                                              <  2:1813329/93‑1 (MQ=255)

                                         gATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCa                                              <  1:9346873/93‑1 (MQ=255)

                                          aTCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCAc                                             >  1:1272998/1‑93 (MQ=255)

                                           tCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACt                                            >  2:4523568/1‑93 (MQ=255)

                                           tCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACt                                            >  2:1186775/1‑93 (MQ=255)

                                           tCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACt                                            >  1:178025/1‑93 (MQ=255)

                                           tCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACt                                            >  1:1117645/1‑93 (MQ=255)

                                            cGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTg                                           <  1:649087/93‑1 (MQ=255)

                                            cGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTg                                           <  1:7235382/93‑1 (MQ=255)

                                             gTATTGAGTTTCTGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTgg                                          <  1:5445379/93‑1 (MQ=255)

                                              tATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGc                                         >  2:2982708/1‑93 (MQ=255)

                                              tATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGc                                         <  1:2604830/93‑1 (MQ=255)

                                              tATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGc                                         >  2:8176818/1‑93 (MQ=255)

                                               aTTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCa                                        >  1:5998777/1‑93 (MQ=255)

                                                 tGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAAt                                      >  1:2938499/1‑93 (MQ=255)

                                                 tGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAAt                                      <  1:5502914/93‑1 (MQ=255)

                                                 tGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACCGGCCAt                                      >  2:5700108/1‑93 (MQ=255)

                                                  gAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATc                                     <  2:5356487/93‑1 (MQ=255)

                                                  gAGTTTCGGGGAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATc                                     >  1:6562669/1‑93 (MQ=255)

                                                   aGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATcc                                    <  2:4807706/93‑1 (MQ=255)

                                                    gTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCa                                   <  2:7538942/93‑1 (MQ=255)

                                                    gTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCa                                   <  2:3154564/93‑1 (MQ=255)

                                                     tttCGGGTAACAACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAg                                  >  2:8528821/1‑93 (MQ=255)

                                                       tCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGt                                <  1:7903907/93‑1 (MQ=255)

                                                          ggTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtcac                             >  1:1991302/1‑92 (MQ=255)

                                                           gTAACCACCAAACCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtcatc                            >  1:8160128/1‑93 (MQ=255)

                                                            tAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtcatca                           <  1:4167066/93‑1 (MQ=255)

                                                            tAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtcatca                           <  1:9501927/93‑1 (MQ=255)

                                                            tAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtcatca                           <  1:920641/93‑1 (MQ=255)

                                                            tAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtcatca                           <  1:8263423/93‑1 (MQ=255)

                                                                caccaACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGgc                       >  1:337231/1‑93 (MQ=255)

                                                                 cccaACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGgcg                      <  2:1439426/92‑1 (MQ=255)

                                                                 accaACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGgcg                      >  1:7477812/1‑93 (MQ=255)

                                                                 accaACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGgcg                      >  2:6506650/1‑93 (MQ=255)

                                                                   caACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGcgcg                    <  2:889532/93‑1 (MQ=255)

                                                                    aaCCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGt                   <  1:2891492/93‑1 (MQ=255)

                                                                      cccAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTaa                 >  2:5916381/1‑93 (MQ=255)

                                                                        cAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTaaaa               >  1:3613564/1‑93 (MQ=255)

                                                                        cAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTaaaa               >  2:3403586/1‑93 (MQ=255)

                                                                        cAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTaaaa               <  1:2222484/93‑1 (MQ=255)

                                                                         aGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAAt              >  1:1178448/1‑93 (MQ=255)

                                                                           gTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATgg            <  1:356721/93‑1 (MQ=255)

                                                                           gTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATgg            <  2:9183456/93‑1 (MQ=255)

                                                                           gTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATgg            >  2:6416494/1‑93 (MQ=255)

                                                                            tCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGa           <  2:7063731/93‑1 (MQ=255)

                                                                              aGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGaaa         <  2:7803523/93‑1 (MQ=255)

                                                                               gCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAAc        <  2:4894861/93‑1 (MQ=255)

                                                                               gCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAAc        <  1:304505/93‑1 (MQ=255)

                                                                                ccAATGCAGCCTATTGCCAAAGTTGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAAcc       >  2:6095513/1‑93 (MQ=255)

                                                                                ccAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAAcc       >  2:279520/1‑93 (MQ=255)

                                                                                  aaTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAAcccc     >  1:606736/1‑93 (MQ=255)

                                                                                  aaTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAAcccc     >  2:8765169/1‑93 (MQ=255)

                                                                                  aaTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAAcccc     <  1:8041189/93‑1 (MQ=255)

                                                                                  aaTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAAcccc     <  1:3194903/93‑1 (MQ=255)

                                                                                   aTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAACCCCg    <  2:2546923/93‑1 (MQ=255)

                                                                                   aTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAACCCCg    <  1:9095642/93‑1 (MQ=255)

                                                                                   aTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAACCCCg    <  2:5422536/93‑1 (MQ=255)

                                                                                    tGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGAAAACCCCgg   >  1:5615885/1‑93 (MQ=255)

ATGTTAAATTTGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGGCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGTAAAATGGAAACCCCGGC  >  NC_000913/2751002‑2751179

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