Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,532,537 G→C 15.1% A216G (GCC→GGC)  yjhP ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,5370GC15.1% 422.1 / 61.1 153A216G (GCC→GGC) yjhPputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (51/77);  new base C (14/9);  total (66/87)
Fisher's exact test for biased strand distribution p-value = 7.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 169 total aligned reads displayed.

TACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTT  >  NC_000913/4532438‑4532636
                                                                                                   |                                                                                                   
tACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggcgg                                                                                                     >  1:1614089/1‑100 (MQ=255)
 aCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCggcggc                                                                                                    >  1:916806/1‑100 (MQ=255)
  ccGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGgcggcg                                                                                                   >  2:1395881/1‑100 (MQ=255)
  ccGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGgcggcg                                                                                                   >  1:1694142/1‑100 (MQ=255)
  ccGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGgcggcg                                                                                                   <  1:3477637/100‑1 (MQ=255)
    gggCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCggcggc                                                                                                    >  2:447816/1‑97 (MQ=255)
    gggCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGaa                                                                                                 >  1:877297/1‑100 (MQ=255)
      gCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGt                                                                                               <  1:3073915/100‑1 (MQ=255)
      gCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGt                                                                                               <  2:332852/100‑1 (MQ=255)
      gCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTCTCGGCGCGATCTTTAGTTCGGCCCTGACGTCGGCGGCGAAGt                                                                                               <  1:2463806/100‑1 (MQ=255)
         aTGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCag                                                                                            <  1:2797808/100‑2 (MQ=255)
           gAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCg                                                                                          <  1:413458/100‑1 (MQ=255)
           gAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCaggg                                                                                          <  2:1925164/100‑4 (MQ=255)
           gAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCaggg                                                                                          <  2:2955239/100‑4 (MQ=255)
           gAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCaggg                                                                                          <  2:2340315/100‑4 (MQ=255)
            aGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                         >  1:2823043/1‑100 (MQ=255)
            aGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                         <  2:877290/100‑1 (MQ=255)
             gTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTc                                                                                        >  1:15753/1‑100 (MQ=255)
               gCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAg                                                                                      <  1:2977423/100‑1 (MQ=255)
               gCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAg                                                                                      >  2:1083360/1‑100 (MQ=255)
                   acacGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGAtt                                                                                  >  2:2688449/1‑100 (MQ=255)
                      cGCCCCGGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGATTCGGCGCGCTGTTTAGTACGCCCCTGACTTCTGCGACGAAGTCAACGTCACGCTTCGc                                                                               >  2:1254192/1‑100 (MQ=255)
                      cGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGc                                                                               >  2:1537475/1‑100 (MQ=255)
                        ccccAGCCAAAGCATTCCCGCGCGTAGGTCACGTCGCGTTTCGGCGCGAGGTTTAGTTCGGCCCTGACTTCCGCGGCGACGCCCTCTTCAGGGTTCGCtt                                                                             >  2:2578702/1‑100 (MQ=255)
                        ccccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCtt                                                                             >  2:1814674/1‑100 (MQ=255)
                         cccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                            >  1:2501995/1‑100 (MQ=255)
                         cccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                            <  2:875053/100‑1 (MQ=255)
                          ccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                           >  2:2209872/1‑100 (MQ=255)
                           cAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcca                                                                          >  2:2241974/1‑100 (MQ=255)
                           cAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCCGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcca                                                                          <  2:738226/100‑1 (MQ=255)
                            aGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccag                                                                         <  1:461065/100‑1 (MQ=255)
                             gCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTCTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGTTTccagc                                                                        <  2:1912828/100‑1 (MQ=255)
                              ccAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagcc                                                                       >  1:2239486/1‑100 (MQ=255)
                              ccAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagcc                                                                       >  1:3276646/1‑100 (MQ=255)
                              atctaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagcc                                                                       <  2:3220174/96‑1 (MQ=255)
                               cAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagcca                                                                      <  1:1409199/100‑1 (MQ=255)
                               cAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagcca                                                                      >  2:2121869/1‑100 (MQ=255)
                                aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagccag                                                                     >  1:2927529/1‑100 (MQ=255)
                                aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagccag                                                                     >  1:3219990/1‑100 (MQ=255)
                                  aGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcg                                                                   >  1:1054933/1‑100 (MQ=255)
                                   gCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGAGGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCcagccagc                                                                    <  2:688727/98‑1 (MQ=255)
                                     aTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcggcg                                                                <  1:1674663/100‑1 (MQ=255)
                                     aTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGGTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcggcg                                                                <  2:1780663/100‑1 (MQ=255)
                                       tCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGCCCCTCACTTCGGCGGCGAAGTCATCGTCACGATTCGCTTCCAGCCAGCCGCGCa                                                              >  2:3362923/1‑100 (MQ=255)
                                        cccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAt                                                             >  2:2571461/1‑100 (MQ=255)
                                        cccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAt                                                             >  2:3623827/1‑100 (MQ=255)
                                        cccGCGCGTAGGTCACGTACCGGTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAt                                                             <  2:2823023/100‑1 (MQ=255)
                                         ccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATg                                                            >  1:2817058/1‑100 (MQ=255)
                                              cgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAg                                                       >  1:992393/1‑100 (MQ=255)
                                               gTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGc                                                      <  1:2060985/100‑1 (MQ=255)
                                                tAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGcc                                                     <  1:2659406/100‑1 (MQ=255)
                                                 aGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCa                                                    <  1:296498/100‑1 (MQ=255)
                                                  ggTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCAGCGCATGGTCAGCCAt                                                   >  2:2133247/1‑100 (MQ=255)
                                                   gTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAtt                                                  <  2:1694134/100‑1 (MQ=255)
                                                    tCACGTACCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAttt                                                 >  1:1907097/1‑100 (MQ=255)
                                                       cGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                              <  2:1097615/100‑1 (MQ=255)
                                                        gTAGCGTTTCGGCGCGATGTTTAGTTGGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGTGGCGCATGGTCAGCCATTTTGCg                                             <  2:978414/100‑1 (MQ=255)
                                                        gTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCg                                             <  1:1081584/100‑1 (MQ=255)
                                                        gTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCg                                             <  1:2590565/100‑1 (MQ=255)
                                                        gTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCg                                             <  1:1500060/100‑1 (MQ=255)
                                                          aGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGc                                           <  1:2124441/100‑1 (MQ=255)
                                                            cGTTTCGGCGCGATCTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                         <  2:2480187/100‑1 (MQ=255)
                                                             gTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTc                                        <  2:2471875/100‑1 (MQ=255)
                                                             gTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTc                                        >  1:2055679/1‑100 (MQ=255)
                                                             gTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTc                                        <  1:718722/100‑1 (MQ=255)
                                                             gTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTc                                        >  1:1771463/1‑100 (MQ=255)
                                                             gTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCAGCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTc                                        <  2:1170802/100‑1 (MQ=255)
                                                              tttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTc                                        <  2:2532660/99‑1 (MQ=255)
                                                              tttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCa                                       <  2:808974/100‑1 (MQ=255)
                                                              tttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCa                                       <  2:91237/100‑1 (MQ=255)
                                                              tttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCa                                       <  2:1483824/100‑1 (MQ=255)
                                                               ttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCAt                                      >  1:2247650/1‑100 (MQ=255)
                                                                  ggCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAcc                                   >  1:2864769/1‑100 (MQ=255)
                                                                    cgcgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCgg                                 >  1:250255/1‑100 (MQ=255)
                                                                          gTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCa                            >  2:1561005/1‑99 (MQ=255)
                                                                          gTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAg                           <  1:1339419/100‑1 (MQ=255)
                                                                          gTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAg                           >  1:1237144/1‑100 (MQ=255)
                                                                          gTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAg                           <  2:1222271/100‑1 (MQ=255)
                                                                           tttAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGc                          <  1:2391058/100‑1 (MQ=255)
                                                                                 ctgacgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                    >  2:2143595/8‑100 (MQ=38)
                                                                                  cGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                   <  2:2511848/100‑1 (MQ=255)
                                                                                  cGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                   <  2:462267/100‑1 (MQ=255)
                                                                                  cGGCCCTGACTTCGGCGGCAAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                   <  1:3275182/100‑1 (MQ=255)
                                                                                   ggCCTTGACTTCGGCGTCGAAGTCATCGTCAGGATTCGCTTCCAGCCAACGGCGCATAGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                   <  2:3410463/99‑1 (MQ=255)
                                                                                   gacgatGACTTCGCCACCGAAGTCATCGTCAGGATTAGCTTCCCGCCAGCGGCGCATGGTCAGCCATTGTGCGGCTTCATACCGGTCCCAGCCTTCCt                    >  2:591474/6‑98 (MQ=255)
                                                                                     cgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCCTGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                >  2:2961113/4‑100 (MQ=255)
                                                                                     cgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGt                 >  1:2168744/4‑99 (MQ=255)
                                                                                     cgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                >  1:348575/4‑100 (MQ=255)
                                                                                     cgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                >  1:1678220/4‑100 (MQ=255)
                                                                                     cccTGACTTCGGCGGCGAAGGCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                >  2:1740839/1‑100 (MQ=255)
                                                                                      gatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt               >  1:1493337/3‑100 (MQ=255)
                                                                                      gatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt               >  1:790048/3‑100 (MQ=255)
                                                                                      ccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt               <  1:340022/100‑1 (MQ=255)
                                                                                      ccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt               >  1:1217134/1‑100 (MQ=255)
                                                                                       cTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTg              >  1:1183097/1‑100 (MQ=255)
                                                                                         gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGcc            <  2:2760296/100‑1 (MQ=255)
                                                                                         gACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGcc            <  1:2143609/100‑1 (MQ=255)
                                                                                            ttCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGt         <  2:2063392/100‑1 (MQ=255)
                                                                                               ggcggcGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTAcc      <  1:957779/100‑1 (MQ=255)
                                                                                                gcggcgAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCa     <  2:1798685/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTT  >  NC_000913/4532438‑4532636

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: