BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,046,716	A→C	intergenic (‑114/‑250)	yihG ← / → polA	inner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,046,716	0	A	C	98.8%	584.9 / ‑4.4	171	intergenic (‑114/‑250)	yihG/polA	inner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease
Reads supporting (aligned to +/- strand): ref base A (1/1); new base C (84/85); total (85/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

AGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAT  >  NC_000913/4046620‑4046812

agCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATCTCTATCCACATTCCAAAAAATCCCCTTa                                                                                               >  2:1339125/1‑100 (MQ=255)

agCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTa                                                                                               >  1:522360/1‑100 (MQ=255)

agCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTa                                                                                               <  1:1567359/100‑1 (MQ=255)

 gcgcAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTAt                                                                                              >  2:1554291/1‑100 (MQ=255)

   gcAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTAtct                                                                                            >  2:753771/1‑100 (MQ=255)

   gcAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTAtct                                                                                            >  2:2432938/1‑100 (MQ=255)

     aaGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCa                                                                                          <  2:3285424/100‑1 (MQ=255)

      aGGCGCAGATTTTAGCAGAAAACGACCCAAATGACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCaa                                                                                         <  1:1009727/100‑1 (MQ=255)

      aGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCaa                                                                                         >  1:591112/1‑100 (MQ=255)

        gcgcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAatt                                                                                       >  1:2904457/1‑100 (MQ=255)

        gcgcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAatt                                                                                       >  1:3900103/1‑100 (MQ=255)

        gcgcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAatt                                                                                       <  2:1240599/100‑1 (MQ=255)

         cgcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAatta                                                                                      >  2:103763/1‑100 (MQ=255)

          gcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAattat                                                                                     >  2:760082/1‑100 (MQ=255)

          gcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAattat                                                                                     >  2:3592079/1‑100 (MQ=255)

          gcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAattat                                                                                     <  1:2109218/100‑1 (MQ=255)

          aaaGATTTTAGAAGAAAACGACCCAAATAACAGAAAAACCTTAAGGAGAAAAATAATTCATATCTAAAAAAAAAAAAAAAAATCCCCTTATCTCAattat                                                                                     <  2:3750896/98‑1 (MQ=255)

              aTTTTAGCAGAAAACGACTCCAATAATGGATGATCCTTAAGGAGAGAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAgg                                                                                 <  1:3819418/100‑1 (MQ=255)

              aTTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCGATTATTAgg                                                                                 <  2:2738877/100‑1 (MQ=255)

              aTTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAgg                                                                                 <  2:2023547/100‑1 (MQ=255)

              aTTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAgg                                                                                 <  2:1942481/100‑1 (MQ=255)

              aTTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAgg                                                                                 <  1:3059337/100‑1 (MQ=255)

              aTTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAgg                                                                                 <  1:3338181/100‑1 (MQ=255)

                    gCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatgga                                                                           <  1:2618404/100‑1 (MQ=255)

                     cAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatggat                                                                          >  2:3940504/1‑100 (MQ=255)

                     cAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatggat                                                                          >  1:3837651/1‑100 (MQ=255)

                     cAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatggat                                                                          >  1:2799417/1‑100 (MQ=255)

                     cAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatggat                                                                          >  2:2179434/1‑100 (MQ=255)

                      aGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACCTTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGatt                                                                         >  2:547296/1‑100 (MQ=255)

                       gAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGattt                                                                        >  1:1970537/1‑100 (MQ=255)

                       gAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCACCTTATCTCAATTATTAGGGATGGattt                                                                        >  2:3368794/1‑100 (MQ=255)

                           aCGACCCACAAAACGGATGATCCGTAAGGAGAAAAATACTTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGAATGGatttattt                                                                    <  2:3854903/100‑1 (MQ=255)

                            cGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTAtttt                                                                   <  2:2023994/100‑1 (MQ=255)

                              aCCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTa                                                                 <  1:1830390/100‑1 (MQ=255)

                              aCCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTa                                                                 <  1:3641610/100‑1 (MQ=255)

                               cccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTaa                                                                >  2:1311554/1‑100 (MQ=255)

                                ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAAc                                                               >  2:1399308/1‑100 (MQ=255)

                                 cAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACt                                                              <  2:284953/100‑1 (MQ=255)

                                 cAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACt                                                              >  1:1858449/1‑100 (MQ=255)

                                 cAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACt                                                              >  2:3696405/1‑100 (MQ=255)

                                  aaaTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTg                                                             >  2:60680/1‑100 (MQ=255)

                                    aaaaaggATCATCCTTAAGGAGAAAAATAATTCATATCTATCAACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCa                                                           <  2:265453/95‑1 (MQ=255)

                                    aTAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCa                                                           <  1:1328563/100‑1 (MQ=255)

                                      aaCGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATg                                                         >  2:2136615/1‑100 (MQ=255)

                                      aaCGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATg                                                         <  1:3398727/100‑1 (MQ=255)

                                        cGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGaa                                                       >  1:1532329/1‑100 (MQ=255)

                                         ggATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTCACTGCATGaaa                                                      >  2:1673107/1‑100 (MQ=255)

                                          gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGaaaa                                                     <  1:2634260/100‑1 (MQ=255)

                                          gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGaaaa                                                     <  1:727343/100‑1 (MQ=255)

                                          gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGaaaa                                                     <  2:2860844/100‑1 (MQ=255)

                                            tgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAAc                                                   <  2:1802960/100‑1 (MQ=255)

                                              atCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACaa                                                 >  1:1450542/1‑100 (MQ=255)

                                                 cTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCACCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGa                                              >  2:815621/1‑100 (MQ=255)

                                                    aaGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACaa                                           <  2:1448277/100‑1 (MQ=255)

                                                      ggAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAc                                         <  1:2321652/100‑1 (MQ=255)

                                                      ggAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAc                                         >  1:3078650/1‑100 (MQ=255)

                                                      ggAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAc                                         >  1:589376/1‑100 (MQ=255)

                                                      ggAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAc                                         <  2:1815950/100‑1 (MQ=255)

                                                         gaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatc                                      >  2:3370304/1‑100 (MQ=255)

                                                         gaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatc                                      >  1:1779671/1‑100 (MQ=255)

                                                          aaaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatca                                     >  1:2428795/1‑100 (MQ=255)

                                                          aaaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatca                                     <  2:1871055/100‑1 (MQ=255)

                                                           aaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatcat                                    <  1:1642282/100‑1 (MQ=255)

                                                           aaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatcat                                    >  2:2509051/1‑100 (MQ=255)

                                                           aaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatcat                                    <  1:2800473/100‑1 (MQ=255)

                                                            aaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGATAAACAAAGACAAACATCCTg                                   <  1:3180308/100‑1 (MQ=255)

                                                             aataatTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGc                                  <  2:797152/100‑1 (MQ=255)

                                                               taatTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTg                                <  1:165386/100‑1 (MQ=255)

                                                               taatTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTg                                <  2:561719/100‑1 (MQ=255)

                                                                 atTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTa                              >  2:2596188/1‑100 (MQ=255)

                                                                  ttCAAATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTATTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaa                             >  2:3385070/1‑100 (MQ=255)

                                                                      tatctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAg                         >  1:3494023/1‑100 (MQ=255)

                                                                       atctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGc                        <  1:1987430/100‑1 (MQ=255)

                                                                       atctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGc                        <  2:495478/100‑1 (MQ=255)

                                                                        tctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCa                       >  1:4000887/1‑100 (MQ=255)

                                                                        tctatcCACATTAAAAAAAATCCTCTTATCTCAATTATTAGGGATGGATTTATTTTTAAATGCATACAAAACAAAGACAATCATTATGCTGGAAAAAGCa                       >  2:2289818/1‑100 (MQ=255)

                                                                         ctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCAt                      <  2:3298896/100‑1 (MQ=255)

                                                                          tatcCACATTAGAAAAAATGCCCTTATCTCAATTATTAGGGATGGATTTATTTCTAACTGCATGAAATACAAAGACAAACATCATGCTGTAAAAAGCATg                     <  1:3618162/100‑1 (MQ=255)

                                                                          tatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATg                     >  1:133488/1‑100 (MQ=255)

                                                                           atcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGa                    >  2:2790435/1‑100 (MQ=255)

                                                                            tcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGat                   <  1:827287/100‑1 (MQ=255)

                                                                            tcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGat                   <  2:3634278/100‑1 (MQ=255)

                                                                              cacaTTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGataa                 >  1:3012955/1‑100 (MQ=255)

                                                                               acaTTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGataat                >  1:444367/1‑100 (MQ=255)

                                                                               acaTTAGAAAAAATCCCCTCATCTTAATTATTAGGCATCTATTTATTTTTAACTGCATGAAAAACCAAGACAAATATCCTGCTTTTTACAGCATGatgtt                <  2:646563/100‑4 (MQ=255)

                                                                                    aGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt           <  2:545540/100‑1 (MQ=255)

                                                                                    aGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt           >  2:3001980/1‑100 (MQ=255)

                                                                                    aGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt           >  1:3443211/1‑100 (MQ=255)

                                                                                    aGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt           <  1:1309416/100‑1 (MQ=255)

                                                                                     naaaaaaTCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTa          >  1:3257358/2‑100 (MQ=255)

                                                                                     gAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTa          >  1:1220066/1‑100 (MQ=255)

                                                                                     gAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTa          >  2:2870409/1‑100 (MQ=255)

                                                                                     gAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTa          <  1:195251/100‑1 (MQ=255)

                                                                                     gAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTa          >  1:1979343/1‑100 (MQ=255)

                                                                                     gAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTa          >  2:2215386/1‑100 (MQ=255)

                                                                                       aaaaaTCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTaaa        >  2:1013144/1‑100 (MQ=255)

                                                                                         aagtCCCCTTATCTCAATTATTAGGGATGTATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAg      <  2:2136378/97‑1 (MQ=255)

                                                                                          aaTCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGc     <  2:2475032/100‑1 (MQ=255)

                                                                                             ccccTTATCTCAATTATTAGTGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAt  >  2:305313/1‑100 (MQ=255)

                                                                                             cccATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAt  <  2:2094944/100‑1 (MQ=255)

AGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAT  >  NC_000913/4046620‑4046812

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