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	98.4%	378.8 / ‑2.6	123	intergenic (+88/‑59)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (71/50); total (73/50)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 131 total aligned reads displayed.

CCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGA  >  NC_000913/780443‑780622

ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGTCGAGCGAAGCGAGTCATCCTGcgcg                                                                                         >  2:4153895/1‑93 (MQ=255)

    ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCGACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGAccc                                                                                     <  1:2423673/93‑1 (MQ=255)

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

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

       tAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACccacc                                                                                  >  1:224339/1‑93 (MQ=255)

          cATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAAt                                                                               >  2:3736329/1‑93 (MQ=255)

           aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATc                                                                              >  2:4304326/1‑93 (MQ=255)

            tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCt                                                                             <  1:3150852/93‑1 (MQ=255)

            tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCt                                                                             >  2:2040924/1‑93 (MQ=255)

            tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCAGCCTGCGCGACCCACCAATCt                                                                             >  1:916887/1‑93 (MQ=255)

            tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGAGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATcc                                                                             >  1:147040/1‑92 (MQ=255)

             cAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCtt                                                                            <  1:3185597/93‑1 (MQ=255)

             cAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCtt                                                                            <  1:3629183/93‑1 (MQ=255)

              aaaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTa                                                                           <  1:1652004/93‑1 (MQ=255)

                acacAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTaaa                                                                         <  2:577192/90‑1 (MQ=255)

                aCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTaaa                                                                         >  2:4073472/1‑93 (MQ=255)

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

                  tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGa                                                                       <  2:1372464/93‑1 (MQ=255)

                  tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGa                                                                       >  2:2525915/1‑93 (MQ=255)

                  tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGa                                                                       >  2:4659271/1‑93 (MQ=255)

                   cAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGAt                                                                      <  1:703319/93‑1 (MQ=255)

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

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

                    aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGAtt                                                                     >  2:2115227/1‑93 (MQ=255)

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

                     gATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCTACCCACCATTCTTAAAGATTg                                                                    <  2:565251/93‑1 (MQ=14)

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

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

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

                       tGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGc                                                                  >  1:592371/1‑93 (MQ=255)

                         ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGccc                                                                >  2:3904425/1‑93 (MQ=255)

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

                             gCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAg                                                            <  1:1292832/93‑1 (MQ=255)

                             ccGTGAAGTATACCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTTGCCCCGAg                                                            <  2:413039/92‑1 (MQ=18)

                                tGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGtca                                                         >  1:866923/1‑91 (MQ=255)

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

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

                                 gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                                        >  1:4572162/1‑93 (MQ=255)

                                  aaGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaa                                                       <  2:3071020/93‑1 (MQ=255)

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

                                   aGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaaa                                                      <  1:310478/93‑1 (MQ=255)

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

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

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

                                    gaaaaaCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                                     <  2:4190021/89‑1 (MQ=255)

                                     gATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                                    <  2:62392/93‑1 (MQ=255)

                                      aTAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTATAGATTGGCCCCGAGTAAAAATCt                                                   >  2:2582911/1‑93 (MQ=255)

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

                                         aCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                <  2:144714/93‑1 (MQ=255)

                                         aCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                >  2:57199/1‑93 (MQ=255)

                                         aCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                >  1:1076651/1‑93 (MQ=255)

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

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

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

                                           gTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                              <  1:3183589/93‑1 (MQ=255)

                                               cgcCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                          <  1:2548456/93‑1 (MQ=255)

                                                gcCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAac                                         >  1:4255767/1‑93 (MQ=255)

                                                 ccAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAaca                                        <  1:2634626/93‑1 (MQ=255)

                                                   aGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTACAGATTGGCCCCCAGTAAAAATCTTTCAGCTAACAcc                                      >  1:1412975/1‑93 (MQ=21)

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

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

                                                      aaCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACCCCCGt                                   >  1:2263600/1‑93 (MQ=255)

                                                      aaCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGt                                   >  1:4605051/1‑93 (MQ=255)

                                                        cGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCTACCCCCCAATCTTAAAGATTGGCCCCGAGTAAACATCTTTCAGGTAAAACCCgcag                                 >  2:4485301/1‑90 (MQ=25)

                                                        cGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCTACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                 >  2:4485303/1‑93 (MQ=255)

                                                        cGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                 <  1:3158128/93‑1 (MQ=255)

                                                         ggCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                                <  1:2402619/93‑1 (MQ=255)

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

                                                              gTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCg                           >  1:2005971/1‑93 (MQ=255)

                                                              gTAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCg                           >  1:542216/1‑93 (MQ=255)

                                                               tAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt                          >  2:401308/1‑93 (MQ=255)

                                                               tAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt                          >  1:1142345/1‑93 (MQ=255)

                                                                aGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGtt                         >  2:3795110/1‑93 (MQ=255)

                                                                 gggCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAACAATCTTTCAGGTAACACCCGTATGGGTCGTTa                        >  2:3255830/1‑93 (MQ=255)

                                                                 gggCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa                        >  2:3255829/1‑93 (MQ=255)

                                                                 gggCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa                        >  2:81781/1‑93 (MQ=255)

                                                                  ggCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAg                       <  1:1630825/93‑1 (MQ=255)

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

                                                                      agggaagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa                   <  1:3378158/90‑1 (MQ=255)

                                                                       gcgaagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACTCCCGTATGGGTCGTTAGCTCAg                  <  2:604411/93‑1 (MQ=255)

                                                                       gcgaagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAg                  >  1:3841152/1‑93 (MQ=255)

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

                                                                         gaagcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGtt                >  2:534527/1‑93 (MQ=255)

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

                                                                           agcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg              <  1:2580270/93‑1 (MQ=255)

                                                                            gcgaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGt             <  2:2459784/93‑1 (MQ=255)

                                                                            acgaGTCATCCAAAGCGACCCACAAATCATGACGATTGGCCCCGAGTCAATAACTTTCAGATAACACCCGTATGGGTCGTTTGCTCAGATGGt             <  1:4222492/92‑1 (MQ=255)

                                                                              gaGTCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag           <  2:1057957/93‑1 (MQ=255)

                                                                              gaGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag           >  2:4010593/1‑93 (MQ=255)

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

                                                                                 tCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGc        <  2:1351190/93‑1 (MQ=255)

                                                                                 tCATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGc        <  2:1351189/93‑1 (MQ=255)

                                                                                  cATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa       >  2:3126541/1‑93 (MQ=255)

                                                                                  cATCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa       >  1:4729127/1‑93 (MQ=255)

                                                                                    tCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt     <  2:4727673/93‑1 (MQ=255)

                                                                                    tCCTGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt     <  1:170419/93‑1 (MQ=255)

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

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

                                                                                       tGCGCGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTTAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa  >  1:4206932/1‑93 (MQ=255)

                                                                                       tGCGCGACCCACCAATCTTAAAGATTGGCCACGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa  >  1:2365516/1‑93 (MQ=255)

CCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGA  >  NC_000913/780443‑780622

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