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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,5340GC18.1% 431.2 / 89.6 171A217G (GCC→GGC) yjhPputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (69/71);  new base C (12/19);  total (81/90)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

Only 100 of 182 total aligned reads displayed.

TGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACC  >  NC_000913/4532435‑4532632
                                                                                                   |                                                                                                  
tgttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcgg                                                                                                    <  2:1475008/100‑1 (MQ=255)
 gttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACTTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggc                                                                                                   <  1:3503603/100‑1 (MQ=255)
   tACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggcgg                                                                                                 <  2:1326439/100‑1 (MQ=255)
   tACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggcgg                                                                                                 <  2:1880048/100‑1 (MQ=255)
   tACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggcgg                                                                                                 >  1:1880051/1‑100 (MQ=255)
        ggCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAg                                                                                            <  2:3362781/100‑1 (MQ=255)
          cGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTc                                                                                          >  2:2837746/1‑100 (MQ=255)
            aTGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCAt                                                                                        >  2:1826416/1‑100 (MQ=255)
            aTGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCag                                                                                        <  1:3994735/100‑2 (MQ=255)
            aTGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCag                                                                                        <  2:3788815/100‑2 (MQ=255)
             tGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagg                                                                                       <  2:2269682/100‑3 (MQ=255)
             tGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagg                                                                                       <  2:1599309/100‑3 (MQ=255)
             tGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagg                                                                                       <  2:3796950/100‑3 (MQ=255)
             tGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagg                                                                                       <  2:4143882/100‑3 (MQ=255)
             tGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagg                                                                                       <  2:757455/100‑3 (MQ=255)
              gAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTTACTTCGCCGCCGAAGTCagg                                                                                       <  2:2563394/99‑3 (MQ=255)
              gAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCaggg                                                                                      <  1:1013657/100‑4 (MQ=255)
              gAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCaggg                                                                                      <  1:1836879/100‑4 (MQ=255)
               aGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                     >  1:3381854/1‑100 (MQ=255)
               aGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagggc                                                                                     <  1:3708440/100‑5 (MQ=255)
               aGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagggc                                                                                     <  1:1745606/100‑5 (MQ=255)
               aGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagggc                                                                                     <  1:2923970/100‑5 (MQ=255)
                gTGCGAACACGCCCCAGCGAAAGCATGCCCGGGCGTACGTCACGTAGCGTTTCGGCGCGATGTCTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagggc                                                                                     <  2:2551781/99‑5 (MQ=255)
                gTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagggcc                                                                                    <  1:3491335/100‑6 (MQ=255)
                 tGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCa                                                                                   >  2:1200681/1‑100 (MQ=255)
                  gCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAg                                                                                  >  1:4036730/1‑100 (MQ=255)
                  gCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAg                                                                                  >  2:928083/1‑100 (MQ=255)
                   cGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAgg                                                                                 >  2:3914791/1‑100 (MQ=255)
                   cGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAgg                                                                                 <  2:784490/100‑1 (MQ=255)
                    gAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGa                                                                                >  1:495593/1‑100 (MQ=255)
                       cacGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTc                                                                             <  1:3451715/100‑1 (MQ=255)
                        acGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCg                                                                            <  1:3942131/100‑1 (MQ=255)
                          gCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCt                                                                          <  1:480422/100‑1 (MQ=255)
                          gCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCt                                                                          >  1:2801682/1‑100 (MQ=255)
                          gCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCt                                                                          <  1:3567306/100‑1 (MQ=255)
                          gCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCt                                                                          <  2:246247/100‑1 (MQ=255)
                           ccccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCtt                                                                         >  1:2045251/1‑100 (MQ=255)
                           ccccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCtt                                                                         >  1:792764/1‑100 (MQ=255)
                           ccccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCtt                                                                         >  2:2337147/1‑100 (MQ=255)
                           ccccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCtt                                                                         >  2:4002403/1‑100 (MQ=255)
                            cccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                        <  1:4165546/100‑1 (MQ=255)
                             ccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTGCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcc                                                                       <  2:1497853/100‑1 (MQ=255)
                               aGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccag                                                                     <  1:4123834/100‑1 (MQ=255)
                                  cAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagcca                                                                  >  1:2332599/1‑100 (MQ=255)
                                  cAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagcca                                                                  >  1:1568898/1‑100 (MQ=255)
                                  cAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagcca                                                                  >  1:2712834/1‑100 (MQ=255)
                                   aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTtcagc                                                                    >  2:2001987/1‑97 (MQ=255)
                                    aaGCATTCCCGCGCGTAGGTCACGTAGCGTTCCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCCGCGAAGTCATCGTCAGGATTCGCTTCcagccagc                                                                <  2:4183314/100‑1 (MQ=255)
                                      gCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcgg                                                              <  1:1970120/100‑1 (MQ=255)
                                       cgttCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGCTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcggc                                                             >  1:2850430/3‑100 (MQ=255)
                                        aTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAgcggcg                                                            <  1:1405199/100‑1 (MQ=255)
                                         ttCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGgcgc                                                           <  2:634898/100‑1 (MQ=255)
                                           cccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAt                                                         >  1:2354512/1‑100 (MQ=255)
                                           cccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAt                                                         >  1:3171269/1‑100 (MQ=255)
                                            ccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATg                                                        >  1:3040696/1‑100 (MQ=255)
                                             cgcgcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATgg                                                       <  1:189506/100‑1 (MQ=255)
                                              gcgcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGt                                                      >  2:3175249/1‑100 (MQ=255)
                                              gcgcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGt                                                      <  2:704963/100‑1 (MQ=255)
                                               cgcgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTc                                                     <  1:866669/100‑1 (MQ=255)
                                                 cgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAg                                                   <  2:3034358/100‑1 (MQ=255)
                                                 cgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAg                                                   <  1:1293512/100‑1 (MQ=255)
                                                 cgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAg                                                   >  1:2608949/1‑100 (MQ=255)
                                                 cgTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGACGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAg                                                   >  1:3796741/1‑100 (MQ=255)
                                                  gTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGc                                                  <  1:2799458/100‑1 (MQ=255)
                                                    aGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCa                                                <  1:671284/100‑1 (MQ=255)
                                                    aGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCa                                                >  1:1589383/1‑100 (MQ=255)
                                                     ggTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAt                                               <  2:3713202/100‑1 (MQ=255)
                                                     ggTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAt                                               >  2:2099993/1‑100 (MQ=255)
                                                       tCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAttt                                             >  2:3846891/1‑100 (MQ=255)
                                                           gTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                          <  1:323615/99‑1 (MQ=255)
                                                           gTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCg                                         <  1:251759/100‑1 (MQ=255)
                                                             aGCGTCTCGGCGCGATGTTGAGTTCGGCACTGACTCCGGCGGCGAAGTCGCCGTCAGCATTCGCTGCCAGCCAGCGGCGCATGGCCAGCCATTTTGCGGc                                       <  2:3030655/100‑1 (MQ=255)
                                                              gCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCt                                      >  2:677171/1‑100 (MQ=255)
                                                                gTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTc                                    >  1:3846465/1‑100 (MQ=255)
                                                                 tttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCa                                   <  2:3391505/100‑1 (MQ=255)
                                                                       cgcgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCgg                             <  2:3312795/100‑1 (MQ=255)
                                                                            tGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCa                        >  1:1192957/1‑100 (MQ=255)
                                                                            tGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCa                        <  2:3550592/100‑1 (MQ=255)
                                                                              tttAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGc                      <  2:564582/100‑1 (MQ=255)
                                                                                 aGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCtt                   <  2:4007060/100‑1 (MQ=255)
                                                                                 aGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGCCCCAGCCtt                   >  2:3284674/1‑100 (MQ=255)
                                                                                  gTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTCCATACCGGTCCCAGCCTTc                  >  2:3158853/1‑100 (MQ=255)
                                                                                   ttCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTcc                 <  2:3895593/100‑1 (MQ=255)
                                                                                   ttCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTcc                 <  1:1761199/100‑1 (MQ=255)
                                                                                    ctgacgatGACTTCGCCGCCGAAGTCATCGTCCGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                >  2:1760995/8‑100 (MQ=38)
                                                                                       gCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGt             >  2:2032926/1‑100 (MQ=255)
                                                                                        cgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc            >  2:3638254/4‑100 (MQ=255)
                                                                                        cgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc            >  1:2335478/4‑100 (MQ=255)
                                                                                        cccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc            >  1:1705857/1‑100 (MQ=255)
                                                                                        cccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc            <  1:3652393/100‑1 (MQ=255)
                                                                                         gatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt           >  1:407169/3‑100 (MQ=255)
                                                                                         ccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt           <  2:2190886/100‑1 (MQ=255)
                                                                                          cTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTg          <  2:1265877/100‑1 (MQ=255)
                                                                                          cTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTg          >  2:623913/1‑100 (MQ=255)
                                                                                            gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGcc        <  2:1683961/100‑1 (MQ=255)
                                                                                            gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGcc        >  1:2181250/1‑100 (MQ=255)
                                                                                             aCTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCa       <  1:534974/100‑1 (MQ=255)
                                                                                                 cgggggCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTAc   <  1:1038944/100‑1 (MQ=255)
                                                                                                 cggcggCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTa    <  2:1152695/99‑1 (MQ=255)
                                                                                                 cggcggCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTAc   >  2:2372323/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACC  >  NC_000913/4532435‑4532632

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: