BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,554,713	T→C	V22A (GTA→GCA)	osmC →	OsmC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,554,713	0	T	C	96.7%	270.5 / ‑0.6	92	V22A (GTA→GCA)	osmC	OsmC
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (46/43); total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTTAAGCCCACAGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGTATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTGAAGAACTGA  >  NC_000913/1554626‑1554802

tttAAGCCCACAGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCa                                                                                      >  1:249799/1‑93 (MQ=255)

 ttAAGCCCACAGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCaa                                                                                     <  1:2403446/93‑2 (MQ=255)

 ttAAGCCCACAGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCAc                                                                                     <  2:2904361/93‑1 (MQ=255)

    aGCCCACAGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCga                                                                                  >  1:2169356/1‑93 (MQ=255)

      cccACAGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCgaga                                                                                >  2:3339247/1‑93 (MQ=255)

        cacaGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGt                                                                              <  1:350567/93‑1 (MQ=255)

        cacaGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGt                                                                              >  1:2865711/1‑93 (MQ=255)

         acaGGAGAGCGACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTg                                                                             <  2:2840229/93‑1 (MQ=255)

         acaGGAGAACAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTg                                                                             <  1:3096814/93‑1 (MQ=255)

          caGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTgg                                                                            >  1:137347/1‑93 (MQ=255)

             gagaGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGt                                                                         <  2:1856637/93‑1 (MQ=255)

               gagCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGc                                                                       <  1:3415450/93‑1 (MQ=255)

               gagCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGc                                                                       <  1:2163063/93‑1 (MQ=255)

               gagCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGc                                                                       <  2:348301/93‑1 (MQ=255)

                agCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCt                                                                      <  2:2091109/93‑1 (MQ=255)

                   aacaaTGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGaa                                                                   >  1:287742/1‑93 (MQ=255)

                   aacaaTGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGaa                                                                   >  1:1523354/1‑93 (MQ=255)

                      aaTGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAAcaa                                                                >  2:1588319/1‑91 (MQ=255)

                      aaTGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCa                                                                >  2:1807398/1‑93 (MQ=255)

                       aTGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCaa                                                               >  1:4423261/1‑93 (MQ=255)

                       aTGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCaa                                                               >  2:2639230/1‑93 (MQ=255)

                        tGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAAc                                                              <  1:2004049/93‑1 (MQ=255)

                         gACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACa                                                             <  1:3013268/93‑1 (MQ=255)

                           cAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGc                                                           >  2:2902851/1‑93 (MQ=255)

                             aTCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCg                                                         >  1:1904763/1‑93 (MQ=255)

                              tCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGt                                                        >  2:3291792/1‑93 (MQ=255)

                              tCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGCAGTGCTGAACCAACAGCCGt                                                        >  2:2428758/1‑93 (MQ=255)

                                cATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTAt                                                      <  2:3808872/93‑1 (MQ=255)

                                cATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTAt                                                      >  2:4357025/1‑93 (MQ=255)

                                cATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTAt                                                      <  1:3063765/93‑1 (MQ=255)

                                 aTAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGTATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATg                                                     >  2:2218222/1‑93 (MQ=255)

                                 aTAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATg                                                     <  1:1346671/93‑1 (MQ=255)

                                 aTAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATg                                                     <  2:2597552/93‑1 (MQ=255)

                                 aTAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGCGGCGTGCTGAACCAACAGCCGTATg                                                     >  2:3295371/1‑93 (MQ=255)

                                   aaGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGa                                                   <  2:4131203/93‑1 (MQ=255)

                                    aGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGAt                                                  <  2:3983185/93‑1 (MQ=255)

                                     gAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGTATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGAtt                                                 >  2:4606734/1‑93 (MQ=255)

                                     gAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGAtt                                                 <  1:3355959/93‑1 (MQ=255)

                                     gAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGAtt                                                 >  2:1552594/1‑93 (MQ=255)

                                     gAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGAtt                                                 >  1:1139471/1‑93 (MQ=255)

                                       aaGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTa                                               <  1:1242635/93‑1 (MQ=255)

                                        aGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTaa                                              <  1:3072651/93‑1 (MQ=255)

                                        aGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTaa                                              <  1:48356/93‑1 (MQ=255)

                                          gTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAaca                                            >  2:1849409/1‑93 (MQ=255)

                                          gTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAaca                                            <  1:3553776/93‑1 (MQ=255)

                                           tCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGTATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAacac                                           <  1:607113/93‑1 (MQ=255)

                                           tCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAacac                                           >  1:3236287/1‑93 (MQ=255)

                                           tCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAacac                                           <  2:41822/93‑1 (MQ=255)

                                              ggCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACAcgcg                                        >  1:763062/1‑93 (MQ=255)

                                                cacaCTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGtt                                      >  2:3477343/1‑93 (MQ=255)

                                                cacaCTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGtt                                      <  1:2034053/93‑1 (MQ=255)

                                                cacaCGGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTAGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGtt                                      <  1:3396100/93‑1 (MQ=255)

                                                 acacTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGttt                                     >  2:209639/1‑93 (MQ=255)

                                                 acacTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGttt                                     >  2:1101533/1‑93 (MQ=255)

                                                 acacTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGttt                                     >  2:1169080/1‑93 (MQ=255)

                                                 acacTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGttt                                     >  1:2074185/1‑93 (MQ=255)

                                                  cacTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGtttt                                    >  1:4293250/1‑93 (MQ=255)

                                                   acTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTg                                   >  1:1133932/1‑93 (MQ=255)

                                                     tGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGaa                                 <  1:2163669/93‑1 (MQ=255)

                                                     tGGGAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCATGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGaa                                 >  1:1295318/1‑93 (MQ=255)

                                                      gggAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAg                                <  1:1100836/93‑1 (MQ=255)

                                                       ggAAGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAgg                               <  2:1939656/93‑1 (MQ=255)

                                                          cggCGATATCAAACGCGGTCAGGGCACAGCATCCACCGACAGCGGCGTGCTTAACCACACGCCGTATGGATTTACACCGCGTTTTGCAGGCgc                            >  2:2328890/2‑92 (MQ=255)

                                                          aGGCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGa                            <  2:234452/93‑1 (MQ=255)

                                                            gCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGaaa                          >  1:1107559/1‑93 (MQ=255)

                                                            gCGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGaaa                          >  1:154955/1‑93 (MQ=255)

                                                             cGATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGaaaa                         <  2:4145246/93‑1 (MQ=255)

                                                              gATATCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGaaaaa                        <  2:4328431/93‑1 (MQ=255)

                                                               atatCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACCGGCGTATGGATTTAACACGCGTTTTGACGGCGAAAAAg                       >  2:3975748/1‑93 (MQ=255)

                                                                 atCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACAAGCGTTTTGAAGGCGAAAAAGGa                     >  2:4860/1‑93 (MQ=255)

                                                                  tCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGaa                    >  2:2859349/1‑93 (MQ=255)

                                                                  tCAAACGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGaa                    <  1:84737/93‑1 (MQ=255)

                                                                     aaCGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGaacca                 >  1:30641/1‑93 (MQ=255)

                                                                     aaCGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGaacca                 >  1:2445210/1‑93 (MQ=255)

                                                                     aaCGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGaacca                 <  1:4318686/93‑1 (MQ=255)

                                                                     aaCGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGaacca                 >  2:1792666/1‑93 (MQ=255)

                                                                      aCGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGaaccaa                <  2:1629041/93‑1 (MQ=255)

                                                                      aCGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGaaccaa                <  2:1629043/93‑1 (MQ=255)

                                                                      aCGCGGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGaaccaa                <  1:3254461/93‑1 (MQ=255)

                                                                         cgGGAAGGGAACAGCATCCACCGAGAGTGGCGTTCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAAccc             <  2:3417682/93‑1 (MQ=255)

                                                                         cgGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAAccc             <  2:567104/93‑1 (MQ=255)

                                                                         cgGGAAGGGAACAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAAccc             <  1:3167162/93‑1 (MQ=255)

                                                                          gggaagggaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCt            >  1:4292973/1‑93 (MQ=255)

                                                                           ggaagggaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTg           >  2:1378920/1‑93 (MQ=255)

                                                                             aagggaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAAccccgga         >  2:3620518/1‑89 (MQ=255)

                                                                              agggaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTgaag        >  2:3305023/1‑93 (MQ=255)

                                                                              agggaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTgaag        >  1:2197057/1‑93 (MQ=255)

                                                                               gggaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTgaaga       >  2:2744749/1‑93 (MQ=255)

                                                                               gggaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGCTTTGAAGGCGAAAAAGGAACCAACCCTgaaga       >  1:2163201/1‑93 (MQ=255)

                                                                                ggaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTgaagaa      <  2:1305260/93‑1 (MQ=255)

                                                                                ggaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTcaagaa      <  2:2970109/93‑1 (MQ=255)

                                                                                 gaaCAGTATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTGAAGAAc     <  1:1035999/93‑1 (MQ=255)

                                                                                 gaaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTGAAGAAc     <  1:410824/93‑1 (MQ=255)

                                                                                  aaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTGAAGAACt    <  1:1955223/93‑1 (MQ=255)

                                                                                  aaCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTGAAGAACt    >  2:20415/1‑93 (MQ=255)

                                                                                   aCAGCATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTGAAGAACTg   <  2:167977/93‑1 (MQ=255)

                                                                                   aCAGCATCCACCGAGAGTGGCATGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTGAAGAACTg   <  1:2508651/93‑1 (MQ=255)

                                                                                    cAGCATCCCCCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTGAAGAACTGa  <  2:1104279/93‑1 (MQ=255)

TTTAAGCCCACAGGAGAGCAACAATGACAATCCATAAGAAAGGTCAGGCACACTGGGAAGGCGATATCAAACGCGGGAAGGGAACAGTATCCACCGAGAGTGGCGTGCTGAACCAACAGCCGTATGGATTTAACACGCGTTTTGAAGGCGAAAAAGGAACCAACCCTGAAGAACTGA  >  NC_000913/1554626‑1554802

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