Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,289,526 A→G 100% K94K (AAA→AAG nrfB → nitrite reductase, formate‑dependent, penta‑heme cytochrome c

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,289,5260AG94.4% 332.2 / 5.6 106K94K (AAA→AAGnrfBnitrite reductase, formate‑dependent, penta‑heme cytochrome c
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (58/42);  total (61/45)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 106 total aligned reads displayed.

AAAACATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTC  >  NC_000913/4289430‑4289622
                                                                                                |                                                                                                
aaaaCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCACCCATCACCGCACCCCCGCGAAGGGGAGAAGGAt                                                                                               >  2:1285906/1‑100 (MQ=255)
  aaCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtgt                                                                                             <  1:219613/100‑1 (MQ=255)
   aCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtgtg                                                                                            <  1:372509/100‑1 (MQ=255)
    cATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGa                                                                                           <  2:799268/100‑1 (MQ=255)
      tGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGAt                                                                                          <  2:1672966/99‑1 (MQ=255)
        cTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATgcg                                                                                       >  2:2119407/1‑100 (MQ=255)
        cTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATgcg                                                                                       >  1:2211541/1‑100 (MQ=255)
        cTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATgcg                                                                                       >  1:2115963/1‑100 (MQ=255)
         ttCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATgcgc                                                                                      <  1:1870879/100‑1 (MQ=255)
          tCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCt                                                                                     <  2:1334214/100‑1 (MQ=255)
           ccGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCtt                                                                                    >  1:1254654/1‑100 (MQ=255)
           ccGTTAGCAACCCGAATGACAAACTGCCCGTCACCTGCACCAACTGGCACGGCCAGCCATCACCGCAACACCGCGAATGGGTGAAGGATGTGATGCGCtt                                                                                    <  2:2559451/100‑1 (MQ=255)
            cGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGCAGGATGTGATGCGCttt                                                                                   >  2:723952/1‑100 (MQ=255)
              ttATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTaa                                                                                 <  1:512465/100‑1 (MQ=255)
               tATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAAc                                                                                <  1:1403251/100‑1 (MQ=255)
               tATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAAc                                                                                <  1:2344490/100‑1 (MQ=255)
               tATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAAc                                                                                >  1:1017195/1‑100 (MQ=255)
                aTCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACg                                                                               >  2:2442617/1‑100 (MQ=255)
                aTCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACCGTCACGGCCAGCCATCACCGCAACACCGCAAAGGGGTGAAGGATGTGATGCGCTTTAACg                                                                               >  2:1072760/1‑100 (MQ=255)
                  cAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAg                                                                             >  2:1853292/1‑100 (MQ=255)
                  cAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAg                                                                             >  2:1848566/1‑100 (MQ=255)
                   aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGc                                                                            >  1:2502944/1‑100 (MQ=255)
                   aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAACCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGc                                                                            >  2:1071235/1‑100 (MQ=255)
                   aaCCCGAATAACAAACTGCCGGGCACCTGCTCCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAAGTGATACGCTTTAACGacc                                                                            >  2:2479742/1‑98 (MQ=255)
                    aCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGcc                                                                           >  2:19862/1‑100 (MQ=255)
                     cccGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCg                                                                          >  1:2106830/1‑100 (MQ=255)
                        gAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATg                                                                       <  1:705316/100‑1 (MQ=255)
                        gAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATg                                                                       <  2:501087/100‑1 (MQ=255)
                        gAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATg                                                                       >  2:968700/1‑100 (MQ=255)
                           tAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTAc                                                                    >  1:531522/1‑100 (MQ=255)
                           tAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTAc                                                                    >  1:491445/1‑100 (MQ=255)
                           tAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTAc                                                                    >  2:1189743/1‑100 (MQ=255)
                            aaCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACa                                                                   <  1:660464/100‑1 (MQ=255)
                             aCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACaa                                                                  >  2:1835177/1‑100 (MQ=255)
                             aCAAACTGCCGGTCACATGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACaa                                                                  >  1:1798170/1‑100 (MQ=255)
                              cAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAg                                                                 >  1:746917/1‑100 (MQ=255)
                              cAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAg                                                                 >  1:2422260/1‑100 (MQ=255)
                              cAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAg                                                                 >  1:213252/1‑100 (MQ=255)
                              cAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAg                                                                 <  2:2262872/100‑1 (MQ=255)
                              cAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAg                                                                 >  1:2040638/1‑100 (MQ=255)
                               aaaCTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAgg                                                                <  2:891172/100‑1 (MQ=255)
                                aaCTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGt                                                               <  2:698911/100‑1 (MQ=255)
                                aaCTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGt                                                               <  1:1853306/100‑1 (MQ=255)
                                  cTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTg                                                             <  2:1115433/100‑1 (MQ=255)
                                  cTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCACAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTg                                                             >  1:1477414/1‑100 (MQ=255)
                                   tGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTgg                                                            <  2:1202057/100‑1 (MQ=255)
                                   tGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTgg                                                            <  1:985529/100‑1 (MQ=255)
                                    gCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGga                                                           >  1:2398763/1‑100 (MQ=255)
                                      cGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGgaga                                                         >  1:2061309/1‑100 (MQ=255)
                                      cGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGgaga                                                         <  2:1051139/100‑1 (MQ=255)
                                       ggTCACCTGCACCAACTGTCACTGCCAGCAATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGagag                                                        <  1:1876559/100‑1 (MQ=255)
                                       ggTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGagag                                                        <  2:27409/100‑1 (MQ=255)
                                        gTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGc                                                       <  2:2211526/100‑1 (MQ=255)
                                         tCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCa                                                      <  2:369100/100‑1 (MQ=255)
                                           aCCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGa                                                    >  2:120796/1‑100 (MQ=255)
                                            ccTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGaa                                                   >  1:668534/1‑100 (MQ=255)
                                              tGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACa                                                 >  1:2151770/1‑100 (MQ=255)
                                                cACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACCGc                                               >  1:2361486/1‑100 (MQ=255)
                                                  ccAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGt                                             <  1:1802753/100‑1 (MQ=255)
                                                   cAACTGGCACGGCCAGCGATCACCGCAACACCGCGAAGGGGTGAAGGGTGTGATGAGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTc                                            <  1:284131/100‑1 (MQ=255)
                                                     aCTGTCACGGCCAGCCATCACCGCAACACTGCGAAGTGGTGAAGGATGTGATGCGCTTTAACGAGCCTATGTACAAGGTTGGAGAGCAGAACAGCGTCt                                           <  2:2407156/99‑1 (MQ=255)
                                                     aCTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTg                                          >  1:954822/1‑100 (MQ=255)
                                                     aCTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTg                                          <  1:2173824/100‑1 (MQ=255)
                                                     aCTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTg                                          <  1:1383143/100‑1 (MQ=255)
                                                     aCTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTg                                          <  1:1586858/100‑1 (MQ=255)
                                                      cTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGt                                         <  2:139874/100‑1 (MQ=255)
                                                        gTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTAt                                       >  2:2018052/1‑100 (MQ=255)
                                                          cagggCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGGGCATGAACTAGCCGATGTACAAGGAAGGAGAGCAGAACAGCGTCCGt                                         >  2:1433230/4‑96 (MQ=255)
                                                          cACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGt                                     >  1:1681404/1‑100 (MQ=255)
                                                              gccagccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTg                                 >  2:2039180/1‑100 (MQ=255)
                                                               ccagccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGt                                >  1:1525562/1‑100 (MQ=255)
                                                               ccagccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGt                                >  1:1966475/1‑100 (MQ=255)
                                                               ccagccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGt                                <  1:1650336/100‑1 (MQ=255)
                                                                cagccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTc                               >  1:2000864/1‑100 (MQ=255)
                                                                  gccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCAt                             <  1:1331343/100‑1 (MQ=255)
                                                                    caTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCt                           >  1:13826/1‑100 (MQ=255)
                                                                       cACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGcc                        >  1:563020/1‑100 (MQ=255)
                                                                       cACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGcc                        <  1:1070921/100‑1 (MQ=255)
                                                                       cACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGcc                        >  2:1211454/1‑100 (MQ=255)
                                                                       cACCGCAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGcc                        <  2:44248/100‑1 (MQ=255)
                                                                        aCCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCg                       >  2:1318775/1‑100 (MQ=255)
                                                                        aCCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCg                       >  2:154834/1‑100 (MQ=255)
                                                                         ccGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCgg                      <  1:968708/100‑1 (MQ=255)
                                                                          cGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGa                     >  1:904557/1‑100 (MQ=255)
                                                                            cAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAAc                   >  1:23212/1‑100 (MQ=255)
                                                                            cAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAAc                   <  1:745357/100‑1 (MQ=255)
                                                                            cAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGGTGGAGAGCAGAACGGCGTCTGTATGTCCTGTCATCTGCCGGAAc                   >  1:117095/1‑100 (MQ=255)
                                                                             aaCACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACATCGTCTGTATGTGCTGTCATCTGCAGGAACa                  >  2:349550/1‑100 (MQ=255)
                                                                              acacCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAg                 >  1:1268507/1‑100 (MQ=255)
                                                                                acCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGt                >  2:2154018/1‑99 (MQ=255)
                                                                                 ccGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTg              >  2:2105564/1‑100 (MQ=255)
                                                                                 ccGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTg              <  2:2502929/100‑1 (MQ=255)
                                                                                  cgcgAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGc             >  2:427157/1‑100 (MQ=255)
                                                                                  cgcgAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGc             >  1:956414/1‑100 (MQ=255)
                                                                                  cgcgAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGc             <  2:1669171/100‑1 (MQ=255)
                                                                                    cgAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCaa           <  2:2119521/100‑1 (MQ=255)
                                                                                     gAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCaaa          >  2:1776049/1‑100 (MQ=255)
                                                                                        ggggTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAg       <  1:378445/100‑1 (MQ=255)
                                                                                           gTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGt    <  1:1587074/100‑1 (MQ=255)
                                                                                             gAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTc  >  2:1838431/1‑100 (MQ=255)
                                                                                                |                                                                                                
AAAACATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTC  >  NC_000913/4289430‑4289622

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: