Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,513 T→C 18.2% 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,5130TC18.2% 310.5 / 64.9 150M1A (GTG→GCG) †nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base T (69/53);  new base C (12/15);  total (82/68)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

Only 100 of 152 total aligned reads displayed.

AGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTG  >  NC_000913/4291416‑4291610
                                                                                                 |                                                                                                 
aGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAATTCGGTTTTGTGGCGTTGTTGTTAAGTCTCGGGGTCAACgtgt                                                                                                 >  1:1888560/1‑100 (MQ=255)
 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                >  2:1297393/1‑100 (MQ=255)
 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                >  2:2628162/1‑100 (MQ=255)
 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                >  2:2628509/1‑100 (MQ=255)
 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                <  1:1297388/100‑1 (MQ=255)
 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                >  1:1496592/1‑100 (MQ=255)
  ttGTCGGTGGCCTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCATTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                <  1:1094193/99‑1 (MQ=255)
  ttGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTg                                                                                               <  2:2758544/100‑1 (MQ=255)
   tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                              >  1:2187238/1‑100 (MQ=255)
      cGGTGGCGTAAGCCAGAAAAGAGAAGGTTTCTGGACGTATTGCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAccc                                                                                           >  1:1081914/1‑100 (MQ=255)
       ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                          <  1:222233/100‑1 (MQ=255)
       ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                          <  2:2891987/100‑1 (MQ=255)
       ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                          <  1:1889732/100‑1 (MQ=255)
        gTGGCGTAAGCCAGAAAAGTGGTGGTTTCTGGACGTATTGCGTCCGGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                          <  1:1911407/99‑1 (MQ=255)
        gTGGCGTAAGCCAGAAAAGCGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCg                                                                                         <  2:1308753/100‑1 (MQ=255)
         tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCg                                                                                         >  1:2371037/1‑99 (MQ=255)
         tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                        <  2:1951636/100‑2 (MQ=255)
          ggCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgag                                                                                       <  2:519674/100‑3 (MQ=255)
           gCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTg                                                                                      <  2:740740/100‑1 (MQ=255)
            cGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCgagac                                                                                     <  2:2187251/100‑5 (MQ=255)
            cGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGa                                                                                     <  2:85648/100‑1 (MQ=255)
             gTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagac                                                                                     <  2:615536/99‑5 (MQ=39)
              tAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACg                                                                                   <  2:2266882/100‑1 (MQ=255)
               aaGCCATAAGAGGTGTGGTTTCTGGACGTATCCCTTCTGGGAGTCGGGTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGTCCCCgagactga                                                                                  <  1:1969926/100‑8 (MQ=255)
               aaGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                  <  2:2915717/100‑1 (MQ=255)
                aGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGc                                                                                 >  1:319404/1‑100 (MQ=255)
                 gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGGAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagacttaa                                                                                 <  1:1376169/99‑9 (MQ=38)
                 gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGcc                                                                                <  1:828115/100‑1 (MQ=255)
                  ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCt                                                                               >  1:1826488/1‑100 (MQ=255)
                  ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagacttaaca                                                                               >  2:1077958/1‑90 (MQ=37)
                   cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCtt                                                                              >  1:477694/1‑100 (MQ=255)
                     gAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTcg                                                                            <  2:1200617/100‑1 (MQ=255)
                     gAAAAGAGGTGGTTTCTGGACGTATTCCTTCCCGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTcg                                                                            >  2:2427284/1‑100 (MQ=255)
                     gAAAAGAGGTGGTTTCTGGACGTATTCCGGCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGc                                                                                 <  1:2416088/95‑1 (MQ=255)
                        aaGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                         >  1:1978015/1‑100 (MQ=255)
                         agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                        <  2:551910/100‑1 (MQ=255)
                         agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                        <  1:2823547/100‑1 (MQ=255)
                           agGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                      <  2:484933/100‑1 (MQ=255)
                            ggtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                     >  1:2101440/1‑100 (MQ=255)
                              tggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                   >  2:1339563/1‑100 (MQ=255)
                                gtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCgg                                                                 <  2:2466112/100‑1 (MQ=255)
                                ctttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCg                                                                  >  1:565974/2‑99 (MQ=255)
                                 tttGTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGGGCCCGTTGACGGCCTTCGCGGGAGGGCGGt                                                                <  1:2504549/100‑1 (MQ=255)
                                  ttCTGGACGTATTCCTTCCTGAAGTCGGTTTTGTGGCGTAGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGTCCTTCGCGGAAGTGCGGtt                                                               >  2:2764320/1‑100 (MQ=255)
                                   tCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTg                                                              <  2:134288/100‑1 (MQ=255)
                                    cTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgc                                                             >  1:169719/1‑100 (MQ=255)
                                     tGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                            >  2:1947345/1‑100 (MQ=255)
                                      ggACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                           >  1:1452797/1‑100 (MQ=255)
                                       gACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCt                                                          <  1:839323/100‑1 (MQ=255)
                                       gACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCt                                                          <  2:649957/100‑1 (MQ=255)
                                        aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                         <  1:1599849/100‑1 (MQ=255)
                                         cGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTgg                                                        >  2:2154316/1‑100 (MQ=255)
                                         cGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGGTGCTAAGACTCGGGGTCAGCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTgg                                                        <  1:2836449/100‑1 (MQ=255)
                                            aTTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                      >  2:284515/1‑99 (MQ=255)
                                            aTTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGc                                                       >  1:2659310/1‑98 (MQ=255)
                                               ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                  >  2:504798/1‑100 (MQ=255)
                                               ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                  <  1:1403302/100‑1 (MQ=255)
                                                cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGCCGGCCTTCGCGGGAGCGCGGTTGCGCTGGACTGcc                                                  >  1:2426764/1‑99 (MQ=255)
                                                  tcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatg                                               <  2:889931/100‑1 (MQ=255)
                                                     tGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                            >  2:574256/1‑100 (MQ=255)
                                                     tGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                            <  2:1673723/100‑1 (MQ=255)
                                                        aggcGGGTTTCTGGTGACGTTGTTAAGTCTCGGGGGCAACGACTTGGCCCCGTTGACGGCCTTCGAGGGAGTTCGGTTGCGCTGACCTGCCatg                                               >  1:1533840/4‑94 (MQ=255)
                                                         gTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                        <  2:214335/100‑1 (MQ=255)
                                                            ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                     >  1:1730166/1‑100 (MQ=255)
                                                            ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCa                                     >  1:2342917/1‑100 (MQ=255)
                                                                ttCTGGCGTTGTTGGTAAGTCTCGGGGTCCACGTGTTGCCCCCGTTGACGGCCTTCGCGGGCGTGCGGTTGCGCTGGCCTGCCATGATGCGCCTCACTTg                                 <  1:616131/100‑1 (MQ=255)
                                                                   tGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                              >  1:2374638/1‑100 (MQ=255)
                                                                    ggCGTTGTTGTTAAGTCTCGGGGTCACCGTGTTGACCCAGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGACATGATGCGACTCACTTGCCTc                             >  1:1237380/1‑100 (MQ=255)
                                                                      cGTTGTTGTTAAGTCTCGTGGTCAACGTGTTGACCCCGTTGACGTCCTTCGCGGGATTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCg                            <  1:2791745/99‑1 (MQ=255)
                                                                      cGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTTACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                           <  1:914178/100‑1 (MQ=255)
                                                                         tgttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCAt                        <  1:958358/100‑1 (MQ=255)
                                                                            tgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCt                     >  2:1908460/1‑100 (MQ=255)
                                                                             gttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATcagatcggaa                    >  2:1635577/1‑91 (MQ=255)
                                                                               tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGTTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                  >  2:529630/1‑100 (MQ=255)
                                                                                ctgTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                 <  2:2218261/98‑1 (MQ=255)
                                                                                aaGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                 <  2:1081111/100‑1 (MQ=255)
                                                                                aaGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                 <  1:710202/100‑1 (MQ=255)
                                                                                aaGTCTCGGGGCCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                 >  1:2029527/1‑100 (MQ=255)
                                                                                 aGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                >  2:320146/1‑100 (MQ=255)
                                                                                  gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgag               >  1:2218248/1‑98 (MQ=255)
                                                                                  gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCa               >  1:1926312/1‑100 (MQ=255)
                                                                                    ctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGt             >  2:1168656/1‑100 (MQ=255)
                                                                                    aacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGt             >  1:555094/3‑100 (MQ=255)
                                                                                    aacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGt             >  2:1996408/3‑100 (MQ=255)
                                                                                     acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt            >  2:2306165/2‑100 (MQ=255)
                                                                                     acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt            >  2:2355022/2‑100 (MQ=255)
                                                                                     acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGGGCGGGTGCGCTGGCCTGCAATGATCCGACCCACTTGCATCGGCATTCTGgc                  >  1:1815095/2‑94 (MQ=255)
                                                                                      cGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc           >  2:1041788/1‑100 (MQ=255)
                                                                                       ggggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcg          <  2:1846757/100‑1 (MQ=255)
                                                                                       ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcg          >  1:235941/1‑100 (MQ=255)
                                                                                        gggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc         >  1:2765730/1‑100 (MQ=255)
                                                                                         ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgcg        >  2:1388487/1‑100 (MQ=255)
                                                                                         ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgcg        <  2:2353752/100‑1 (MQ=255)
                                                                                         ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgcg        <  1:945980/100‑1 (MQ=255)
                                                                                          gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCgcgc       <  2:1013187/100‑1 (MQ=255)
                                                                                            cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc     >  1:1328288/1‑100 (MQ=255)
                                                                                            cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc     >  1:1192717/1‑100 (MQ=255)
                                                                                             aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTcc    >  2:1590369/1‑100 (MQ=255)
                                                                                               cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctg  >  1:312920/1‑100 (MQ=255)
                                                                                               cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctg  <  1:504795/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
AGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTG  >  NC_000913/4291416‑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: