BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,458,704	C→T	S70S (AGC→AGT)	b1396 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,458,704	0	C	T	99.1%	397.5 / NA	111	S70S (AGC→AGT)	b1396	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (54/56); minor base G (1/0); total (55/56)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

CACTGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGG  >  NC_000913/1458614‑1458792

cACTGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCa                                                                                        >  2:3837247/1‑93 (MQ=255)

 aCTGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAg                                                                                       <  2:1598522/93‑1 (MQ=255)

 aCTGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAg                                                                                       <  1:2120940/93‑1 (MQ=255)

  cTGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAgg                                                                                      <  2:1300601/93‑1 (MQ=255)

  cTGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAgg                                                                                      >  2:3270159/1‑93 (MQ=255)

  cTGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGGCAgg                                                                                      >  2:2411858/1‑93 (MQ=255)

   tGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAggg                                                                                     <  1:1007508/93‑1 (MQ=255)

   tGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAggg                                                                                     <  2:4068560/93‑1 (MQ=255)

   tGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAggg                                                                                     <  2:2737000/93‑1 (MQ=255)

    gCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAgggg                                                                                    >  1:3649875/1‑93 (MQ=255)

    gCACAAATGCTTAACGGTCATCAAAGGTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAgggg                                                                                    <  2:1426276/93‑1 (MQ=255)

     cacaAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGc                                                                                   <  2:3140356/93‑1 (MQ=255)

     cacaAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGc                                                                                   <  1:3928671/93‑1 (MQ=255)

      acaAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTGCGCCTGCAATAGTCAGGGGCt                                                                                  <  1:1582995/93‑1 (MQ=255)

      acaAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCt                                                                                  <  2:3512698/93‑1 (MQ=255)

         aaTGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGc                                                                               >  2:3255486/1‑93 (MQ=255)

           tGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGccg                                                                             >  2:1039681/1‑91 (MQ=255)

           tGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAg                                                                             >  1:1804930/1‑93 (MQ=255)

            gCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGc                                                                            >  2:1693473/1‑93 (MQ=255)

              ttAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCg                                                                          >  1:2280641/1‑93 (MQ=255)

                aaCGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTc                                                                        <  1:2075907/93‑1 (MQ=255)

                   ggTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCt                                                                     <  2:1193235/93‑1 (MQ=255)

                   ggTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCt                                                                     >  2:4443797/1‑93 (MQ=255)

                    gTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCtt                                                                    >  1:3216227/1‑93 (MQ=255)

                    gTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCtt                                                                    <  1:1341650/93‑1 (MQ=255)

                      catcaAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTct                                                                  <  1:2772114/93‑1 (MQ=255)

                       atcaAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTctg                                                                 >  2:3909403/1‑93 (MQ=255)

                       atcaAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTctg                                                                 >  1:1812349/1‑93 (MQ=255)

                       atcaAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTctg                                                                 >  2:2850612/1‑93 (MQ=255)

                        tcaAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTctgc                                                                >  1:2333405/1‑93 (MQ=255)

                         caAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTctgcc                                                               <  2:4666140/93‑1 (MQ=255)

                         caAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTctgcc                                                               <  1:2200529/93‑1 (MQ=255)

                         caAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTctgcc                                                               <  1:2424627/93‑1 (MQ=255)

                           aaGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTctgcctg                                                             >  1:1497411/1‑93 (MQ=255)

                             gTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCa                                                           <  2:2595769/93‑1 (MQ=255)

                             gTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCa                                                           >  1:2267319/1‑93 (MQ=255)

                              ttGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCAc                                                          <  2:562386/93‑1 (MQ=255)

                              ttGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCAc                                                          <  2:568079/93‑1 (MQ=255)

                                 ccATGGCGGGGAGCTATTTTCACTGCCTGATACTGCCTTTTCGTACGCCTGCAATGGTCAGGGGCTGGCAGCAGTCGCTTCTCCCTACACgat                                                       <  2:2835026/93‑1 (MQ=255)

                                 ccACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACgat                                                       <  1:2562499/93‑1 (MQ=255)

                                     ggcggGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACgattgat                                                   <  1:3375599/93‑1 (MQ=255)

                                     ggcggGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACgattgat                                                   >  1:4004458/1‑93 (MQ=255)

                                     ggcggGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACgattgat                                                   <  1:1751928/93‑1 (MQ=255)

                                        gggCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGAtttt                                                <  2:4350674/93‑1 (MQ=255)

                                        gggCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGAtttt                                                <  2:1377888/93‑1 (MQ=255)

                                          gCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGAtttttt                                              <  2:3053184/93‑1 (MQ=255)

                                          gCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGAtttttt                                              <  2:4233068/93‑1 (MQ=255)

                                          gCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGAtttttt                                              <  1:1187259/93‑1 (MQ=255)

                                            aGCTATTTTCACTGGCTGATTCTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGc                                            <  1:2871048/93‑1 (MQ=255)

                                            aGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGCTTTTTTGc                                            <  2:4558248/93‑1 (MQ=255)

                                            aGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGc                                            <  1:4106573/93‑1 (MQ=255)

                                             gCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGGTTTTTTGCg                                           >  1:3213721/1‑93 (MQ=255)

                                             gCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGCTTTTTTGCg                                           >  1:4531905/1‑93 (MQ=255)

                                               tATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCgcc                                         >  1:3742822/1‑91 (MQ=255)

                                                aTTTTCACTTGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGccc                                        <  2:605129/93‑1 (MQ=255)

                                                aTTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGccc                                        <  1:752518/93‑1 (MQ=255)

                                                aTTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGccc                                        >  1:3387157/1‑93 (MQ=255)

                                                  tttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAg                                      <  1:2290202/93‑1 (MQ=255)

                                                    tCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGc                                    >  1:584242/1‑93 (MQ=255)

                                                    tCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGc                                    <  1:2117480/93‑1 (MQ=255)

                                                       cTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCttt                                 >  2:4537491/1‑93 (MQ=255)

                                                         ggCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGc                               <  2:2625070/93‑1 (MQ=255)

                                                         ggCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGc                               >  2:2630758/1‑93 (MQ=255)

                                                          gCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTTGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGcc                              >  2:2476565/1‑93 (MQ=255)

                                                            tGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCgg                            >  1:303473/1‑93 (MQ=255)

                                                             gttgcTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCGGCCGTCGCTGCTGCCTGCGCGATTGATTTTTTGCGCCCACGCTTTGCCGga                           <  2:2991245/89‑1 (MQ=255)

                                                             gATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGga                           <  1:2827218/93‑1 (MQ=255)

                                                             gATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGga                           >  1:852707/1‑93 (MQ=255)

                                                              aTACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGgag                          >  1:1769867/1‑93 (MQ=255)

                                                               tACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGgaga                         >  2:3760237/1‑93 (MQ=255)

                                                                aCTGCCTTTGCCTACGCCTGCAATCGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGac                        <  1:1532816/93‑1 (MQ=255)

                                                                 cTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGaca                       <  1:4318508/93‑1 (MQ=255)

                                                                  tGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGacac                      <  2:1440143/93‑1 (MQ=255)

                                                                  tGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGacac                      >  1:1669270/1‑93 (MQ=255)

                                                                  tGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGacac                      <  1:444226/93‑1 (MQ=255)

                                                                   gCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACAcc                     >  2:3067389/1‑93 (MQ=255)

                                                                   gCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACAcc                     <  1:1155528/93‑1 (MQ=255)

                                                                    ccTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCt                    >  1:637984/1‑93 (MQ=255)

                                                                    ccTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCt                    >  1:2778649/1‑93 (MQ=255)

                                                                     cTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCtt                   >  1:193693/1‑93 (MQ=255)

                                                                     cTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCtt                   >  1:2441814/1‑93 (MQ=255)

                                                                      tttGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTa                  >  2:3459686/1‑93 (MQ=255)

                                                                         cccTCCGCCTGCAATAGTCAGGGGCTGGCACCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACt               <  1:2653/92‑1 (MQ=255)

                                                                          ccTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTg              <  2:492983/93‑1 (MQ=255)

                                                                            tACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt            >  1:1934706/1‑93 (MQ=255)

                                                                            tACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt            >  2:3659145/1‑93 (MQ=255)

                                                                            tACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt            >  1:4634045/1‑93 (MQ=255)

                                                                            tACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt            >  1:2717768/1‑93 (MQ=255)

                                                                             aCGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTa           <  2:3138954/93‑1 (MQ=255)

                                                                              cGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc          <  1:502220/93‑1 (MQ=255)

                                                                              cGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc          <  1:1751967/93‑1 (MQ=255)

                                                                               gCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACt         >  2:2148662/1‑93 (MQ=255)

                                                                                ccTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTg        >  2:1575668/1‑93 (MQ=255)

                                                                                 cTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc       >  1:309978/1‑93 (MQ=255)

                                                                                 cTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc       >  2:3224942/1‑93 (MQ=255)

                                                                                  tGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgcg      >  2:80437/1‑93 (MQ=255)

                                                                                   gCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgcgc     >  1:97960/1‑93 (MQ=255)

                                                                                     aaTAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAg   <  1:4239995/93‑1 (MQ=255)

                                                                                      aTAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAgg  >  2:4559693/1‑93 (MQ=255)

                                                                                      aTAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAgg  >  2:698694/1‑93 (MQ=255)

CACTGCACAAATGCTTAACGGTCATCAAAGTTGCCACGGCGGGCAGCTATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGG  >  NC_000913/1458614‑1458792

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