Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,513 T→C 17.0% M1A (GTG→GCG) †‡ 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,5130TC17.0% 475.3 / 85.3 216M1A (GTG→GCG) †nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base T (89/90);  new base C (21/16);  total (110/106)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 221 total aligned reads displayed.

CAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCT  >  NC_000913/4291415‑4291612
                                                                                                  |                                                                                                   
cAGTTTGCGGTGGCGTACGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACgtg                                                                                                    <  1:4288533/100‑1 (MQ=255)
  gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                  >  1:901563/1‑100 (MQ=255)
  gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                  >  2:3033850/1‑100 (MQ=255)
    tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                                >  2:87637/1‑100 (MQ=255)
    tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                                <  2:594077/100‑1 (MQ=255)
    tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                                >  2:4260319/1‑100 (MQ=255)
    tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                                >  1:4230839/1‑100 (MQ=255)
       cGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAccc                                                                                             >  2:1448559/1‑100 (MQ=255)
       cGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAccc                                                                                             <  2:3427296/100‑1 (MQ=255)
        ggTGGCGTAAGCCAGAAAAGCGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAcccc                                                                                            <  1:2156630/100‑1 (MQ=255)
        ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAcccc                                                                                            <  2:2820966/100‑1 (MQ=255)
        ggTGGCGTAAGCCAGAAAAGAGGTGGGTTCTGGACGTATTCCTTCCTGAGGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                            <  1:136006/100‑1 (MQ=255)
         gTGGCGTCAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                            <  1:275435/99‑1 (MQ=255)
         gTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCg                                                                                           <  2:1857135/100‑1 (MQ=255)
            gCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTg                                                                                        <  2:163933/100‑1 (MQ=255)
             cGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagac                                                                                       <  2:1981486/100‑5 (MQ=255)
              gTATGCCATAGAAGAGGTTGTTTCTGTACTTATTCCTGCCTGAAGTCGGTTTGTTGGCGTTGTTGTGGAGTCTCGGGGTCAACACGTTGACCCCgag                                                                                         <  1:1508315/97‑3 (MQ=255)
                aaGTCAGAAAAGATGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                    <  2:2237401/100‑1 (MQ=255)
                aaGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                    <  2:2924260/100‑1 (MQ=255)
                  gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGcc                                                                                  <  1:3512263/100‑1 (MQ=255)
                   gcagcagagagGAGGGTTGTGGACGTATTCCTTCCTGCAGTCGGTTTTCTGGCGATGTAGTTACGTCACGAGGGCAACGTGTTGGCCCCGTTGACGGc                                                                                   >  1:2378610/8‑98 (MQ=255)
                   ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCt                                                                                 >  1:2830132/1‑100 (MQ=255)
                   ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCt                                                                                 >  1:3991858/1‑100 (MQ=255)
                    cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCtt                                                                                >  2:2754358/1‑100 (MQ=255)
                     aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTATGGCGTTGTTGTTAAGTCTCGGGGTCCACGTGTTGACCCCGTTGAAGGCCTTc                                                                               >  1:1927401/1‑100 (MQ=255)
                      gAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTcg                                                                              >  2:4071177/1‑100 (MQ=255)
                          agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                          <  2:268755/100‑1 (MQ=255)
                          agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                          <  2:920323/100‑1 (MQ=255)
                          agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGAGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                          >  1:1857268/1‑100 (MQ=255)
                            tggtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                        <  1:444865/99‑1 (MQ=255)
                             ggtggtTTCTGGACGTATGCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                       >  1:648384/1‑100 (MQ=255)
                              gtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTg                                                                      >  2:1404563/1‑100 (MQ=255)
                               ttgtTTCTGGACGTACTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                     <  1:940673/98‑1 (MQ=255)
                                 gtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCgg                                                                   >  2:2467195/1‑100 (MQ=255)
                                  tttCTGGAGCTAGTCGTTTCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTGGGGGTCAACGTTTTTACCCCGTTGACGGCCTTCGCGGGAGTGCgg                                                                   <  1:2894398/99‑1 (MQ=255)
                                  tttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                  <  2:422871/100‑1 (MQ=255)
                                  tttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                  <  1:1031572/100‑1 (MQ=255)
                                    tCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGtt                                                                 >  1:716402/1‑99 (MQ=255)
                                     cTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgc                                                               <  1:4070870/100‑1 (MQ=255)
                                     cTGGACGTATTACTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgc                                                               <  1:313491/100‑1 (MQ=255)
                                      tGGACGTATTCCGTGCTGAAGTCGGTTTTCGGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                              <  1:3364580/100‑1 (MQ=255)
                                         aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                           >  2:2190101/1‑100 (MQ=255)
                                         aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                           <  1:863003/100‑1 (MQ=255)
                                         aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                           <  2:756939/100‑1 (MQ=255)
                                          cGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTgg                                                          >  2:4174416/1‑100 (MQ=255)
                                            tATTCCTTCCTGCAGTCGGTTTTCTGGCGTTGTTGTGAGGTCTCGGAGTCAAGGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGc                                                         <  1:1905083/99‑1 (MQ=255)
                                              ttccttccTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTg                                                      <  2:1796554/100‑1 (MQ=255)
                                               tccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCt                                                       >  1:2794188/1‑98 (MQ=255)
                                                ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                    <  2:1010844/100‑1 (MQ=255)
                                                ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                    <  2:3807620/100‑1 (MQ=255)
                                                  ttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCat                                                  <  1:3571023/100‑1 (MQ=255)
                                                   tcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatg                                                 >  2:1277669/1‑100 (MQ=255)
                                                       gAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGTCTGCCATGATGc                                             <  1:3254813/100‑1 (MQ=255)
                                                         aGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGcag                                           >  1:2879521/1‑98 (MQ=255)
                                                           tCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                          <  1:1431657/99‑1 (MQ=255)
                                                             ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                       <  2:3356040/100‑1 (MQ=255)
                                                              gTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGGCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCAc                                      <  1:303214/100‑1 (MQ=255)
                                                              gTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCAc                                      <  2:4403093/100‑1 (MQ=255)
                                                                tttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                    <  2:178464/100‑1 (MQ=255)
                                                                 ttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTg                                   <  2:1121308/100‑1 (MQ=255)
                                                                    tGGCGTTGTTGTTGAGTGTCTGGGTCAGCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGGGCGGTTGCGCTGGCCTGCCATGATGCGAc                                          <  1:4024231/90‑1 (MQ=255)
                                                                    tGGCGTTGTTGTTGAGTCTCGGGGTCAACGCGTTGACCACGTTGACGGCCTTCGCTGGAGTGCGGTTGCGCTTGCCTGCCATTATGCGACTCACTTGCa                                 <  1:1463898/99‑1 (MQ=255)
                                                                     ggCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGGCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATc                               <  2:3222388/100‑1 (MQ=255)
                                                                     ggCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATc                               >  2:237978/1‑100 (MQ=255)
                                                                     ggCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTAGACGGCCTTTGAGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGTATc                               <  1:437651/100‑1 (MQ=255)
                                                                      gCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCg                              <  2:2130341/100‑1 (MQ=255)
                                                                       cGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                             >  1:1743879/1‑100 (MQ=255)
                                                                       cGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                             <  1:2475623/100‑1 (MQ=255)
                                                                       cGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                             <  2:4155565/100‑1 (MQ=255)
                                                                           gttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCAtt                         >  1:191116/1‑100 (MQ=255)
                                                                            ttgttAAGTCTCGGGGTAAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGACTGCCATGATGCGACTAACTTGAATCGGCAtt                         >  1:1096896/1‑99 (MQ=255)
                                                                              gttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTg                      >  1:723355/1‑100 (MQ=255)
                                                                               ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                     <  2:723356/100‑1 (MQ=255)
                                                                               ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                     <  2:3965908/100‑1 (MQ=255)
                                                                                tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                    >  1:3698134/1‑100 (MQ=255)
                                                                                tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                    >  1:3975566/1‑100 (MQ=255)
                                                                                 aaGTCTCGCGGTCAACGTGTTGACCCCTTTGACGGCCTTCGCGGGCGTGTGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                   <  1:544469/100‑1 (MQ=255)
                                                                                   gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCa                 <  1:1404557/100‑1 (MQ=255)
                                                                                   gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCa                 >  1:2281128/1‑100 (MQ=255)
                                                                                   gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCa                 >  1:1559504/1‑100 (MQ=255)
                                                                                    tctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAg                <  2:1587164/100‑1 (MQ=255)
                                                                                     aacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGt               >  2:352924/3‑100 (MQ=255)
                                                                                      tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGt               >  1:2460758/1‑99 (MQ=255)
                                                                                      acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt              >  2:1713787/2‑100 (MQ=255)
                                                                                      acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt              >  1:44769/2‑100 (MQ=255)
                                                                                      acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt              >  2:3263159/2‑100 (MQ=255)
                                                                                       cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCATTTc             >  2:128668/1‑100 (MQ=255)
                                                                                       cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc             >  2:2412468/1‑100 (MQ=255)
                                                                                       cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc             >  2:288379/1‑100 (MQ=255)
                                                                                       cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc             >  2:600982/1‑100 (MQ=255)
                                                                                         gggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc           >  2:4320253/1‑100 (MQ=255)
                                                                                          ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgcg          <  2:4040140/100‑1 (MQ=255)
                                                                                          ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgcg          <  2:484239/100‑1 (MQ=255)
                                                                                           gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCgcgc         <  1:1281886/100‑1 (MQ=255)
                                                                                            tCAACGTGTTGCCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCt        <  2:648385/100‑1 (MQ=255)
                                                                                             cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc       >  2:2117623/1‑100 (MQ=255)
                                                                                             cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc       >  1:1246761/1‑100 (MQ=255)
                                                                                               aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct     >  1:2038102/1‑100 (MQ=255)
                                                                                                cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctg    >  1:3435461/1‑100 (MQ=255)
                                                                                                  tgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgct  >  2:1919961/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCT  >  NC_000913/4291415‑4291612

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: