BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	440,252	A→G	L25P (CTG→CCG)	ispA ←	geranyltranstransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	440,252	0	A	G	99.6%	852.3 / ‑4.6	247	L25P (CTG→CCG)	ispA	geranyltranstransferase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (114/132); total (115/132)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

Only 100 of 252 total aligned reads displayed.

GTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGT  >  NC_000913/440162‑440340

gTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGt                                                                                        >  2:3204562/1‑93 (MQ=255)

  ggCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTgg                                                                                      <  2:5888289/93‑1 (MQ=255)

   gCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTggg                                                                                     <  1:5921681/93‑1 (MQ=255)

   gCATAAACCAGGAAAGGTCGCAGGAGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTggg                                                                                     >  1:2611935/1‑93 (MQ=255)

      tAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCg                                                                                  >  1:5784338/1‑93 (MQ=255)

       aaaCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGa                                                                                 >  1:6455985/1‑93 (MQ=255)

       aaaCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGa                                                                                 <  1:855436/93‑1 (MQ=255)

        aaCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGAt                                                                                <  1:4612186/93‑1 (MQ=255)

        aaCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGCGTTCTGAAAGGGCGGTGGGGCGAt                                                                                >  2:737448/1‑93 (MQ=255)

         aCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATa                                                                               >  2:4691883/1‑93 (MQ=255)

         aCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATa                                                                               >  2:369255/1‑93 (MQ=255)

         aCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATa                                                                               >  1:911928/1‑93 (MQ=255)

           cAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATaaa                                                                             >  1:2287885/1‑93 (MQ=255)

           cAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGACTGTTCTGAAAGGGCGGCGGGGCGATaaa                                                                             >  1:336795/1‑93 (MQ=255)

            aGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATaaaa                                                                            >  2:6315900/1‑93 (MQ=255)

            aGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATaaaa                                                                            <  2:5845214/93‑1 (MQ=255)

             ggAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAAc                                                                           <  1:5306798/93‑1 (MQ=255)

             ggAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAAc                                                                           >  2:3756165/1‑93 (MQ=255)

             ggAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAAc                                                                           <  2:1262351/93‑1 (MQ=255)

             ggAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACAGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAAc                                                                           <  1:4423858/93‑1 (MQ=255)

              gAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACg                                                                          >  2:6091760/1‑93 (MQ=255)

               aaaGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACgg                                                                         >  2:6387977/1‑93 (MQ=255)

                aaGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGc                                                                        >  2:5050647/1‑93 (MQ=255)

                aaGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGc                                                                        <  2:4963433/93‑1 (MQ=255)

                   gTCGCAGGCGCTTACCACCTAATAATGCGCCATACCGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCa                                                                     >  2:164787/1‑93 (MQ=255)

                    tCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAg                                                                    <  1:2162710/93‑1 (MQ=255)

                    tCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAg                                                                    >  2:1756207/1‑93 (MQ=255)

                    tCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAg                                                                    <  2:4452260/93‑1 (MQ=255)

                    tCGCAGGCGCTTACCACCTAAAAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAg                                                                    >  1:4092501/1‑93 (MQ=255)

                      gCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAgcg                                                                  <  1:4242175/93‑1 (MQ=255)

                        aGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGcc                                                                >  1:804806/1‑93 (MQ=255)

                        aGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGcc                                                                <  1:1030703/93‑1 (MQ=255)

                         ggCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCt                                                               >  2:6205711/1‑93 (MQ=255)

                         ggCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCt                                                               <  1:3429813/93‑1 (MQ=255)

                         ggCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCt                                                               <  2:3425346/93‑1 (MQ=255)

                         ggCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCt                                                               <  1:3858345/93‑1 (MQ=255)

                         ggCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCt                                                               <  1:3952837/93‑1 (MQ=255)

                          gcgcTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGTGATAAAACGGCTCAGCGcccg                                                              >  2:1984316/1‑91 (MQ=255)

                           cgcTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTgg                                                             <  2:871353/93‑1 (MQ=255)

                             cTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGtt                                                           <  2:1007050/93‑1 (MQ=255)

                              ttACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTg                                                          <  2:4864619/93‑1 (MQ=255)

                              ttACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTg                                                          <  2:421632/93‑1 (MQ=255)

                              ttACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTg                                                          <  2:5232535/93‑1 (MQ=255)

                                 ccaccTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGcc                                                       <  1:6033217/93‑1 (MQ=255)

                                   accTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTg                                                     >  2:5730184/1‑93 (MQ=255)

                                    ccTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAGCGGCTCAGCGCCTGGTTGGCCTGc                                                    <  1:6153540/93‑1 (MQ=255)

                                    ccTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGc                                                    <  2:3672624/93‑1 (MQ=255)

                                    ccTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGc                                                    >  2:4987197/1‑93 (MQ=255)

                                     cTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCt                                                   >  1:6048467/1‑93 (MQ=255)

                                     cTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCGGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCt                                                   <  2:2936916/93‑1 (MQ=255)

                                      taataaTGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCtt                                                  <  2:1937710/93‑1 (MQ=255)

                                      taataaTGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCtt                                                  >  1:1675946/1‑93 (MQ=255)

                                       aataatGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTa                                                 >  1:5019861/1‑93 (MQ=255)

                                          aatGCGCCATACTGCAAGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACg                                              >  2:2231290/1‑93 (MQ=255)

                                           atGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGc                                             <  1:691035/93‑1 (MQ=255)

                                               gcCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGAAGGc                                         >  2:4771570/1‑93 (MQ=255)

                                                ccATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCt                                        >  1:2904101/1‑93 (MQ=255)

                                                ccATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCt                                        <  1:3719163/93‑1 (MQ=255)

                                                ccATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCt                                        <  2:5291947/93‑1 (MQ=255)

                                                 cATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCtt                                       <  1:6073039/93‑1 (MQ=255)

                                                 cATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCtt                                       <  1:1970808/93‑1 (MQ=255)

                                                 cATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGATGGGGCGATAAAACGGCTCAGCGCCAGGTTGGCCTGCTTAACGCAGGCtt                                       <  1:6472538/93‑1 (MQ=255)

                                                  aTACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTc                                      <  1:5591442/93‑1 (MQ=255)

                                                  aTACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTc                                      >  2:2827357/1‑93 (MQ=255)

                                                  aTACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTc                                      <  2:2166436/93‑1 (MQ=255)

                                                  aTACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTc                                      >  1:486670/1‑93 (MQ=255)

                                                   tACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCg                                     <  1:3730274/93‑1 (MQ=255)

                                                   tACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCg                                     <  1:3596757/93‑1 (MQ=255)

                                                    aCTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGa                                    <  2:3621902/93‑1 (MQ=255)

                                                    aCTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGa                                    >  2:824769/1‑93 (MQ=255)

                                                    aCTGCATGGTTCCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGa                                    <  2:3883590/93‑1 (MQ=255)

                                                     cTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAg                                   >  2:3936656/1‑93 (MQ=255)

                                                        cATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtg                                <  2:4587525/93‑1 (MQ=255)

                                                         aTGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgc                               <  2:245299/93‑1 (MQ=255)

                                                          tGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgct                              <  2:5600331/93‑1 (MQ=255)

                                                          tGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgct                              >  2:549721/1‑93 (MQ=255)

                                                          tGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgct                              >  1:3025137/1‑93 (MQ=255)

                                                          tGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgct                              <  1:5096384/93‑1 (MQ=255)

                                                          tGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgct                              >  1:1356728/1‑93 (MQ=255)

                                                           ggTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgctg                             >  2:3992598/1‑93 (MQ=255)

                                                           ggTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgctg                             <  2:4050395/93‑1 (MQ=255)

                                                                 gACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGaaa                       <  1:1008397/93‑1 (MQ=255)

                                                                  aCCACGGGTGTGTTCTGAATGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGCAAAg                      >  1:6312853/1‑93 (MQ=255)

                                                                  aCCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAg                      >  1:5566049/1‑93 (MQ=255)

                                                                   ccACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGt                     >  1:335208/1‑93 (MQ=255)

                                                                     aCGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTcc                   <  1:4456763/93‑1 (MQ=255)

                                                                     aCGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTcc                   >  1:1803669/1‑93 (MQ=255)

                                                                     aCGGCAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTcc                   >  1:1615421/1‑93 (MQ=255)

                                                                      cGGGAGTGTGCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                  <  1:3185871/93‑1 (MQ=255)

                                                                        ggAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCAtt                >  1:3679456/1‑93 (MQ=255)

                                                                          aGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTAc              >  1:232991/1‑93 (MQ=255)

                                                                            tgtTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTc            <  2:2347533/93‑1 (MQ=255)

                                                                            tgtTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTc            <  1:1227720/93‑1 (MQ=255)

                                                                               tCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCAtt         >  1:5397653/1‑93 (MQ=255)

                                                                                cTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTg        <  2:6189287/93‑1 (MQ=255)

                                                                                 tGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGt       >  1:4043022/1‑93 (MQ=255)

                                                                                 tGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGt       >  1:2516854/1‑93 (MQ=255)

                                                                                  gAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTc      <  1:1745774/93‑1 (MQ=255)

                                                                                   aaaGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTcc     >  2:5157261/1‑93 (MQ=255)

                                                                                     aGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCgg   <  2:5430518/93‑1 (MQ=255)

GTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGT  >  NC_000913/440162‑440340

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