Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,713,585 T→C 6.4% T1162A (ACC→GCC) ‡ metH ← cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,713,5850TC 6.4% 630.9 / 20.9 201T1162A (ACC→GCC) ‡metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (88/100);  new base C (10/3);  total (98/103)
Fisher's exact test for biased strand distribution p-value = 4.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 203 total aligned reads displayed.

ACGCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCAGTGCTCTTTA  >  NC_002947/2713436‑2713733
                                                                                                                                                     |                                                                                                                                                    
aCGCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGt                                                                                                                                                      <  1:2661592/150‑1 (MQ=255)
aCGCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGt                                                                                                                                                      <  2:3270824/150‑1 (MQ=255)
 cGCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTg                                                                                                                                                     <  1:4874230/150‑1 (MQ=255)
      gATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATc                                                                                                                                                >  2:740273/1‑150 (MQ=255)
        tCTGGTCCTTGTCGCCCTTGCCTACGGCAAAATTCTTCGCCTGCGGGGGGGCGTAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCACCGCCATCgg                                                                                                                                              <  1:609828/150‑2 (MQ=255)
          tGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGc                                                                                                                                            >  2:1554617/1‑150 (MQ=255)
             tCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAgg                                                                                                                                         >  1:2062241/1‑150 (MQ=255)
               cTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCg                                                                                                                                       >  1:1154148/1‑150 (MQ=255)
               cTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCg                                                                                                                                       <  1:555928/150‑1 (MQ=255)
                ttGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCgg                                                                                                                                      >  1:3577671/1‑150 (MQ=255)
                    cGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACa                                                                                                                                  <  1:4008794/150‑1 (MQ=255)
                     gACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAAcgg                                                                                                                                 >  2:1374422/1‑148 (MQ=255)
                     gACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACag                                                                                                                                 <  1:4810929/150‑1 (MQ=255)
                      aCCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACaga                                                                                                                                <  1:1900249/150‑1 (MQ=255)
                      aCCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACaga                                                                                                                                <  1:4036301/150‑1 (MQ=255)
                       ccTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACagag                                                                                                                               >  1:88249/1‑150 (MQ=255)
                         ttGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGtt                                                                                                                             >  2:3148397/1‑150 (MQ=255)
                           gCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTc                                                                                                                           <  2:1690991/150‑1 (MQ=255)
                             cTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCt                                                                                                                         <  1:4072402/150‑1 (MQ=255)
                              tACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCtt                                                                                                                        <  2:3005265/150‑1 (MQ=255)
                                 ggCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTc                                                                                                                     <  1:2946255/150‑1 (MQ=255)
                                    aaaaTACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAgt                                                                                                                  >  2:447075/1‑150 (MQ=255)
                                    aaaaTACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAgt                                                                                                                  >  2:4766829/1‑150 (MQ=255)
                                    aaaaTACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAgt                                                                                                                  >  2:5048174/1‑150 (MQ=255)
                                    aaaaTACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAgt                                                                                                                  >  1:4200426/1‑150 (MQ=255)
                                           ttCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCg                                                                                                           >  2:2517831/1‑150 (MQ=255)
                                           ttCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGATTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCg                                                                                                           >  2:1572934/1‑150 (MQ=255)
                                            tCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCgg                                                                                                          <  2:850997/150‑1 (MQ=255)
                                                 tGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCa                                                                                                                    >  1:917772/1‑135 (MQ=255)
                                                 tGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCa                                                                                                                    <  2:917772/135‑1 (MQ=255)
                                                  gCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAgg                                                                                                    >  1:1411365/1‑150 (MQ=255)
                                                      gTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCgg                                                                                                <  1:319098/150‑1 (MQ=255)
                                                         ggCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATa                                                                                             >  2:4071070/1‑150 (MQ=255)
                                                          gCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATaa                                                                                            <  1:2685863/150‑1 (MQ=255)
                                                          gCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATaa                                                                                            <  1:4539139/150‑1 (MQ=255)
                                                          gCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATaa                                                                                            >  1:2253585/1‑150 (MQ=255)
                                                          gCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGGGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTGTCAGTGTGGTCGGGGCAGGCCGGATaa                                                                                            >  2:3214665/1‑150 (MQ=255)
                                                            gAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAAcc                                                                                          >  1:3789934/1‑150 (MQ=255)
                                                                  ccAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGgcgc                                                                                    >  2:577698/1‑150 (MQ=255)
                                                                  ccAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGgcgc                                                                                    >  2:1596676/1‑150 (MQ=255)
                                                                     gCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAgg                                                                                 <  2:1423394/150‑1 (MQ=255)
                                                                     gCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAgg                                                                                 <  2:2650674/150‑1 (MQ=255)
                                                                       cGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGc                                                                               >  2:239882/1‑150 (MQ=255)
                                                                       cGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGc                                                                               <  2:203540/150‑1 (MQ=255)
                                                                       cGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGc                                                                               <  1:3306067/150‑1 (MQ=255)
                                                                                   ccgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGt                                                                   <  2:2100182/150‑1 (MQ=255)
                                                                                    cgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTa                                                                  <  2:2721465/150‑1 (MQ=255)
                                                                                    cgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTa                                                                  <  2:2451540/150‑1 (MQ=255)
                                                                                     gccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTAc                                                                 <  1:244488/150‑1 (MQ=255)
                                                                                        gggAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGt                                                              <  1:255051/150‑1 (MQ=255)
                                                                                          gAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCt                                                                     <  1:2530299/141‑1 (MQ=255)
                                                                                           aaCATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTcc                                                           <  1:2660027/150‑1 (MQ=255)
                                                                                           aaCATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTcc                                                           <  2:427174/150‑1 (MQ=255)
                                                                                            aCATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCt                                                          >  1:1433684/1‑150 (MQ=255)
                                                                                             cATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCTGGGCGGATGCCGCTGTACTGTTCCtt                                                         <  1:3487009/150‑1 (MQ=255)
                                                                                             cATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCtt                                                         <  1:2000229/150‑1 (MQ=255)
                                                                                              aTCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTg                                                        <  1:4975142/150‑1 (MQ=255)
                                                                                               tCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGcc                                                                                                  <  1:478574/107‑1 (MQ=255)
                                                                                               tCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGa                                                       >  2:2120007/1‑150 (MQ=255)
                                                                                               tCGCGAAGTGCTCGCTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGcc                                                                                                  >  2:478574/1‑107 (MQ=255)
                                                                                                 gcgAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATc                                                     >  2:1770363/1‑150 (MQ=255)
                                                                                                         cTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTc                                             <  1:3939159/150‑1 (MQ=255)
                                                                                                          tCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCg                                            >  1:1461226/1‑150 (MQ=255)
                                                                                                          tCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCg                                            >  1:4511071/1‑150 (MQ=255)
                                                                                                          tCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGAGGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGTTCGGGGCAGGCCGGATAACCCGGCGCAGGTCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCg                                            <  2:1328726/150‑1 (MQ=255)
                                                                                                            ggTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGtt                                          <  2:3942497/150‑1 (MQ=255)
                                                                                                               gAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTc                                       <  1:1054106/150‑1 (MQ=255)
                                                                                                                aGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTcc                                      <  1:4284713/150‑1 (MQ=255)
                                                                                                                 ggTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCa                                     <  2:3582058/150‑1 (MQ=255)
                                                                                                                 ggTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCa                                     <  1:810915/150‑1 (MQ=255)
                                                                                                                  gTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAg                                    <  2:2580934/150‑1 (MQ=255)
                                                                                                                  gTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAg                                    >  2:121926/1‑150 (MQ=255)
                                                                                                                  gTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAg                                    >  1:3703145/1‑150 (MQ=255)
                                                                                                                   tACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATACCCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAgg                                   <  2:3734481/150‑1 (MQ=255)
                                                                                                                    cccccccGCTCGGCCCGGTTTCACCGCTGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTGTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGt                                  <  2:4586961/147‑1 (MQ=255)
                                                                                                                    acacCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGt                                  <  2:817119/150‑1 (MQ=255)
                                                                                                                        cccGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCtc                              >  1:933642/1‑150 (MQ=255)
                                                                                                                             tCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTcgc                         <  1:3148397/150‑1 (MQ=255)
                                                                                                                             tCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTcgc                         <  1:3058803/150‑1 (MQ=255)
                                                                                                                               ggCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCgg                       <  2:1866570/150‑1 (MQ=255)
                                                                                                                               ggCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCgg                       <  1:849577/150‑1 (MQ=255)
                                                                                                                                gCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCggg                      <  1:3731843/150‑1 (MQ=255)
                                                                                                                                  ccGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCa                    >  1:2325817/1‑150 (MQ=255)
                                                                                                                                    ggTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATa                  <  2:4952836/150‑1 (MQ=255)
                                                                                                                                     gTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAg                 >  1:4733476/1‑150 (MQ=255)
                                                                                                                                        tgctcGATGGCACCGCCATCGGGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGccc              >  1:5047193/5‑150 (MQ=255)
                                                                                                                                        tgctcGATGGCACCGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTCCTGTTCCTTGATCAGCGCCGCGTTGTCCAGGTGCTCGTCGCGGGCATAGccc              >  2:36635/5‑150 (MQ=255)
                                                                                                                                        tCACCGATGGCGGTTCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGccc              <  2:2193608/150‑1 (MQ=255)
                                                                                                                                         gctcGATGGCACCGCCATCGGGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGcccc             >  1:1221728/4‑150 (MQ=255)
                                                                                                                                         cACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGcccc             <  1:1554617/150‑1 (MQ=255)
                                                                                                                                          aCCGATGGCGGTGCCATCGAGCAGGCGGAACAGTGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCa            >  2:4228654/1‑150 (MQ=255)
                                                                                                                                          aCCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCa            >  2:3889828/1‑150 (MQ=255)
                                                                                                                                          aCCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCa            >  1:968292/1‑150 (MQ=255)
                                                                                                                                          aCCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCa            >  2:3890082/1‑150 (MQ=255)
                                                                                                                                          aCCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCa            <  2:2322980/150‑1 (MQ=255)
                                                                                                                                          aCCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCa            >  1:3013771/1‑150 (MQ=255)
                                                                                                                                           tcGATGGCACCGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCAg           >  2:2563371/2‑150 (MQ=255)
                                                                                                                                            cGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCAGt          <  1:1318986/150‑1 (MQ=255)
                                                                                                                                             gATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCAGTg         >  2:1544497/1‑150 (MQ=255)
                                                                                                                                                   ggTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCAGTGCTCttt   <  2:3703145/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCAGTGCTCTTTA  >  NC_002947/2713436‑2713733

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