Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,711,618 G→A 11.1% T1162T (ACC→ACT) ‡ metH ← cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,711,6180GA11.1% 241.4 / 15.0 91T1162T (ACC→ACT) ‡metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base G (40/40);  new base A (4/6);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

ACGCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCAG  >  NC_002947/2711471‑2711759
                                                                                                                                                    |                                                                                                                                              
aCGCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCg                                                                                                                                                <  7:170516/148‑1 (MQ=255)
 cGCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGt                                                                                                                                              <  7:184877/149‑1 (MQ=255)
  gCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTg                                                                                                                                             >  3:163830/1‑149 (MQ=255)
   cTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCGGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGc                                                                                                                                            >  7:262026/1‑149 (MQ=255)
    tCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACACCGCTCGGCCCGGTTTCACCGATGGCGGTGc                                                                                                                                            >  3:249186/1‑148 (MQ=255)
     cGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCa                                                                                                                                          >  1:105534/1‑149 (MQ=255)
       aTCTGGTCCTTGTAGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATc                                                                                                                                        >  6:335586/1‑149 (MQ=255)
          tGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAg                                                                                                                                     >  6:308238/1‑149 (MQ=255)
            ggccTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGGGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCACCGCCATCgg                                                                                                                                      <  8:355115/144‑2 (MQ=255)
            gTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCACCGCCATCggt                                                                                                                                     <  3:187219/147‑3 (MQ=255)
            gTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCACCGCCATCgg                                                                                                                                      <  2:35614/146‑2 (MQ=255)
            gTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCACCGCCATCgg                                                                                                                                      <  2:29592/146‑2 (MQ=255)
            gTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCACCGCCATCgg                                                                                                                                      <  3:237526/146‑2 (MQ=255)
              ccTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATc                                                                                                                                        <  7:156692/142‑1 (MQ=255)
              ccTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATc                                                                                                                                        >  8:156692/1‑142 (MQ=255)
              ccTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAgg                                                                                                                                 >  3:275543/1‑149 (MQ=255)
               cTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGc                                                                                                                                <  3:260349/149‑1 (MQ=255)
                      aCCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACag                                                                                                                        <  2:379930/149‑1 (MQ=255)
                            ccTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGttt                                                                                                                   <  8:293102/148‑1 (MQ=255)
                            ccTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTc                                                                                                                  <  7:15587/149‑1 (MQ=255)
                             cTACGGCAAAATACTTCGCCTGCGGGTGGGCGGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTc                                                                                                                  >  1:41283/1‑149 (MQ=255)
                             cTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTcc                                                                                                                 >  4:80412/1‑149 (MQ=255)
                              tACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCACCCCCTGACCCCCGCCCCCGGCAACATCCCGAAGGGCACGGTGAGGTCCACCCCGCTCGGCCCGGTTAAACCGATGGCGGTGCCAGCGAGCCGGCGGAACAGAGTTTCCt                                                                                                                >  8:159911/1‑149 (MQ=255)
                                cGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCTGTTTCACCGATGGCGGTGCCATCGAGCAGGCGTAACAGAGTTTCCt                                                                                                                <  7:58960/147‑1 (MQ=255)
                                 ggCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCtttt                                                                                                             <  5:335586/149‑1 (MQ=255)
                                 ggCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCtttt                                                                                                             >  1:212624/1‑149 (MQ=255)
                                   cAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCtttt                                                                                                             >  2:303505/1‑147 (MQ=255)
                                   cAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCa                                                                                                           <  5:301145/149‑1 (MQ=255)
                                     aaaTACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAgt                                                                                                         <  6:400071/149‑1 (MQ=255)
                                       aTACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAgtgt                                                                                                       >  7:263895/1‑149 (MQ=255)
                                        tACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGtgtg                                                                                                      <  4:275543/149‑1 (MQ=255)
                                             cGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTc                                                                                                   <  5:193847/147‑1 (MQ=255)
                                                     ggTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGcc                                                                                         <  7:278917/149‑1 (MQ=255)
                                                      gTGGGCGAAGTACCAGACGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCg                                                                                        >  8:228159/1‑149 (MQ=255)
                                                         ggCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGa                                                                                      >  8:199116/1‑148 (MQ=255)
                                                                   ccAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGgcg                                                                            >  2:355436/1‑149 (MQ=255)
                                                                   ccAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGgcg                                                                            >  7:370124/1‑149 (MQ=255)
                                                                     aGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCa                                                                          <  4:163830/149‑1 (MQ=255)
                                                                        cGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAggg                                                                       >  3:206284/1‑149 (MQ=255)
                                                                              ccgccgccgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATg                                                                 <  7:228159/149‑1 (MQ=255)
                                                                               cgccgccgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTc                                                                                                                  <  1:52301/99‑1 (MQ=255)
                                                                               cgccgccgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTc                                                                                                                  >  2:52301/1‑99 (MQ=255)
                                                                               cgccgccgccgGGAACATCGCGAAGTGATCGGTGAGGTACACCCCGATCGGCACGGTTTCACCGATGGCCGTGCCATCGAGCAGGCGGAACAGAGTTACCTTTTCAGAGTGGTCGGGGAAGGCCGGATAACCCGGCGCAGGGCGGAAGc                                                                >  2:64027/1‑149 (MQ=255)
                                                                                 ccgccgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCg                                                              <  2:41283/149‑1 (MQ=255)
                                                                                   gccgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCt                                                            <  3:80412/149‑1 (MQ=255)
                                                                                    ccgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTg                                                           <  4:206284/149‑1 (MQ=255)
                                                                                     cgccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTg                                                           >  7:74202/1‑148 (MQ=255)
                                                                                     cgccgGGAACATCGCGAAGTGCTCGGGGAGGTACACCCCGCTCGGCCCGGTCTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGTCCGGATAACCCGGCGCAGGGCGGATGCCGCt                                                            <  1:64027/147‑1 (MQ=255)
                                                                                       ccgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTAc                                                        >  8:241396/1‑149 (MQ=255)
                                                                                        cgGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACt                                                       <  1:355436/149‑1 (MQ=255)
                                                                                         gggAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTg                                                      >  7:207718/1‑149 (MQ=255)
                                                                                             aCATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTcc                                                  <  7:199116/149‑1 (MQ=255)
                                                                                               aTCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCAGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTACTTTTCAGTGTGGTCGGGGCAGGCCGGAAAACCCGCCGCAGCGCGGATGCCGCTCTACTGTTCCtt                                                >  8:23169/1‑149 (MQ=255)
                                                                                                tCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATg                                                                 >  1:185206/1‑131 (MQ=255)
                                                                                                tCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATg                                                                 <  2:185206/131‑1 (MQ=255)
                                                                                                  gcgAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGAt                                             <  2:105534/149‑1 (MQ=255)
                                                                                                  gcgAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGAt                                             >  3:282130/1‑149 (MQ=255)
                                                                                                      aGTGCTCGGTGAGGTACACCCCGCGCGGCCCGGTTTCACCGATTGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGTGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAgc                                         <  4:249186/149‑1 (MQ=255)
                                                                                                          cTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTAc                                                        <  1:511/130‑1 (MQ=255)
                                                                                                          cTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTAc                                                        >  2:511/1‑130 (MQ=255)
                                                                                                          cTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCt                                     >  2:35830/1‑149 (MQ=255)
                                                                                                              gTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGtt                                 <  5:308238/149‑1 (MQ=255)
                                                                                                                  ggTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTcc                             <  1:303505/149‑1 (MQ=255)
                                                                                                                   gTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTc                              >  1:298542/1‑147 (MQ=255)
                                                                                                                     acacCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCGTCGTTGTCCAgg                          >  6:74523/1‑149 (MQ=255)
                                                                                                                      cacCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGt                         <  4:2581/149‑1 (MQ=255)
                                                                                                                       acCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTc                                    >  6:109222/1‑137 (MQ=255)
                                                                                                                       acCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTc                                    <  5:109222/137‑1 (MQ=255)
                                                                                                                         cccGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCt                      >  2:230037/1‑149 (MQ=255)
                                                                                                                           cGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTc                              <  5:286723/139‑1 (MQ=255)
                                                                                                                           cGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTc                              >  6:286723/1‑139 (MQ=255)
                                                                                                                                 gCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCgg              <  7:241396/149‑1 (MQ=255)
                                                                                                                                  cccGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTg                                               >  2:21370/1‑115 (MQ=255)
                                                                                                                                  cccGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTg                                               <  1:21370/115‑1 (MQ=255)
                                                                                                                                      gTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATa         <  8:263895/149‑1 (MQ=255)
                                                                                                                                      gTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATa         <  8:207718/149‑1 (MQ=255)
                                                                                                                                       tttCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGGCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAg        <  8:370124/149‑1 (MQ=255)
                                                                                                                                       cctgctcGATGGCACCGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAg        >  8:37723/7‑149 (MQ=255)
                                                                                                                                         tgctcGATGGCACCGCCATCGAGCAGGCGGAACGAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCt                      >  1:252193/5‑134 (MQ=255)
                                                                                                                                         tgctcGATGGCACCGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCt                      <  2:252193/129‑1 (MQ=255)
                                                                                                                                         tgctcGATGGCACCGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGcc      >  6:19902/5‑149 (MQ=255)
                                                                                                                                         tgcccGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCgg              >  6:287975/4‑141 (MQ=255)
                                                                                                                                         tgcccGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCgg              <  5:287975/138‑1 (MQ=255)
                                                                                                                                          cACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGGCGGATAACCCGGCGCAGGGCGGATGCCGCTGTAGTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGcc      <  7:23169/148‑1 (MQ=255)
                                                                                                                                           ctcGATGGCACCGCCATCGGGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGccc     >  1:75744/3‑148 (MQ=255)
                                                                                                                                           aCCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGccc     >  5:283664/1‑148 (MQ=255)
                                                                                                                                            ccGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGtt                                 >  8:331392/1‑119 (MQ=255)
                                                                                                                                            ccGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGtt                                 <  7:331392/119‑1 (MQ=255)
                                                                                                                                            ccGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGccc     >  1:92160/1‑147 (MQ=255)
                                                                                                                                             cGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCt                                                 >  1:388834/1‑102 (MQ=255)
                                                                                                                                             cGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCt                                                 <  2:388834/102‑1 (MQ=255)
                                                                                                                                             cGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCAg  >  1:2332/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ACGCTCGATCTGGTCCTTGTCGACCTTGCCTACGGCAAAATACTTCGCCTGCGGGTGGGCGAAGTACCAGCCGCTGACCGCCGCCGCCGGGAACATCGCGAAGTGCTCGGTGAGGTACACCCCGCTCGGCCCGGTTTCACCGATGGCGGTGCCATCGAGCAGGCGGAACAGAGTTTCCTTTTCAGTGTGGTCGGGGCAGGCCGGATAACCCGGCGCAGGGCGGATGCCGCTGTACTGTTCCTTGATCAGCGCCTCGTTGTCCAGGTGCTCGTCGCGGGCATAGCCCCAG  >  NC_002947/2711471‑2711759

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