Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,512 G→A 16.6% 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,5120GA16.6% 525.6 / 77.8 216M1M (GTG→ATG) †nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (89/91);  new base A (21/15);  total (110/106)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 221 total aligned reads displayed.

GACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTG  >  NC_000913/4291413‑4291610
                                                                                                   |                                                                                                  
gACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACg                                                                                                    <  2:505159/100‑1 (MQ=255)
  cAGTTTGCGGTGGCGTACGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACgtg                                                                                                  <  1:4288533/100‑1 (MQ=255)
    gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                >  1:901563/1‑100 (MQ=255)
      tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                              >  2:87637/1‑100 (MQ=255)
      tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                              <  2:594077/100‑1 (MQ=255)
      tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                              >  1:4230839/1‑100 (MQ=255)
         cGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAccc                                                                                           >  2:1448559/1‑100 (MQ=255)
          ggTGGCGTAAGCCAGAAAAGCGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAcccc                                                                                          <  1:2156630/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)
           gTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCg                                                                                         <  2:3226769/100‑1 (MQ=255)
           gTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCATTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCg                                                                                         >  2:4200738/1‑100 (MQ=255)
              gCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTg                                                                                      <  2:163933/100‑1 (MQ=255)
                gTATGCCATAGAAGAGGTTGTTTCTGTACTTATTCCTGCCTGAAGTCGGTTTGTTGGCGTTGTTGTGGAGTCTCGGGGTCAACACGTTGACCCCgag                                                                                       <  1:1508315/97‑3 (MQ=255)
                    gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGcc                                                                                <  1:3512263/100‑1 (MQ=255)
                    gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGcc                                                                                >  2:2119049/1‑100 (MQ=255)
                    gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTCACCCCGTTGACGGcc                                                                                >  2:4352255/1‑100 (MQ=255)
                     gcagcagagagGAGGGTTGTGGACGTATTCCTTCCTGCAGTCGGTTTTCTGGCGATGTAGTTACGTCACGAGGGCAACGTGTTGGCCCCGTTGACGGc                                                                                 >  1:2378610/8‑98 (MQ=255)
                     ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCt                                                                               >  1:3991858/1‑100 (MQ=255)
                     ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCt                                                                               >  1:2830132/1‑100 (MQ=255)
                     ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCt                                                                               >  2:3138451/1‑100 (MQ=255)
                       aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTATGGCGTTGTTGTTAAGTCTCGGGGTCCACGTGTTGACCCCGTTGAAGGCCTTc                                                                             >  1:1927401/1‑100 (MQ=255)
                            agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                        <  2:920323/100‑1 (MQ=255)
                            agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                        <  2:2081076/100‑1 (MQ=255)
                            agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGAGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                        >  1:1857268/1‑100 (MQ=255)
                            agagGTGGTTTCGGGACGTATTCCTTCCTGAAGTCGGTTGTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGGCGGCCTTCGCggg                                                                        <  2:2177439/100‑1 (MQ=255)
                              tggtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                      <  1:444865/99‑1 (MQ=255)
                               ggtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                      >  2:444866/1‑99 (MQ=255)
                               ggtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                     >  2:3983489/1‑100 (MQ=255)
                               ggtggtTTCTGGACGTATGCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                     >  1:648384/1‑100 (MQ=255)
                                gtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTg                                                                    >  2:1404563/1‑100 (MQ=255)
                                 ttgtTTCTGGACGTACTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                   <  1:940673/98‑1 (MQ=255)
                                  ggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCg                                                                  >  2:4159258/1‑100 (MQ=255)
                                    tttCTGGAGCTAGTCGTTTCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTGGGGGTCAACGTTTTTACCCCGTTGACGGCCTTCGCGGGAGTGCgg                                                                 <  1:2894398/99‑1 (MQ=255)
                                    tttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                <  1:1031572/100‑1 (MQ=255)
                                     ttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGATTGCGGtt                                                               <  2:2312224/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)
                                         ggACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                           >  2:2449728/1‑100 (MQ=255)
                                           aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                         <  2:756939/100‑1 (MQ=255)
                                           aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                         <  1:863003/100‑1 (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)
                                                    ttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCat                                                <  1:3571023/100‑1 (MQ=255)
                                                     tcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatg                                               >  2:1277669/1‑100 (MQ=255)
                                                     tcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatg                                               <  2:4094954/100‑1 (MQ=255)
                                                         gAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGTCTGCCATGATGc                                           <  1:3254813/100‑1 (MQ=255)
                                                         gAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatcc                                           <  2:2991764/100‑3 (MQ=255)
                                                         gAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGc                                           >  2:3468895/1‑100 (MQ=255)
                                                         ctaGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGc                                           <  2:2879538/98‑1 (MQ=255)
                                                           aGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGcag                                         >  1:2879521/1‑98 (MQ=255)
                                                             tCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                        <  1:1431657/99‑1 (MQ=255)
                                                             tCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACt                                       >  2:4043581/1‑100 (MQ=255)
                                                               ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                     >  2:389416/1‑100 (MQ=255)
                                                               ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                     >  2:2385616/1‑100 (MQ=255)
                                                                gTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGGCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCAc                                    <  1:303214/100‑1 (MQ=255)
                                                                gTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCAc                                    >  2:326620/1‑100 (MQ=255)
                                                                  tttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                  <  2:178464/100‑1 (MQ=255)
                                                                   ttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTg                                 <  2:1121308/100‑1 (MQ=255)
                                                                   ttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGAAGCGACTCACTTg                                 <  2:2890922/100‑1 (MQ=255)
                                                                      tGGCGTTGTTGTTGAGTGTCTGGGTCAGCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGGGCGGTTGCGCTGGCCTGCCATGATGCGAc                                        <  1:4024231/90‑1 (MQ=255)
                                                                      tGGCGTTGTTGTTGAGTCTCGGGGTCAACGCGTTGACCACGTTGACGGCCTTCGCTGGAGTGCGGTTGCGCTTGCCTGCCATTATGCGACTCACTTGCa                               <  1:1463898/99‑1 (MQ=255)
                                                                       ggCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTAGACGGCCTTTGAGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGTATc                             <  1:437651/100‑1 (MQ=255)
                                                                         cGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGTGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                           >  2:2195390/1‑100 (MQ=255)
                                                                         cGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                           <  1:2475623/100‑1 (MQ=255)
                                                                         cGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                           >  1:1743879/1‑100 (MQ=255)
                                                                          gttgttgttAAGTCTTGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGc                          <  2:3776086/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)
                                                                                gttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTg                    >  2:3395892/1‑100 (MQ=255)
                                                                                 ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                   <  2:723356/100‑1 (MQ=255)
                                                                                 ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                   <  2:4230870/100‑1 (MQ=255)
                                                                                  tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                  >  1:3975566/1‑100 (MQ=255)
                                                                                  tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                  >  1:3698134/1‑100 (MQ=255)
                                                                                   aaGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGTGACTCACTTGCATCGGCATTCTGgcg                 <  2:2703019/100‑1 (MQ=255)
                                                                                   aaGTCTCGCGGTCAACGTGTTGACCCCTTTGACGGCCTTCGCGGGCGTGTGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                 <  1:544469/100‑1 (MQ=255)
                                                                                    gtcaacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                >  2:2788290/6‑100 (MQ=39)
                                                                                     gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCa               >  1:1559504/1‑100 (MQ=255)
                                                                                     gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCa               <  1:1404557/100‑1 (MQ=255)
                                                                                     gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCa               >  1:2281128/1‑100 (MQ=255)
                                                                                      tctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAg              <  2:1587164/100‑1 (MQ=255)
                                                                                        tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGt             >  1:2460758/1‑99 (MQ=255)
                                                                                        acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt            >  1:44769/2‑100 (MQ=255)
                                                                                        acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt            >  2:1713787/2‑100 (MQ=255)
                                                                                         cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCATTTc           >  2:128668/1‑100 (MQ=255)
                                                                                         cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc           >  2:600982/1‑100 (MQ=255)
                                                                                           gggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc         >  2:4320017/1‑100 (MQ=255)
                                                                                             gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCgcgc       <  1:1281886/100‑1 (MQ=255)
                                                                                              tCAACGTGTTGCCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCt      <  2:648385/100‑1 (MQ=255)
                                                                                               cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc     >  2:2502003/1‑100 (MQ=255)
                                                                                               cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc     >  1:1246761/1‑100 (MQ=255)
                                                                                                aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTcc    <  2:1964288/100‑1 (MQ=255)
                                                                                                 aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct   >  1:2038102/1‑100 (MQ=255)
                                                                                                  cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctg  >  1:3435461/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTG  >  NC_000913/4291413‑4291610

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: