Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,287,458 G→A 19.5% intergenic (‑87/‑306) acs ← / → nrfA acetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,4580GA19.5% 297.8 / 61.1 140intergenic (‑87/‑306)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base G (65/47);  new base A (19/8);  total (84/56)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-02

Only 100 of 141 total aligned reads displayed.

GAATTTGGCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTT  >  NC_000913/4287359‑4287555
                                                                                                   |                                                                                                 
gAATTTGGCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCg                                                                                                   <  2:1690466/100‑1 (MQ=255)
 aaTTTGGCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGa                                                                                                  >  1:251300/1‑100 (MQ=255)
 aaTTTGGCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGa                                                                                                  >  1:539150/1‑100 (MQ=255)
   tttGGCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACa                                                                                                <  1:827383/100‑1 (MQ=255)
   ttaggCTCATGCTTTTGTTCTCCTTGTATGATGTTAATAATATGTGGTATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTTTTTAGTATTTGGGCGACa                                                                                                <  2:48154/97‑1 (MQ=255)
    ttGGCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAg                                                                                               >  1:630179/1‑100 (MQ=255)
      ggCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGa                                                                                              >  2:441176/1‑99 (MQ=255)
      ggCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGAt                                                                                             >  2:509430/1‑100 (MQ=255)
       gCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATc                                                                                            >  1:1344471/1‑100 (MQ=255)
        cTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCa                                                                                           >  2:1316884/1‑100 (MQ=255)
          cATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACg                                                                                         >  2:1882093/1‑100 (MQ=255)
           aTGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGc                                                                                        <  2:384343/100‑1 (MQ=255)
           aTGCTTTTGTTCTCCTTGAAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGACAGTTTTGTTTAGTATTTGGGCGAAAGATCACg                                                                                         >  2:239820/1‑99 (MQ=255)
             gCTTTTGTTCTCCTTGTAGGAAGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATAGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCaa                                                                                      >  2:803513/1‑100 (MQ=255)
              cTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCaaa                                                                                     >  1:1852389/1‑100 (MQ=255)
              cTTTTGTTCTCCTTGCAGGATGTTAATAATATGTGGCATACGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCa                                                                                       >  2:1606611/1‑98 (MQ=255)
                tttGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAg                                                                                   <  2:779057/100‑1 (MQ=255)
                tttGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAg                                                                                   <  1:477248/100‑1 (MQ=255)
                 ttGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGt                                                                                  <  2:1968695/100‑1 (MQ=255)
                    ttCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGa                                                                               <  1:1026225/100‑1 (MQ=255)
                      ctcCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAAt                                                                             >  2:16249/1‑100 (MQ=255)
                        ncTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAAtt                                                                            >  1:238504/2‑99 (MQ=255)
                         cTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATtgt                                                                          >  2:245669/1‑100 (MQ=255)
                          ttGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATtgtg                                                                         <  1:2028482/100‑1 (MQ=255)
                          ttGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATtgtg                                                                         <  1:1448369/100‑1 (MQ=255)
                           tGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGc                                                                        >  1:567771/1‑100 (MQ=255)
                           tGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGc                                                                        >  2:2121544/1‑100 (MQ=255)
                              aGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCaaa                                                                     >  2:1822389/1‑100 (MQ=255)
                               ggATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAAt                                                                    >  2:628070/1‑100 (MQ=255)
                               ggATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAAt                                                                    >  1:2245048/1‑100 (MQ=255)
                               ggATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAAt                                                                    <  1:1092925/100‑1 (MQ=255)
                               ggATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAAt                                                                    >  2:662710/1‑100 (MQ=255)
                                gATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATa                                                                   <  2:2204026/100‑1 (MQ=255)
                                 aTGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGACGACAGATCACGCAGAAGTAGAATTGTGCAAAt                                                                    <  2:741570/98‑1 (MQ=255)
                                  tGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGGCAGATCACGCAAAAGTAGAATTGTGCAAATaaa                                                                 >  2:1525605/1‑100 (MQ=255)
                                  tGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGATTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATaaa                                                                 >  2:2236346/1‑100 (MQ=255)
                                    ttAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACg                                                               <  2:622876/100‑1 (MQ=255)
                                     taataaTATGTGTGATAAGCGTTAAATTTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACg                                                               <  2:1915844/99‑1 (MQ=255)
                                     taataaTATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACgg                                                              >  1:2127719/1‑100 (MQ=255)
                                       ataataTGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCa                                                            >  2:689214/1‑100 (MQ=255)
                                              gtgGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAAt                                                     >  2:1647911/1‑100 (MQ=255)
                                               tgGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAAtt                                                    >  1:24114/1‑100 (MQ=255)
                                                ggCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAAttt                                                   >  2:275999/1‑100 (MQ=255)
                                                  cATAAGCGTCAACTGTAGGGGTATCGGAAGTATTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAAttttt                                                 >  2:1667993/1‑100 (MQ=255)
                                                   aTAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTg                                                <  1:639057/100‑1 (MQ=255)
                                                    tAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGa                                               <  1:215372/100‑1 (MQ=255)
                                                    tAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGa                                               <  1:2390095/100‑1 (MQ=255)
                                                     aaGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGaa                                              >  2:1163096/1‑100 (MQ=255)
                                                     aaGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGaa                                              <  2:2293102/100‑1 (MQ=255)
                                                      aGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCCGGGTAATTTTTGaa                                              >  1:263571/1‑99 (MQ=255)
                                                       gCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAgg                                            <  1:939610/100‑1 (MQ=255)
                                                       gCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAgg                                            <  1:1947315/100‑1 (MQ=255)
                                                       gCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGAGCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAgg                                            <  1:18466/100‑1 (MQ=255)
                                                         gTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTc                                          >  2:1493053/1‑100 (MQ=255)
                                                            aaaTGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGc                                       >  2:899784/1‑100 (MQ=255)
                                                            aaaTGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGc                                       <  2:1488912/100‑1 (MQ=255)
                                                            aaaTGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGc                                       >  1:1705880/1‑100 (MQ=255)
                                                            aaaTGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGc                                       >  1:1890272/1‑100 (MQ=255)
                                                              aTGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTAAAGGTCAGcaa                                     >  1:1266810/1‑100 (MQ=255)
                                                               tGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGcaac                                    <  2:1802806/100‑1 (MQ=255)
                                                                 tAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGcaacaa                                  >  1:103596/1‑100 (MQ=255)
                                                                     ggTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAg                               >  2:394337/1‑99 (MQ=255)
                                                                     ggTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAACAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGt                              <  1:1769704/100‑1 (MQ=255)
                                                                      gTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGtt                             <  1:2277736/100‑1 (MQ=255)
                                                                         ttGGTAGTTTTGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGAt                          <  1:703952/100‑1 (MQ=255)
                                                                         ttGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGAt                          <  2:1503354/100‑1 (MQ=255)
                                                                         ttGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGAt                          >  1:1215170/1‑100 (MQ=255)
                                                                           ggCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTa                        >  2:934963/1‑100 (MQ=255)
                                                                           ggCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTa                        >  1:808979/1‑100 (MQ=255)
                                                                               gTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTc                    >  1:1588237/1‑100 (MQ=255)
                                                                               gTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTc                    >  1:1505398/1‑100 (MQ=255)
                                                                               gTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTc                    <  1:2336325/100‑1 (MQ=255)
                                                                                 tttGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCtt                  <  1:887169/100‑1 (MQ=255)
                                                                                  ttGTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCttt                 >  1:2291628/1‑100 (MQ=255)
                                                                                   tGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTg                <  1:1176125/100‑1 (MQ=255)
                                                                                    gTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCttt                 >  1:404690/1‑98 (MQ=255)
                                                                                    gTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGa               >  2:916859/1‑100 (MQ=255)
                                                                                    gTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGa               <  1:2053534/100‑1 (MQ=255)
                                                                                    gTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGa               >  2:42393/1‑100 (MQ=255)
                                                                                    gTTTAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTATTTGa               >  2:897829/1‑100 (MQ=255)
                                                                                     tttAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAg              <  1:475686/100‑1 (MQ=255)
                                                                                     tttAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAg              <  2:567768/100‑1 (MQ=255)
                                                                                     tttAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAg              <  1:2119431/100‑1 (MQ=255)
                                                                                     tttAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAg              >  1:1398770/1‑100 (MQ=255)
                                                                                       tAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGa               >  2:966920/1‑97 (MQ=255)
                                                                                       tAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGa            >  2:1174595/1‑100 (MQ=255)
                                                                                       tAGTATTTGGGCAACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGa            <  1:1532197/100‑1 (MQ=255)
                                                                                        aGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGaa           >  1:1967188/1‑100 (MQ=255)
                                                                                        aGTATTTGGGCCACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCATCAACAAAAGTTGATTAATTCTTTGAGGa            <  1:966931/99‑1 (MQ=255)
                                                                                         gTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAAc          >  2:957070/1‑100 (MQ=255)
                                                                                          tATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACa         <  2:27337/100‑1 (MQ=255)
                                                                                          tATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACa         >  1:2184272/1‑100 (MQ=255)
                                                                                          tATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACa         <  1:16250/100‑1 (MQ=255)
                                                                                          tATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGACAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACa         >  2:978536/1‑100 (MQ=255)
                                                                                             ttGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGc      >  2:2018115/1‑100 (MQ=255)
                                                                                              tGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCa     >  2:2389670/1‑100 (MQ=255)
                                                                                               gggCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAg    <  2:1567123/100‑1 (MQ=255)
                                                                                                ggCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGt   <  1:394340/100‑1 (MQ=255)
                                                                                                ggCGACAGATCACGCAAAAGTAGAATTGTGAAAATAAACGGCAGGGTAAGTTTTGAAGGACAGCAACAAAAGTGGAATAATTCTTTGAGGAACATGCa     >  1:600138/1‑98 (MQ=255)
                                                                                                 gCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGtt  >  2:1146814/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
GAATTTGGCTCATGCTTTTGTTCTCCTTGTAGGATGTTAATAATATGTGGCATAAGCGTTAAATGTAGGGGTATTGGCAGTTTTGTTTAGTATTTGGGCGACAGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTT  >  NC_000913/4287359‑4287555

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: