Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,194,270 G→C 6.3% L458V (CTG→GTG)  PP_2799 ← class III aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,194,2700GC 6.3% 1039.4 / 40.3 311L458V (CTG→GTG) PP_2799class III aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (152/139);  new base C (10/10);  total (162/149)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

ACACTTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGA  >  NC_002947/3194121‑3194418
                                                                                                                                                     |                                                                                                                                                    
   attGGGCCATGCTGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTg                                                                                                                                                   <  2:3645509/149‑1 (MQ=255)
    ttGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGc                                                                                                                                                  >  1:6617620/1‑150 (MQ=255)
    ttGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGc                                                                                                                                                  >  2:1520894/1‑150 (MQ=255)
       ggCCATGATTGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCTTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtgt                                                                                                                                               <  1:2573893/150‑1 (MQ=255)
       ggCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTCCGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtgt                                                                                                                                               <  2:2760012/150‑1 (MQ=255)
       ggCCATGAGGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTCGCTCACAGCACCTGCtgt                                                                                                                                               <  1:2254632/150‑1 (MQ=255)
        gCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtgtg                                                                                                                                              >  1:5135466/1‑150 (MQ=255)
         ccATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGa                                                                                                                                             <  1:7950908/150‑1 (MQ=255)
         ccATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGa                                                                                                                                             <  1:3055135/150‑1 (MQ=255)
          cATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAg                                                                                                                                            >  1:800934/1‑150 (MQ=255)
           atgatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGAgc                                                                                                                                           <  2:1126065/150‑1 (MQ=255)
             gatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgc                                                                                                                                         <  2:7437644/150‑1 (MQ=255)
             gatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgc                                                                                                                                         <  1:1299507/150‑1 (MQ=255)
             gatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGAgcaa                                                                                                                                         <  1:7579351/150‑3 (MQ=255)
                ggTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGGgcaaagg                                                                                                                                      <  2:2084300/150‑6 (MQ=255)
                 gTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCAGGATTGGTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGccc                                                                                                                                     >  2:4449443/1‑150 (MQ=255)
                   tCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCaa                                                                                                                                   <  1:692754/150‑1 (MQ=255)
                      cATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATcc                                                                                                                                >  1:1671287/1‑150 (MQ=255)
                       aTGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCg                                                                                                                               >  1:8771302/1‑150 (MQ=255)
                          gCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTcg                                                                                                                            >  1:2512802/1‑150 (MQ=255)
                            aCCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTCTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTcgcg                                                                                                                          >  2:3325074/1‑150 (MQ=255)
                            aCCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTATTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTcgcg                                                                                                                          >  2:254536/1‑150 (MQ=255)
                             ccGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGt                                                                                                                         >  1:5530522/1‑150 (MQ=255)
                                     gATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagccagc                                                                                                                 <  1:3853550/150‑1 (MQ=255)
                                     gATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagccagc                                                                                                                 <  2:5535093/150‑1 (MQ=255)
                                     gATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagccagc                                                                                                                 <  2:7188487/150‑1 (MQ=255)
                                       tCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGt                                                                                                               <  2:8170204/150‑1 (MQ=255)
                                          cGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCg                                                                                                            >  2:1753819/1‑150 (MQ=255)
                                           ggAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCTTGATTGTTGTTTTACGTGAGTAGCCATCTGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGc                                                                                                           <  2:2449832/150‑1 (MQ=255)
                                           ggAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGc                                                                                                           >  2:3129846/1‑150 (MQ=255)
                                           ggAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGc                                                                                                           >  1:6471372/1‑150 (MQ=255)
                                                cGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGaa                                                                                                      <  1:1583620/150‑1 (MQ=255)
                                                 gTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAg                                                                                                     >  2:6294758/1‑150 (MQ=255)
                                                    gAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGgcg                                                                                                  <  1:930954/150‑1 (MQ=255)
                                                     aaCTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGgcgc                                                                                                 <  2:5825707/150‑1 (MQ=255)
                                                        tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGa                                                                                              <  1:4787611/150‑1 (MQ=255)
                                                         cGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGaa                                                                                             <  2:4558923/150‑1 (MQ=255)
                                                          gAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAg                                                                                            >  1:641221/1‑150 (MQ=255)
                                                          gAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAg                                                                                            >  1:7283208/1‑150 (MQ=255)
                                                           aaGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGGTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGa                                                                                           >  1:3658380/1‑150 (MQ=255)
                                                            aGCTTTCATCACTGTGTACACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGAt                                                                                          <  2:6093618/150‑1 (MQ=255)
                                                               tttCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGc                                                                                       >  1:581767/1‑150 (MQ=255)
                                                                   atcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCg                                                                                   >  1:5243593/1‑150 (MQ=255)
                                                                   atcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGAGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCg                                                                                   >  2:7936066/1‑150 (MQ=255)
                                                                    tcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGa                                                                                  >  2:5382612/1‑150 (MQ=255)
                                                                      aCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGata                                                                                >  2:3026661/1‑150 (MQ=255)
                                                                        tgtgTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATc                                                                              >  1:8130292/1‑150 (MQ=255)
                                                                         gtgtTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTATCCATCTGATGCAGCCTTTGCTCAAAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTAGATATCg                                                                             <  2:3730781/150‑1 (MQ=255)
                                                                           gtTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTc                                                                           <  2:6032539/150‑1 (MQ=255)
                                                                            ttCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTcc                                                                          >  2:8202952/1‑150 (MQ=255)
                                                                                 gcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggt                                                                     >  2:1261725/1‑150 (MQ=255)
                                                                                  cacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggtg                                                                    <  2:407447/150‑1 (MQ=255)
                                                                                   acgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggtgg                                                                   <  1:521197/150‑1 (MQ=255)
                                                                                   acgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggtgg                                                                   >  2:4707973/1‑150 (MQ=255)
                                                                                       tCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAAt                                                               >  1:7209758/1‑150 (MQ=255)
                                                                                       tCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAAt                                                               >  2:6519914/1‑150 (MQ=255)
                                                                                         aaaCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCa                                                             >  1:434020/1‑150 (MQ=255)
                                                                                           aCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAAc                                                           <  1:4191700/150‑1 (MQ=255)
                                                                                             ccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgg                                                         >  1:4095734/1‑150 (MQ=255)
                                                                                              cTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                        >  1:1107250/1‑150 (MQ=255)
                                                                                              cTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                        <  1:1728985/150‑1 (MQ=255)
                                                                                               tCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggag                                                       >  2:1565765/1‑150 (MQ=255)
                                                                                               tCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggag                                                       <  1:3262033/150‑1 (MQ=255)
                                                                                                  tGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGAc                                                    >  2:7289572/1‑150 (MQ=255)
                                                                                                    aTTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAg                                                  <  2:3992058/150‑1 (MQ=255)
                                                                                                      tgttgtTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGaa                                                <  2:1784435/150‑1 (MQ=255)
                                                                                                       gttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAg                                               <  2:6926607/150‑1 (MQ=255)
                                                                                                           ttttACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGa                                           >  2:4748509/1‑150 (MQ=255)
                                                                                                           ttttACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGa                                           >  2:2933854/1‑150 (MQ=255)
                                                                                                                    aGTAGCCATCTGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGTCAGAAGGCGATCGTGTCGt                                  <  2:1470125/150‑1 (MQ=255)
                                                                                                                     gTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTa                                 <  2:6757908/150‑1 (MQ=255)
                                                                                                                      tAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGACGATCGCCTCGATCTCGTCCTGGGTGGTAATCAACGGCGGACCGAAGGCGATCGGGTCGTAg                                >  2:6179046/1‑150 (MQ=255)
                                                                                                                        gCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGAt                              >  2:3882228/1‑150 (MQ=255)
                                                                                                                         ccATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATg                             <  2:5907152/150‑1 (MQ=255)
                                                                                                                          cATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATgg                            >  1:5392708/1‑150 (MQ=255)
                                                                                                                          cATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATgg                            >  1:4487719/1‑150 (MQ=255)
                                                                                                                           aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGgc                           >  1:7393838/1‑150 (MQ=255)
                                                                                                                                aTGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGaa                      <  2:4844528/150‑1 (MQ=255)
                                                                                                                                 tGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAAt                     <  1:6177921/150‑1 (MQ=255)
                                                                                                                                 tGCAGCCTTTGCTCACAGCAGGTGCGGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCACGATATCGTCCTGGGTGGTAATCAAAGCAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAAt                     >  1:3145068/1‑150 (MQ=255)
                                                                                                                                   cAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGa                   >  1:4942532/1‑150 (MQ=255)
                                                                                                                                   cAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGa                   >  1:7110959/1‑150 (MQ=255)
                                                                                                                                      ccTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCa                <  2:5201297/150‑1 (MQ=255)
                                                                                                                                       gcgcgCTCACAGCACCTGCTGTGAGCACGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                        <  1:3514882/105‑1 (MQ=39)
                                                                                                                                       gcgcgCTCACAGCACCTGCTGTGAGCACGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                        >  2:3514882/5‑109 (MQ=39)
                                                                                                                                        tttGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGc              <  2:6617620/150‑1 (MQ=255)
                                                                                                                                           gCTCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCa           >  1:6961106/1‑150 (MQ=255)
                                                                                                                                             tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATg         >  1:8226197/1‑150 (MQ=255)
                                                                                                                                             tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATg         >  1:7788893/1‑150 (MQ=255)
                                                                                                                                             tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATg         >  1:3601271/1‑150 (MQ=255)
                                                                                                                                              cacaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGc        <  1:1965816/150‑1 (MQ=255)
                                                                                                                                              cacaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGc        <  2:7765546/150‑1 (MQ=255)
                                                                                                                                               acaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCt       >  2:8140448/1‑150 (MQ=255)
                                                                                                                                                ccgcagGGGCTGTGAGCGCGCCCACGCCGCCGCGTCAGCCCGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGCTATCGTCCTGGGTGGTAATCACCGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGATTGATCAGCCCATGCtt      <  2:8649027/148‑1 (MQ=255)
                                                                                                                                                cagcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCtt      <  1:4449443/150‑1 (MQ=255)
                                                                                                                                                  gcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgg                                                         <  2:8514592/97‑1 (MQ=255)
                                                                                                                                                  gcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTtgt    <  1:2642040/150‑1 (MQ=255)
                                                                                                                                                    agGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGa  <  1:8675402/150‑1 (MQ=255)
                                                                                                                                                    agGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGa  <  1:7670424/150‑1 (MQ=255)
                                                                                                                                                    agGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGa  <  1:1172560/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACACTTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGA  >  NC_002947/3194121‑3194418

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