BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,117,607	A→G	G31G (GGT→GGC)	yceI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,117,607	0	A	G	99.7%	1035.1 / NA	322	G31G (GGT→GGC)	yceI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (142/179); minor base T (0/1); total (142/180)
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.99e-01

Only 100 of 331 total aligned reads displayed.

AAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGACCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCTTTTTTT  >  NC_000913/1117517‑1117696

aaaGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGcc                                                                                         <  1:10037338/93‑1 (MQ=255)

  aGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCctt                                                                                       >  2:3567459/1‑93 (MQ=255)

    gTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCcttct                                                                                     <  1:4846397/93‑1 (MQ=255)

     tAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAGCAAAGGCGTGCTGGCcttctt                                                                                    >  2:4370654/1‑93 (MQ=255)

     tAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCcttctt                                                                                    <  1:4917873/93‑1 (MQ=255)

     tAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCcttctt                                                                                    >  1:7916345/1‑93 (MQ=255)

     tAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCcttctt                                                                                    <  1:7655893/93‑1 (MQ=255)

          gTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCa                                                                               <  2:6262437/93‑1 (MQ=255)

          gTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCa                                                                               <  1:4374582/93‑1 (MQ=255)

           tACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCaa                                                                              <  2:378443/93‑1 (MQ=255)

           tACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCaa                                                                              <  2:5207072/93‑1 (MQ=255)

           tACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCaa                                                                              <  2:6961826/93‑1 (MQ=255)

           tACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCaa                                                                              <  1:5348560/93‑1 (MQ=255)

             ccGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAAtt                                                                            <  2:8744850/93‑1 (MQ=255)

               gTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAAtttt                                                                          >  2:2608642/1‑93 (MQ=255)

               gTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAAtttt                                                                          >  1:5214052/1‑93 (MQ=255)

                 cGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGt                                                                        <  1:2746854/93‑1 (MQ=255)

                    aaTCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAAt                                                                     >  2:6632632/1‑93 (MQ=255)

                    aaTCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAAt                                                                     >  2:9545731/1‑93 (MQ=255)

                    aaTCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAAt                                                                     <  1:249876/93‑1 (MQ=255)

                     aTCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATc                                                                    <  2:285139/93‑1 (MQ=255)

                     aTCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATc                                                                    >  1:6811898/1‑93 (MQ=255)

                      tCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATcg                                                                   <  1:6400033/93‑1 (MQ=255)

                      tCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATcg                                                                   >  1:397397/1‑93 (MQ=255)

                      tCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATcg                                                                   >  1:1778186/1‑93 (MQ=255)

                      tCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATcg                                                                   <  1:731260/93‑1 (MQ=255)

                       cTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATctg                                                                  >  1:1555205/1‑91 (MQ=255)

                       cTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATcgg                                                                  <  2:234653/93‑1 (MQ=255)

                              ggTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAAc                                                           >  1:1325228/1‑93 (MQ=255)

                               gTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAAcc                                                          >  2:7989252/1‑93 (MQ=255)

                                tGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCg                                                         >  1:7045963/1‑93 (MQ=255)

                                  ccGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCt                                                       >  1:7424254/1‑93 (MQ=255)

                                    gTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGa                                                     >  2:1335910/1‑93 (MQ=255)

                                     tataACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGaa                                                    >  1:8999077/1‑93 (MQ=255)

                                     tataACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGaa                                                    >  1:9375827/1‑93 (MQ=255)

                                      ataACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAAc                                                   <  2:4205999/93‑1 (MQ=255)

                                       taACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAAcc                                                  >  2:4488513/1‑93 (MQ=255)

                                       taACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAAcc                                                  >  1:8500424/1‑93 (MQ=255)

                                        aaCCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCg                                                 <  2:7470717/93‑1 (MQ=255)

                                        aaCCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCg                                                 >  1:5624874/1‑93 (MQ=255)

                                         aCCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCgg                                                <  2:329477/93‑1 (MQ=255)

                                         aCCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCgg                                                <  2:309050/93‑1 (MQ=255)

                                          ccAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGc                                               <  2:5885964/93‑1 (MQ=255)

                                            aGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCag                                             >  2:7336348/1‑93 (MQ=255)

                                            aGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCag                                             <  1:1451184/93‑1 (MQ=255)

                                             gCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCaga                                            >  2:2986618/1‑93 (MQ=255)

                                               tataGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAgaga                                          >  2:1522219/1‑93 (MQ=255)

                                                ataGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGaa                                         <  1:9131097/93‑1 (MQ=255)

                                                ataGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGaa                                         >  2:9835948/1‑93 (MQ=255)

                                                 taGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAAc                                        >  2:6777792/1‑93 (MQ=255)

                                                  aGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACa                                       <  2:841123/93‑1 (MQ=255)

                                                  aGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACa                                       <  2:1426404/93‑1 (MQ=255)

                                                   gCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACAt                                      <  1:8426610/93‑1 (MQ=255)

                                                    ccAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATc                                     <  2:4457944/93‑1 (MQ=255)

                                                     cAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCa                                    <  2:2172225/93‑1 (MQ=255)

                                                      aaGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAg                                   >  2:5380823/1‑93 (MQ=255)

                                                      aaGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAg                                   <  2:8153331/93‑1 (MQ=255)

                                                        ggTGCTGGATGCGGAAATTAACAGAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAggg                                 <  2:5520921/93‑1 (MQ=255)

                                                          tGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGAc                               >  1:3723709/1‑93 (MQ=255)

                                                          tGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGAc                               >  1:8277662/1‑93 (MQ=255)

                                                          tGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGAc                               >  2:7664995/1‑93 (MQ=255)

                                                          tGCTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGAc                               <  1:1978415/93‑1 (MQ=255)

                                                            cTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGAcgc                             <  2:5968637/93‑1 (MQ=255)

                                                            cTGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGAcgc                             <  1:9650506/93‑1 (MQ=255)

                                                             tGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTTAAGCGGCAGAGAACATCAGGGAcgcg                            >  2:4791926/1‑93 (MQ=255)

                                                             tGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGAcgcg                            >  2:5126717/1‑93 (MQ=255)

                                                             tGGATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGAcgcg                            >  2:10018301/1‑93 (MQ=255)

                                                              ggATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGa                           >  2:8220168/1‑93 (MQ=255)

                                                               gATGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGAAACCGCTGAACCGGCAGAGAACATCAGGGACGCGaa                          >  2:8905087/1‑93 (MQ=255)

                                                                 tGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAgg                        <  2:8547597/93‑1 (MQ=255)

                                                                 tGCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAgg                        <  1:9733824/93‑1 (MQ=255)

                                                                  gCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGt                       <  2:389772/93‑1 (MQ=255)

                                                                  gCGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGt                       <  1:7184201/93‑1 (MQ=255)

                                                                   cGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGtt                      <  2:1977444/93‑1 (MQ=255)

                                                                   cGGAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGtt                      <  2:5817887/93‑1 (MQ=255)

                                                                    ggAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTa                     <  1:3277276/93‑1 (MQ=255)

                                                                     gAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTaa                    >  1:886324/1‑93 (MQ=255)

                                                                     gAAATTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGAAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTaa                    <  1:7116298/93‑1 (MQ=255)

                                                                       aaTTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAAc                  >  2:9293770/1‑93 (MQ=255)

                                                                       aaTTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAAc                  <  1:3084272/93‑1 (MQ=255)

                                                                       aaTTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAAc                  >  1:4622163/1‑93 (MQ=255)

                                                                        aTTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAAcc                 >  2:9093156/1‑93 (MQ=255)

                                                                        aTTAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAAcc                 <  1:612037/93‑1 (MQ=255)

                                                                         ttAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCa                <  2:9969384/93‑1 (MQ=255)

                                                                          tAACAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCaa               <  1:4158846/93‑1 (MQ=255)

                                                                            aacaaGGCGGGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGc             <  2:7900052/90‑1 (MQ=255)

                                                                            aCAAAGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGc             <  1:7991092/93‑1 (MQ=255)

                                                                              aaaGGCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAg           <  2:633230/93‑1 (MQ=255)

                                                                                 ggCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAggat        >  1:3409907/1‑91 (MQ=255)

                                                                                 ggCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCt        <  1:6714821/93‑1 (MQ=255)

                                                                                  gCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCtt       >  1:3522529/1‑93 (MQ=255)

                                                                                  gCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCtt       <  1:6663576/93‑1 (MQ=255)

                                                                                  gCGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCtt       <  1:8817795/93‑1 (MQ=255)

                                                                                   tgTGCTGGGCTTCTTTGTCAATTTTGTAAGCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCttt      <  1:5509367/92‑1 (MQ=255)

                                                                                   cGTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCttt      >  2:1677653/1‑93 (MQ=255)

                                                                                    gTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCtttt     >  1:5057079/1‑93 (MQ=255)

                                                                                    gTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCtttt     >  1:5922925/1‑93 (MQ=255)

                                                                                    gTGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCtttt     >  2:9693412/1‑93 (MQ=255)

                                                                                     tGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTATACCAAGCAGGCttttt    >  1:9942408/1‑93 (MQ=255)

                                                                                     tGCTGGCCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCttttt    >  1:2110839/1‑93 (MQ=255)

AAAGGTAAAAGTACCGTCGAAATCTTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTGGATGCGGAAATTAACAAAGGCGTGCTGACCTTCTTTGTCAATTTTGTAATCGGCGGCAACCGCTGAACCGGCAGAGAACATCAGGGACGCGAAGGTTAAACCAAGCAGGCTTTTTTT  >  NC_000913/1117517‑1117696

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: `‑`