BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,533	A→G	intergenic (+88/‑59)	lysY → / → lysZ	tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,533	0	A	G	99.6%	819.3 / ‑5.3	269	intergenic (+88/‑59)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (147/120); total (149/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 288 total aligned reads displayed.

ACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGA  >  NC_000913/780442‑780622

 ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGcgcg                                                                                         >  1:4350328/1‑93 (MQ=255)

    gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGTGCGAcc                                                                                      >  2:1600497/1‑93 (MQ=255)

    gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGAcc                                                                                      >  1:6947813/1‑93 (MQ=255)

     ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGAccc                                                                                     >  1:1089982/1‑93 (MQ=255)

      tttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGAcccc                                                                                    >  2:2139181/1‑92 (MQ=255)

      tttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACcca                                                                                    >  2:2705295/1‑93 (MQ=255)

       ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACccac                                                                                   >  2:2462292/1‑93 (MQ=255)

       ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACccac                                                                                   >  1:3260742/1‑93 (MQ=255)

         aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACccacca                                                                                 >  1:7847591/1‑93 (MQ=255)

         aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACccacca                                                                                 >  2:1160351/1‑93 (MQ=255)

         aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACccacca                                                                                 >  2:7350986/1‑93 (MQ=255)

         aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACccacca                                                                                 >  2:4107638/1‑93 (MQ=255)

          aCATCAACCTAAGATGTTAAGCGTGAAGGATAACGTTGCACCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCaa                                                                                <  2:2370692/93‑1 (MQ=14)

          aCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCACCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCaa                                                                                <  2:3042915/93‑1 (MQ=255)

          aCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCaa                                                                                >  1:6273165/1‑93 (MQ=255)

             tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGACATCCTGCGCGACCCACCAATCt                                                                             >  2:344719/1‑93 (MQ=255)

              cAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCtt                                                                            >  2:537413/1‑93 (MQ=255)

               aaaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTa                                                                           >  2:6201093/1‑93 (MQ=255)

                aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTaa                                                                          >  1:1491112/1‑93 (MQ=255)

                 agtCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTaaa                                                                         <  2:2823716/91‑1 (MQ=255)

                 aCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTaaa                                                                         <  1:2874256/93‑1 (MQ=255)

                 aCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTaaa                                                                         <  1:601207/93‑1 (MQ=255)

                 aCTCAGATGTTAAGCGTGAAGGATAACGTGGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTaaa                                                                         <  1:6156420/93‑1 (MQ=255)

                  cTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAg                                                                        >  1:7442321/1‑93 (MQ=255)

                  cTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAg                                                                        >  1:2920780/1‑93 (MQ=255)

                  cTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAg                                                                        >  1:1886344/1‑93 (MQ=255)

                  cTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAg                                                                        >  2:7548517/1‑93 (MQ=255)

                    cAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGAt                                                                      >  1:4465032/1‑93 (MQ=255)

                     aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGAtt                                                                     <  1:2399486/93‑1 (MQ=255)

                     aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGAtt                                                                     <  1:1053658/93‑1 (MQ=255)

                      gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTg                                                                    >  1:5196735/1‑93 (MQ=255)

                      gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTg                                                                    >  1:1060636/1‑93 (MQ=255)

                       aTGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTgg                                                                   >  2:5109375/1‑93 (MQ=255)

                        tGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGc                                                                  <  2:3380529/93‑1 (MQ=255)

                         gTTAAGCGTGAAGGATAACGCTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTACAGATTGGcc                                                                 >  2:7613019/1‑93 (MQ=255)

                           tAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGcccc                                                               >  1:5707314/1‑93 (MQ=255)

                           tAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAATTTGGcccg                                                               >  2:3750595/1‑92 (MQ=17)

                            aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCg                                                              <  1:3749012/93‑1 (MQ=255)

                             aGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGa                                                             >  2:7873901/1‑93 (MQ=255)

                             aGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGa                                                             >  1:7064844/1‑93 (MQ=255)

                              gCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCCCGACCCACCAATCTTAAAGATTGGCCCCGAg                                                            >  1:1757147/1‑93 (MQ=255)

                                gTGAGGGATACCGTGGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTa                                                          <  2:4224449/93‑1 (MQ=21)

                                gTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATATTAAAGATTGGCCCCGAGTa                                                          >  2:6029127/1‑93 (MQ=255)

                                 tGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTaa                                                         >  1:2124359/1‑93 (MQ=255)

                                   aaGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaa                                                       >  2:5742694/1‑93 (MQ=255)

                                     ggATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                                     >  1:3876780/1‑93 (MQ=255)

                                     ggATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                                     >  1:5520621/1‑93 (MQ=255)

                                      gATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGCGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                                    >  2:2221125/1‑93 (MQ=255)

                                      gATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                                    >  1:1657577/1‑93 (MQ=255)

                                      gATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                                    >  2:1627577/1‑93 (MQ=255)

                                        tAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                                  >  2:6814909/1‑93 (MQ=255)

                                        tAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                                  >  2:5667088/1‑93 (MQ=255)

                                         aaCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                                 >  1:5731268/1‑93 (MQ=255)

                                           cGTTGCACCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                               >  2:5468700/1‑93 (MQ=255)

                                            gTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                              >  2:5157073/1‑93 (MQ=255)

                                             ttGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAgg                                             >  1:1456102/1‑93 (MQ=255)

                                              tGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                                            <  2:4139596/93‑1 (MQ=255)

                                               gcgcCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                                           >  1:1602051/1‑93 (MQ=255)

                                               gcgcCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                                           >  1:5312699/1‑93 (MQ=255)

                                                 gcCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAac                                         <  1:2560889/93‑1 (MQ=255)

                                                   cAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAacac                                       >  1:7300004/1‑93 (MQ=255)

                                                    aGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAcc                                      >  2:4761937/1‑93 (MQ=255)

                                                     gCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAACATCTTTCAGGTAACAccc                                     >  2:4932970/1‑93 (MQ=255)

                                                     gCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAccc                                     >  2:7778631/1‑93 (MQ=255)

                                                      cAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCg                                    <  1:6586903/93‑1 (MQ=255)

                                                        aCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGGCCCGAGTAAAAATCTTTCAGGTAACACCCGTa                                  >  2:383456/1‑93 (MQ=255)

                                                         cGGCCCGTAGGGCGAGAGAAGCGAGTCATCCTGCGCGCCCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                 <  2:4639432/93‑1 (MQ=255)

                                                           gCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                               >  1:263108/1‑93 (MQ=255)

                                                            cccGTAGGTCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                              <  2:6922229/93‑1 (MQ=255)

                                                             ccGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                             <  1:4829137/93‑1 (MQ=255)

                                                             ccGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                             >  2:4704878/1‑93 (MQ=255)

                                                              cGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTc                            >  1:1212969/1‑93 (MQ=255)

                                                                  ggtcgagcgaAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa                        <  2:6478936/90‑1 (MQ=255)

                                                                  gggCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAACATCTTTCAGGTAACACCCGTATGGGTCGTTa                        >  1:7575777/1‑93 (MQ=255)

                                                                  gggCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa                        <  1:3494192/93‑1 (MQ=255)

                                                                   ggCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAg                       >  2:1976619/1‑93 (MQ=255)

                                                                    gcgagcgaAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGc                      >  1:7756322/1‑93 (MQ=255)

                                                                    gcgagcgaAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGc                      <  2:5996877/93‑1 (MQ=255)

                                                                    gcgagcgaAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGc                      <  1:4933613/93‑1 (MQ=255)

                                                                    gcgagcgaAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGc                      <  2:5410595/93‑1 (MQ=255)

                                                                      gagcgaAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc                    >  1:358281/1‑93 (MQ=255)

                                                                       agcgaagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa                   <  2:6506151/93‑1 (MQ=255)

                                                                         cgaagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt                 <  2:1031588/93‑1 (MQ=255)

                                                                          gaagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGtt                <  1:7149723/93‑1 (MQ=255)

                                                                          gaagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGtt                <  1:6363964/93‑1 (MQ=255)

                                                                           aagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTg               >  1:5014188/1‑93 (MQ=255)

                                                                           aagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTg               >  2:3070915/1‑93 (MQ=255)

                                                                           aagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTg               <  1:7682609/93‑1 (MQ=255)

                                                                             gcgaGTCATCCTGCGCGCCCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGt             <  2:7150864/93‑1 (MQ=255)

                                                                             gcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGt             >  2:1283706/1‑93 (MQ=255)

                                                                              cgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa            >  1:6833312/1‑93 (MQ=255)

                                                                              cgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa            <  2:6704337/93‑1 (MQ=255)

                                                                               gaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag           >  1:4421469/1‑93 (MQ=255)

                                                                                aGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTaga          >  2:1392238/1‑93 (MQ=255)

                                                                                aGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTaga          <  2:3356596/93‑1 (MQ=255)

                                                                                aGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTaga          <  1:3079138/93‑1 (MQ=255)

                                                                                    aTCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAg      <  2:685508/93‑1 (MQ=255)

                                                                                     tCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt     >  1:2327726/1‑93 (MQ=255)

                                                                                      ccTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGtt    <  1:911096/93‑1 (MQ=255)

                                                                                      ccTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGtt    <  1:2034231/93‑1 (MQ=255)

ACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGA  >  NC_000913/780442‑780622

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