Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,512 G→A 17.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,5120GA17.8% 346.7 / 55.2 150M1M (GTG→ATG) †nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (70/53);  new base A (12/15);  total (82/68)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

GACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTG  >  NC_000913/4291413‑4291610
                                                                                                   |                                                                                                  
   aGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAATTCGGTTTTGTGGCGTTGTTGTTAAGTCTCGGGGTCAACgtgt                                                                                                 >  1:1888560/1‑100 (MQ=255)
    gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                <  1:1297388/100‑1 (MQ=255)
    gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                >  2:1297393/1‑100 (MQ=255)
    gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                >  2:2628509/1‑100 (MQ=255)
    gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                >  1:1496592/1‑100 (MQ=255)
     ttGTCGGTGGCCTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCATTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                <  1:1094193/99‑1 (MQ=255)
      tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                              >  1:2187238/1‑100 (MQ=255)
         cGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAccc                                                                                           >  2:2593116/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:2726382/100‑1 (MQ=255)
          ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                          <  2:233737/100‑1 (MQ=255)
          ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                          <  1:5504/100‑1 (MQ=255)
          ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGGTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAcccc                                                                                          <  1:959416/100‑1 (MQ=255)
           gTGGCGTAAGCCAGAAAAGTGGTGGTTTCTGGACGTATTGCGTCCGGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                          <  1:1911407/99‑1 (MQ=255)
            tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCg                                                                                         >  1:2371037/1‑99 (MQ=255)
            tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                        <  2:1951636/100‑2 (MQ=255)
             ggCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgag                                                                                       <  2:519674/100‑3 (MQ=255)
             ggCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgag                                                                                       <  2:1826206/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)
                gTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGAc                                                                                    <  2:1072125/100‑1 (MQ=255)
                gTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagac                                                                                     <  2:615536/99‑5 (MQ=39)
                 tAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACg                                                                                   <  2:1107665/100‑1 (MQ=255)
                 tAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACg                                                                                   <  2:2665822/100‑1 (MQ=255)
                  aaGCCATAAGAGGTGTGGTTTCTGGACGTATCCCTTCTGGGAGTCGGGTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGTCCCCgagactga                                                                                  <  1:1969926/100‑8 (MQ=255)
                  aaGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                  <  2:2915717/100‑1 (MQ=255)
                  aaGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                  <  2:1387640/100‑1 (MQ=255)
                    gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGGAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagacttaa                                                                                 <  1:1376169/99‑9 (MQ=38)
                    gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGcc                                                                                <  1:828115/100‑1 (MQ=255)
                     ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCt                                                                               >  1:1826488/1‑100 (MQ=255)
                     ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCt                                                                               >  1:889928/1‑100 (MQ=255)
                     ccAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCgagacttaaca                                                                               >  2:1077958/1‑90 (MQ=37)
                      cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCtt                                                                              >  1:477694/1‑100 (MQ=255)
                        gAAAAGAGGTGGTTTCTGGACGTATTCCGGCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGc                                                                                 <  1:2416088/95‑1 (MQ=255)
                            agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                        <  1:2823547/100‑1 (MQ=255)
                            agagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                        <  2:551910/100‑1 (MQ=255)
                              agGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                      <  2:484933/100‑1 (MQ=255)
                              agGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                      <  2:161339/100‑1 (MQ=255)
                               ggtggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGt                                                                     >  1:2888219/1‑100 (MQ=255)
                                 tggtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                   >  2:1339563/1‑100 (MQ=255)
                                   gtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCgg                                                                 <  2:2466112/100‑1 (MQ=255)
                                    tttGTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGGGCCCGTTGACGGCCTTCGCGGGAGGGCGGt                                                                <  1:2504549/100‑1 (MQ=255)
                                    tttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                <  2:2259007/100‑1 (MQ=255)
                                     ttCTGGACGTATTCCTTCCTGAAGTCGGTTTTGTGGCGTAGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGTCCTTCGCGGAAGTGCGGtt                                                               >  2:2764320/1‑100 (MQ=255)
                                       cTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgc                                                             >  1:169719/1‑100 (MQ=255)
                                        tGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                            >  2:1947345/1‑100 (MQ=255)
                                         ggACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                           >  1:1452797/1‑100 (MQ=255)
                                          gACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCt                                                          <  2:649957/100‑1 (MQ=255)
                                           aCGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTg                                                         <  1:1599849/100‑1 (MQ=255)
                                            cGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTgg                                                        >  2:2154316/1‑100 (MQ=255)
                                              tATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                      >  1:810554/1‑100 (MQ=255)
                                               aTTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                      >  2:284515/1‑99 (MQ=255)
                                                  ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                  <  1:1898937/100‑1 (MQ=255)
                                                  ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                  <  1:1403302/100‑1 (MQ=255)
                                                  ccttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                  >  2:504798/1‑100 (MQ=255)
                                                   cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGCCGGCCTTCGCGGGAGCGCGGTTGCGCTGGACTGcc                                                  >  1:2426764/1‑99 (MQ=255)
                                                     tcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatg                                               <  2:889931/100‑1 (MQ=255)
                                                        tGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCctgat                                             >  1:1936453/1‑99 (MQ=255)
                                                        tGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                            >  2:574256/1‑100 (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)
                                                                         cGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTTACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                           <  1:914178/100‑1 (MQ=255)
                                                                               tgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCt                     >  2:1908460/1‑100 (MQ=255)
                                                                               tgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCt                     >  1:342131/1‑100 (MQ=255)
                                                                                gttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATcagatcggaa                    >  2:1635577/1‑91 (MQ=255)
                                                                                 ttAAGTCTCGGGGTCAACGTGTTGGCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                   <  1:1984276/100‑1 (MQ=255)
                                                                                  tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGTTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgc                  >  2:529630/1‑100 (MQ=255)
                                                                                   ctgTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                 <  2:2218261/98‑1 (MQ=255)
                                                                                   aaGTCTCGGGGCCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                 >  1:2029527/1‑100 (MQ=255)
                                                                                    aGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                >  2:320146/1‑100 (MQ=255)
                                                                                     gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgag               >  1:2218248/1‑98 (MQ=255)
                                                                                       ctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGt             >  1:2795355/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:2355022/2‑100 (MQ=255)
                                                                                        acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt            >  2:2306165/2‑100 (MQ=255)
                                                                                        acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGGGCGGGTGCGCTGGCCTGCAATGATCCGACCCACTTGCATCGGCATTCTGgc                  >  1:1815095/2‑94 (MQ=255)
                                                                                         cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCCCTTGCATCGGCATTCTGGCGCAGTTc           >  2:2376830/1‑100 (MQ=255)
                                                                                         cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc           >  1:2134781/1‑100 (MQ=255)
                                                                                          ggggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcg          <  2:1846757/100‑1 (MQ=255)
                                                                                          ggggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcg          >  1:235941/1‑100 (MQ=255)
                                                                                           gggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc         >  2:1253756/1‑100 (MQ=255)
                                                                                           gggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc         >  1:2765730/1‑100 (MQ=255)
                                                                                            ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgcg        <  1:945980/100‑1 (MQ=255)
                                                                                            ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgcg        >  2:1388487/1‑100 (MQ=255)
                                                                                             gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCgcgc       <  2:1013187/100‑1 (MQ=255)
                                                                                               cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc     >  1:1192717/1‑100 (MQ=255)
                                                                                               cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc     <  2:1331921/100‑1 (MQ=255)
                                                                                               cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTc     >  1:1328288/1‑100 (MQ=255)
                                                                                                  cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctg  >  1:312920/1‑100 (MQ=255)
                                                                                                  cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctg  >  1:58370/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: