Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,512 G→A 18.8% M1M (GTG→ATG) †‡ nrfE → heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,291,5120GA18.8% 452.2 / 81.2 200M1M (GTG→ATG) †nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (80/82);  new base A (15/23);  total (95/105)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

GACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGC  >  NC_000913/4291413‑4291611
                                                                                                   |                                                                                                   
 aCAGTTGACGTTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGTAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACg                                                                                                     <  1:4235907/99‑1 (MQ=255)
    gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTGCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTGAACGTGtt                                                                                                 >  1:1019962/1‑100 (MQ=255)
    gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                 >  2:2785678/1‑100 (MQ=255)
    gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                 >  1:1561356/1‑100 (MQ=255)
     ttGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTg                                                                                                <  2:649170/100‑1 (MQ=255)
     ttGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTg                                                                                                <  2:293221/100‑1 (MQ=255)
      tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                               >  1:2636657/1‑100 (MQ=255)
      tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                               <  2:29685/100‑1 (MQ=255)
       gACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGCCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGGTGGTGTTAAGTCTCGGGGTCAACATGTTGAc                                                                                              >  1:787336/1‑100 (MQ=255)
       gACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAc                                                                                              <  2:794864/100‑1 (MQ=255)
         cGGTTGCGCAACCCAGAAACGGGGTGGTTTCTGGTCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAccc                                                                                            <  1:2189574/100‑1 (MQ=255)
          ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                           <  2:771593/100‑1 (MQ=255)
          ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                           <  2:1296699/100‑1 (MQ=255)
           gTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCg                                                                                          <  1:3509127/100‑1 (MQ=255)
           gTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCg                                                                                          <  2:3889659/100‑1 (MQ=255)
             ggCTTAAGACAGAAAAGAGGTGGTTTTTGGACGTATTCCTGCCGGACGTCGGTTTGCTGGCCTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgag                                                                                        <  1:708628/100‑3 (MQ=255)
             ggCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGtt                                                                                        >  1:3384345/1‑100 (MQ=255)
             ggCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGtt                                                                                        <  2:3263133/100‑1 (MQ=255)
              gCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTg                                                                                       <  2:4166774/100‑1 (MQ=255)
               cGTAAGCCAGATAAGAGGTGTTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTGGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGa                                                                                      >  1:1086403/1‑100 (MQ=255)
               cGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagac                                                                                      <  2:2599550/100‑5 (MQ=39)
               cGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagac                                                                                      <  2:2246540/100‑5 (MQ=39)
               cGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGGCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGa                                                                                      <  2:1531794/100‑1 (MQ=255)
                 tagtccAGCAAAGAGGTGGTTTCTGGTTGTATTCTTTCCTGAAGTCCGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagac                                                                                      <  1:3803405/94‑5 (MQ=25)
                  aaGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                   >  2:2075615/1‑100 (MQ=255)
                  aaGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                   <  2:1697573/100‑1 (MQ=255)
                  aaGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                   >  1:3286547/1‑100 (MQ=255)
                  aaGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                   <  2:853476/100‑1 (MQ=255)
                   aGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGc                                                                                  >  2:2743732/1‑100 (MQ=255)
                      cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCtt                                                                               >  1:852702/1‑100 (MQ=255)
                       aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTc                                                                              >  1:2740569/1‑100 (MQ=255)
                        gAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTcg                                                                             <  2:3293502/100‑1 (MQ=255)
                          aaaGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTcgcg                                                                           >  1:1812885/1‑100 (MQ=255)
                           aaGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                          <  2:3476354/100‑1 (MQ=255)
                           aaGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTGGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                          >  1:745248/1‑100 (MQ=255)
                            ttagGTGGTTTCTTGATGTATTCGTTCCTGAAGTCAGTTTTCTGGCGTTGTTGTTGAGTCTCGGGGTCAACGTGTTAACCCCGTTGACgg                                                                                   <  1:169648/88‑1 (MQ=255)
                              agGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                       <  1:3026691/100‑1 (MQ=255)
                               ggtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                      >  1:1495974/1‑100 (MQ=255)
                                 tggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                    <  2:4018962/100‑1 (MQ=255)
                                 tggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                    <  2:884026/100‑1 (MQ=255)
                                 tggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                    >  1:1779458/1‑100 (MQ=255)
                                  ggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCg                                                                   >  1:909254/1‑100 (MQ=255)
                                     ttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                 >  1:3478229/1‑99 (MQ=255)
                                      tCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTg                                                               >  1:2925709/1‑100 (MQ=255)
                                       cTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgc                                                              <  1:4153737/100‑1 (MQ=255)
                                         ggACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                            >  2:733511/1‑100 (MQ=255)
                                           aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                          <  2:4253228/100‑1 (MQ=255)
                                           aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                          <  1:3946832/100‑1 (MQ=255)
                                             gTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTGACGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTgg                                                         <  2:2140836/99‑1 (MQ=255)
                                              tATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                       <  1:1591783/100‑1 (MQ=255)
                                              tATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                       <  1:1718529/100‑1 (MQ=255)
                                              tATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                       <  2:1054908/100‑1 (MQ=255)
                                              tATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                       <  2:1495979/100‑1 (MQ=255)
                                                 tccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGc                                                    >  1:2064092/1‑100 (MQ=255)
                                                 tccttcctGAAGTCGGTTGTCTGGCGTTGTTGTTAAGTCTCGGGGTCACCGTGTTGACGCCGTTGACGGCCTTCGCGGGAGTGCGGTGGCGCTGGCTTGc                                                    >  1:200739/1‑100 (MQ=255)
                                                  ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGc                                                        <  1:50965/95‑1 (MQ=255)
                                                   cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCa                                                  <  1:3243452/100‑1 (MQ=255)
                                                     tcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatg                                                <  2:3144450/100‑1 (MQ=255)
                                                       ctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgat                                              >  2:3841584/1‑100 (MQ=255)
                                                       ctGAAGTCGGTTTTCTGGCGTTGTTGTTAAAACTCAGCACCAACGTGTTGACCCCGTTGACGACCTTCGCGGGAGTTCGGTTTCGCTGGCCTCCCatta                                               >  1:1951497/1‑97 (MQ=255)
                                                        tGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                             <  2:552026/100‑1 (MQ=255)
                                                         gAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGc                                            <  1:420872/100‑1 (MQ=255)
                                                         gAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGc                                            <  1:2252237/100‑1 (MQ=255)
                                                            gTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                         <  2:4383016/100‑1 (MQ=255)
                                                            gTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                         <  1:1795523/100‑1 (MQ=255)
                                                            gTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                         <  2:103133/100‑1 (MQ=255)
                                                              cGGTTTTCGGGGGTTGTTGTTGAGGCTCGGGGTCAACGTGTTGGCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCAGGATGCGACt                                        <  1:3785754/99‑1 (MQ=255)
                                                               ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                      >  1:1184304/1‑100 (MQ=255)
                                                               ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                      >  2:605614/1‑100 (MQ=255)
                                                               ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                      >  1:1438714/1‑100 (MQ=255)
                                                                gTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTAAc                                     >  1:3829031/1‑100 (MQ=255)
                                                                  tttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACtaaatt                                   >  2:2039617/1‑95 (MQ=255)
                                                                   ttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTg                                  <  2:492700/100‑1 (MQ=255)
                                                                    tCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGc                                 >  2:1845769/1‑100 (MQ=255)
                                                                        gCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCg                             <  2:3777646/100‑1 (MQ=255)
                                                                            tgttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCCGCCATGATGCGACTCACTTGCATCGGCAt                         <  2:3277099/100‑1 (MQ=255)
                                                                               tgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCt                      >  1:1305143/1‑100 (MQ=255)
                                                                               tgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCt                      >  1:1261225/1‑100 (MQ=255)
                                                                               tgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCt                      >  2:1311005/1‑100 (MQ=255)
                                                                                gttAAGTCTCGGGGGCAACGAGTTGACTCGGTTGAAGGCCTTCGCGGGAGGACGATTGCACTGGCCTGACATGATGCGACTCACTGGCaggggaa                          >  1:2671461/1‑89 (MQ=255)
                                                                                  tAAGTCTCGGGGTCAACGTGTTGTCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                   <  1:3147680/100‑1 (MQ=255)
                                                                                    gtcaacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                 >  2:2650227/6‑100 (MQ=39)
                                                                                    gtcaacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                 >  2:3713435/6‑100 (MQ=39)
                                                                                    aGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                 >  2:3678856/1‑100 (MQ=255)
                                                                                       ctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGt              >  1:2165285/1‑100 (MQ=255)
                                                                                       aacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGt              >  2:255974/3‑100 (MQ=255)
                                                                                        tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt             <  2:3522946/100‑1 (MQ=255)
                                                                                        acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt             >  1:115599/2‑100 (MQ=255)
                                                                                         cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc            >  2:3638160/1‑100 (MQ=255)
                                                                                         cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGATTACGCTGGCCTGCCATGATGAGACTCACTTGCATCGGCATTCTGGAGCAGTTc            >  1:190311/1‑100 (MQ=255)
                                                                                          ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcg           >  2:4002371/1‑100 (MQ=255)
                                                                                             gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCgcgc        >  2:3113119/1‑100 (MQ=255)
                                                                                               cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc      >  1:467998/1‑100 (MQ=255)
                                                                                               cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc      >  1:2980944/1‑100 (MQ=255)
                                                                                               cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCAATCTGGCGAAGTTCACGCta      >  2:1461233/1‑99 (MQ=255)
                                                                                               cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGCGCGGTTGCGCTGGCCTGCGATGATGCGACTCGCTTGGATCGGCCGTCTGgcg                  >  1:3629344/1‑88 (MQ=255)
                                                                                                 aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct    <  1:369936/100‑1 (MQ=255)
                                                                                                   gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgc  <  1:1845763/100‑1 (MQ=255)
                                                                                                   gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgc  <  2:1215766/100‑1 (MQ=255)
                                                                                                   gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgc  >  2:1466007/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGC  >  NC_000913/4291413‑4291611

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