Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,532,537 G→C 5.3% 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,5370GC 5.3% 186.5 / 2.9 75A216G (GCC→GGC) yjhPputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (30/40);  new base C (2/2);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

AAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGG  >  NC_000913/4532389‑4532682
                                                                                                                                                    |                                                                                                                                                 
aaagGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTTGGCCCTGACTTcggcgg                                                                                                                                                   <  1:546228/149‑1 (MQ=255)
aaagGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGATCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggcgg                                                                                                                                                   >  1:730127/1‑149 (MQ=255)
    gTAATGCCTCGATATTCCCAGCCAGAATTCCCGGGAACCTGTTGTTACCGGGCGATGAGTGCGTACACGCCACATCCAAAGCATTCCCGCGCGGAGGTCACGTAGCGCTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGaa                                                                                                                                               <  1:570484/149‑1 (MQ=255)
    gTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGaa                                                                                                                                               >  1:273396/1‑149 (MQ=255)
     tAAAGCCTCGATATTTCCAGCCAGAAATACCGGTAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAg                                                                                                                                              <  2:142858/149‑1 (MQ=255)
         gCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCTGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCag                                                                                                                                          <  2:324267/149‑2 (MQ=255)
         gCCGCGATATTTCCAGCCAGAAATACCGGGAACATGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTTGGTCACGAAGCGTTGCGGCGCGATGTTTAGTTCGGCCCTTACGACGGCGGCGGAGGCaa                                                                                                                                          >  2:690119/1‑148 (MQ=255)
            tCGATATATCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                       <  1:126043/149‑1 (MQ=255)
            gcGATATTTCCCCACTGAAATCCCCGGCACCGGTCGTCACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGAAGCGATGCGGCGCGATGTTTAGTGCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                       <  2:584255/148‑1 (MQ=255)
               ataaTTCCAGCCAGAAATACTGGGAACCTGTTGTTACCGGGCGATGAGGGCCCACACGCCCCAGCCAAAGCTTTCCCGCGCGTAGGGCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAg                                                                                                                                    <  2:306709/149‑1 (MQ=255)
                   ttCCAGCGAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCCGCGGCGAAGTCATCGTCAggcta                                                                                                                                >  2:492954/1‑146 (MQ=255)
                   ttCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGAtt                                                                                                                                >  2:894884/1‑149 (MQ=255)
                    tCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTc                                                                                                                               <  2:335380/149‑1 (MQ=255)
                      cagccagAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGc                                                                                                                             <  2:334334/149‑1 (MQ=255)
                      aagccagAAATACCGAGATCCTGTAGATACCGGGCGAGGAGTGCGAACACGCCCCAGCCAAAGCATTCCCTCGCGTAGGTCAGGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGTCGGCTAAGTCATCGTCAGGATTCGc                                                                                                                             <  1:317433/148‑1 (MQ=255)
                        gacagAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCGAAATCATTCCCGCGCGTAGGTCACGTATCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCtt                                                                                                                           <  2:567796/147‑1 (MQ=255)
                         ccagAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGCACACGCCCCAGCCCAAGCATACCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                                                                          >  2:590724/1‑149 (MQ=255)
                           agAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcca                                                                                                                        >  2:409562/1‑149 (MQ=255)
                           agAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTcca                                                                                                                        <  1:771326/149‑1 (MQ=255)
                           agAAATACCGGGAACCTGTTGTCACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagg                                                                                                                                         <  2:296249/132‑3 (MQ=255)
                           agAAATACCGGGAACCTGTTGTCACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCagg                                                                                                                                         >  1:296249/1‑130 (MQ=255)
                             aaaTACCGGGACCCTGTTGTTCACGTGAGCTGAGCGCGAACAAGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagc                                                                                                                      <  1:731145/149‑1 (MQ=255)
                              aaTACCGGGAACCTGTTGTTACCGGGCGATTAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagcc                                                                                                                     <  1:987604/149‑1 (MQ=255)
                                       aaCCTGTTGTTACCGGGCGATGAGTGAGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCa                                                                                                            <  2:143620/149‑1 (MQ=255)
                                            gttgttACCGGGCGATGAGTGCGCACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGTCCCTGACATCGGCGGCGAAGTCATAGTCAGGATTCGGATCCAGCCAGCCGCGCATgtta                                                                                                       >  2:866794/1‑146 (MQ=255)
                                             ttgttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCACTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCa                                                                                                      >  2:1007597/1‑149 (MQ=255)
                                                ttACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGGAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGcc                                                                                                   <  1:469272/149‑1 (MQ=255)
                                                 tACCGGGCGATGAGTGCGAACACGCCACAGCCAAAGCATTAACGCGCGTAGGTCACGTATAGTTTAGGCGAGAAGTTAAGTTCGGACCTGACTTCGGCGGCGAAGTCTTAGTCAGGTTTCTCTTCCTGCCAGCGGCGCGTGGTAAGCCa                                                                                                  >  2:663897/1‑149 (MQ=255)
                                                  aCCGGGCGATGAGTGCGAACACGCCCCCGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGGCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCTGCCag                                                                                                 >  1:585302/1‑148 (MQ=255)
                                                  aCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTACCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACATCGTCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCACCGGCGCATTGTCAGCCAt                                                                                                 >  2:782377/1‑149 (MQ=255)
                                                  aCCGGGCGATGAGTGCGAAAACGCCCCACACAAAGCATTCCCGCGCGTAGGTCAAGTAGAGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGACATCGTCAGGAATCGCATCCAGCCAGCGGCGCATGGTCAGCCaa                                                                                                 >  2:479098/1‑148 (MQ=255)
                                                    cGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAttt                                                                                               >  2:20328/1‑149 (MQ=255)
                                                      ggCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTg                                                                                             <  2:314328/149‑1 (MQ=255)
                                                       gctttGAGTGCTAACTCGCCCCAGCCAAAGCATTCCGGCGCGGAGGGCACGTAGGGTTTCGGCGCGTTGTTTAGTTCGGCCCTGACTTCGGGGGCGTAGTGATCGTCAGGATTCGCTTCCAGCCATCGGCGCATGGTCTGCCATTTTGc                                                                                            <  1:782377/145‑1 (MQ=255)
                                                       gCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                                                                            <  1:974655/149‑1 (MQ=255)
                                                             aGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGGCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTc                                                                                      <  1:490465/149‑1 (MQ=255)
                                                              gTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCa                                                                                     >  2:546284/1‑149 (MQ=255)
                                                              gTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCa                                                                                     <  1:890006/149‑1 (MQ=255)
                                                                  gAACACGCCCCAGCCTAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGATGTCATCGGGAGGATTCGCTTCCCGCCAGCGGCGCATGGTCAGCAATTTTGCGGGTGCATTcc                                                                                 >  2:280674/1‑149 (MQ=255)
                                                                  gAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAcc                                                                                 <  1:660078/149‑1 (MQ=255)
                                                                    acaaGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGGTTCGGCGCGATGTTTAGTTAGGCCCTGACTTCGGCGGCGAAGTCATCGGCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCgg                                                                               <  2:734571/149‑1 (MQ=255)
                                                                      acGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTc                                                                             <  2:553905/149‑1 (MQ=255)
                                                                        gCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTATGCGGCTTCATACCGGTccc                                                                           >  2:334139/1‑149 (MQ=255)
                                                                           ccAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGAGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGc                                                                        >  2:343523/1‑149 (MQ=255)
                                                                               ccAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCTTACCGGTCCCAGCCTTc                                                                    >  2:429155/1‑149 (MQ=255)
                                                                                cAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTcc                                                                   <  1:639483/149‑1 (MQ=255)
                                                                                 aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGGGTTGTTCGGCCCTGACTTCGGCGTCGATGTCATCGTCAGGAGTCGCTTCCTGCCAGCGGCGCATGGTCAGCCTTTTTGCGGGTTCATACCCGTCCCCGCGTGCCt                                                                  >  2:733805/1‑149 (MQ=255)
                                                                                         cccGCGCGTAGGTCACGTAGCGTGTGGGCGCGATGTTTAGTTCGGCCCTGTATTCGGCGGCGTAGTCATCGTCAGGATTCGCTTCCAGCCTGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGGCTGcc                                                          <  1:128476/149‑1 (MQ=255)
                                                                                          ccGCGCCTAGGTCACGTAGCGTTTCTGCGCGATGTTTAGTTCGGCCCTGACTACGGCGGCGAAGTCATCGTCATGATTCGCTTCCAGCCAGCCGAGCATGTTCATCCATTTTGCGGATTCTTACCGGTACCAGACATACTGGACTGACa                                                         >  2:1028106/1‑149 (MQ=255)
                                                                                           cgcgcgTAGGTCACGTAGCGTATCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCTGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCGTGATGCCGGTCCCAGCCTTCCTTGTCTGGCAg                                                        >  2:180820/1‑149 (MQ=255)
                                                                                                gTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAGTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCa                                                   >  2:823101/1‑149 (MQ=255)
                                                                                                    gTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTc                                               >  2:143640/1‑149 (MQ=255)
                                                                                                    gTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTc                                               <  1:409562/149‑1 (MQ=255)
                                                                                                    gTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGACCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCCGCGCATGGTCAGCCATTTTGCGGCGTCATACCGGTCCCAGACATCCTGGTCTGCCAGTACCAGTTc                                               >  1:565341/1‑149 (MQ=255)
                                                                                                     tCACGTAGCGTTTCGGCGCGATGTTTAGCTCGGCCCTGACTTCGGCGGCGAAGTCATCTTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGGCTGCCAGTACCATTTcc                                              <  1:994706/149‑1 (MQ=255)
                                                                                                       aCGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTccac                                            <  2:385575/149‑1 (MQ=255)
                                                                                                         gTAGCGTTTCGGCGCGATGATTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCAAGCCTTCCGGGTCTGCCAGGACCATTTccccca                                          >  1:703313/1‑149 (MQ=255)
                                                                                                           aGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTATGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACg                                        <  1:334139/149‑1 (MQ=255)
                                                                                                             cGTTTCGGCGCTATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTc                                      <  1:38556/149‑1 (MQ=255)
                                                                                                             cGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGGGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTc                                      <  1:108996/149‑1 (MQ=255)
                                                                                                                ttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCGGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGGACCTTTTCCGCCACGTCaga                                   >  2:772604/1‑147 (MQ=255)
                                                                                                                 tCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAg                                  <  1:454593/149‑1 (MQ=255)
                                                                                                                    gcgcGATGTCTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGGCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGGACCATTTCCACCACGTCATAGCCg                               <  2:919833/149‑1 (MQ=255)
                                                                                                                     cgcgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCga                              <  1:10935/149‑1 (MQ=255)
                                                                                                                     cgcgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGGCTGCCAGTACCATTTCCACCACGTCATAGCCga                              <  1:1022098/149‑1 (MQ=255)
                                                                                                                       cgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCgaga                            <  2:273396/149‑1 (MQ=255)
                                                                                                                       cgATGCTTAGTTCGGCCCTGACGTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCgaga                            <  1:690119/149‑1 (MQ=255)
                                                                                                                              tAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGa                     <  1:998812/149‑1 (MQ=255)
                                                                                                                                gTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAAc                   <  1:854285/149‑1 (MQ=255)
                                                                                                                                  tCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTGCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGc                 <  2:638111/149‑1 (MQ=255)
                                                                                                                                   tgacgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTCCCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCt                >  2:106154/7‑149 (MQ=255)
                                                                                                                                      cccTGACTTCGGCGGCGAAGTCATGGTCAGGATTCTATTCCAGCCAGCGGCGCATGGTCATCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCAGGGTCGTCCAGAACCATTTCAACCACGTCATCGCCGAGATCGTCGCAAGCTcat             >  2:769130/1‑147 (MQ=255)
                                                                                                                                        cTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAGCGTCGAACGCTCCGAc           >  1:912078/1‑149 (MQ=255)
                                                                                                                                          gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGCTCGTCGAACGCTCCGACaa         >  2:518724/1‑149 (MQ=255)
                                                                                                                                          gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACaa         >  1:835507/1‑149 (MQ=255)
                                                                                                                                           aCTTCGTCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCATCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAg        <  1:179103/149‑1 (MQ=255)
                                                                                                                                           aCTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGGCAGTACCATTTCAACCCCGTCAGAGACTAGATCGTCGAACGCTCCGACAAg        >  2:68122/1‑149 (MQ=255)
                                                                                                                                                 gcggcgAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTgg  <  1:520205/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGG  >  NC_000913/4532389‑4532682

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑