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	98.9%	866.3 / ‑2.3	260	L25P (CTG→CCG)	ispA	geranyltranstransferase
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (117/140); total (117/143)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

GTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTG  >  NC_000913/440162‑440341

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

gTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGt                                                                                         <  2:4592694/93‑1 (MQ=255)

gTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGt                                                                                         <  1:5340515/93‑1 (MQ=255)

 tGGCATAAACCAGGGAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTg                                                                                        <  1:1626517/93‑1 (MQ=255)

 tGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTg                                                                                        >  1:3268242/1‑93 (MQ=255)

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

  ggCATAAACCAGGAAAGGTCGCACGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTgg                                                                                       <  1:2674120/93‑1 (MQ=255)

   gCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTggg                                                                                      >  1:1955464/1‑93 (MQ=255)

    cATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTgggg                                                                                     <  2:3851364/93‑1 (MQ=255)

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

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

        aaCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGAt                                                                                 >  1:70175/1‑93 (MQ=255)

          ccAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATaa                                                                               <  2:5104746/93‑1 (MQ=255)

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

              gAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACg                                                                           <  2:3487549/93‑1 (MQ=255)

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

                 aGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCt                                                                        >  2:1796295/1‑93 (MQ=255)

                 aGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCt                                                                        >  2:4279922/1‑93 (MQ=255)

                  ggTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTc                                                                       <  1:7249345/93‑1 (MQ=255)

                  ggTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTc                                                                       <  1:484281/93‑1 (MQ=255)

                  ggTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTc                                                                       >  2:284062/1‑93 (MQ=255)

                   gTCGCAGGCGCTTACCACCTAATATTGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCa                                                                      <  2:6310085/93‑1 (MQ=255)

                   gTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCa                                                                      >  1:228762/1‑93 (MQ=255)

                     cGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAgc                                                                    <  1:417222/93‑1 (MQ=255)

                     cGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAgc                                                                    <  2:4164083/93‑1 (MQ=255)

                     cGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAgc                                                                    <  1:2040490/93‑1 (MQ=255)

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

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

                        aGGCGCTTACCACCGAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGcc                                                                 >  2:5193324/1‑93 (MQ=255)

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

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

                          gcgcTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTg                                                               <  2:6078713/93‑1 (MQ=255)

                              ttACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTg                                                           >  2:6286546/1‑93 (MQ=255)

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

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

                               tACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTgg                                                          <  2:2258542/93‑1 (MQ=255)

                                accaccTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGc                                                         <  1:48684/93‑1 (MQ=255)

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

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

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

                                     cTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCt                                                    <  2:3424776/93‑1 (MQ=255)

                                     cTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCt                                                    <  2:1364845/93‑1 (MQ=255)

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

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

                                       aataatGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTa                                                  <  1:1519604/93‑1 (MQ=255)

                                        ataatGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTaa                                                 >  2:4968651/1‑93 (MQ=255)

                                        ataatGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTaa                                                 <  2:7193903/93‑1 (MQ=255)

                                         taatGCGCCATACTGCATGGTTTCGACGACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAAc                                                >  2:6673658/1‑93 (MQ=255)

                                         taatGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAAc                                                <  1:4969281/93‑1 (MQ=255)

                                           atGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGc                                              >  1:6263760/1‑93 (MQ=255)

                                            tGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCa                                             <  2:2769600/93‑1 (MQ=255)

                                            tGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCa                                             <  2:485399/93‑1 (MQ=255)

                                               gcCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGc                                          <  2:4800100/93‑1 (MQ=255)

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

                                                 cATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCtt                                        >  1:6879216/1‑93 (MQ=255)

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

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

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

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

                                                   tACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCg                                      >  1:3518443/1‑93 (MQ=255)

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

                                                      tGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGt                                   >  1:1823621/1‑93 (MQ=255)

                                                       gCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGtt                                  <  2:5348139/93‑1 (MQ=255)

                                                       gCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGtt                                  >  2:5730465/1‑93 (MQ=255)

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

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

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

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

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

                                                            gTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgctgc                             <  1:3101281/93‑1 (MQ=255)

                                                            gTTTCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgctgc                             >  1:7171141/1‑93 (MQ=255)

                                                             tttCGACCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCg                            <  1:3862859/93‑1 (MQ=255)

                                                                  aCCACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAg                       <  1:4526273/93‑1 (MQ=255)

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

                                                                    cACGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTc                     >  2:6795095/1‑93 (MQ=255)

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

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

                                                                      cGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                   <  2:2371245/93‑1 (MQ=255)

                                                                      cGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                   <  2:1003761/93‑1 (MQ=255)

                                                                      cGGGAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                   <  1:682381/93‑1 (MQ=255)

                                                                       gggAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCAt                  >  2:1162648/1‑93 (MQ=255)

                                                                        ggAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCAtt                 <  2:5421631/93‑1 (MQ=255)

                                                                         gAGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTa                <  1:4587590/93‑1 (MQ=255)

                                                                          aGTGTTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTAc               <  2:6204746/93‑1 (MQ=255)

                                                                           gtgtTCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACt              >  2:7317396/1‑93 (MQ=255)

                                                                           gtgtTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACt              <  2:5608586/93‑1 (MQ=255)

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

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

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

                                                                               tCTGAAAGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCAtt          <  1:5759499/93‑1 (MQ=255)

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

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

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

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

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

                                                                                   aaaGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTcc      <  1:394195/93‑1 (MQ=255)

                                                                                     aGGGCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCgg    >  2:1917449/1‑93 (MQ=255)

                                                                                      gggCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGgt   >  2:2514588/1‑93 (MQ=255)

                                                                                       ggCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGgtg  <  2:1122183/93‑1 (MQ=255)

                                                                                       ggCGGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGgtg  >  2:3211413/1‑93 (MQ=255)

GTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTG  >  NC_000913/440162‑440341

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