Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,194,271 G→C 5.9% H457Q (CAC→CAG PP_2799 ← class III aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,194,2710GC 5.9% 1127.7 / 39.5 335H457Q (CAC→CAGPP_2799class III aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (156/159);  new base C (10/10);  total (166/169)
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 339 total aligned reads displayed.

ACTTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGAGC  >  NC_002947/3194123‑3194420
                                                                                                                                                    |                                                                                                                                                     
acTTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGt                                                                                                                                                      >  2:7451525/1‑150 (MQ=255)
  ttGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGc                                                                                                                                                    >  1:8785121/1‑150 (MQ=255)
    gggCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtg                                                                                                                                                  <  2:2157177/150‑1 (MQ=255)
     ggCCATGATTGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCTTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtgt                                                                                                                                                 <  1:2573893/150‑1 (MQ=255)
     ggCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtgt                                                                                                                                                 <  2:5987374/150‑1 (MQ=255)
     ggCCATGAGGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTCGCTCACAGCACCTGCtgt                                                                                                                                                 <  1:2254632/150‑1 (MQ=255)
      gCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtgtg                                                                                                                                                >  1:5135466/1‑150 (MQ=255)
       ccATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGa                                                                                                                                               <  1:7951021/150‑1 (MQ=255)
       ccATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGa                                                                                                                                               <  1:5227031/150‑1 (MQ=255)
       ccATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGa                                                                                                                                               <  1:3055135/150‑1 (MQ=255)
        cATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAg                                                                                                                                              >  1:800934/1‑150 (MQ=255)
           gatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgc                                                                                                                                           <  1:1299507/150‑1 (MQ=255)
           gatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGAgcaa                                                                                                                                           <  1:7579351/150‑3 (MQ=255)
            atgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTCGCTCACAGCACCTGCTGTGAgcaaa                                                                                                                                          <  2:3805341/150‑4 (MQ=255)
               gTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCAGGATTGGTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGccc                                                                                                                                       >  2:4449443/1‑150 (MQ=255)
                 tCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCaa                                                                                                                                     <  1:692754/150‑1 (MQ=255)
                     aTGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCg                                                                                                                                 >  2:7049616/1‑150 (MQ=255)
                        gCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTcg                                                                                                                              >  1:2512802/1‑150 (MQ=255)
                           ccGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGt                                                                                                                           >  1:5530522/1‑150 (MQ=255)
                                 ccGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagcca                                                                                                                     >  1:6636379/1‑150 (MQ=255)
                                   gATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagccagc                                                                                                                   <  2:3087868/150‑1 (MQ=255)
                                   gATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagccagc                                                                                                                   <  1:3853550/150‑1 (MQ=255)
                                   gATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagccagc                                                                                                                   >  2:1172560/1‑150 (MQ=255)
                                        cGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCg                                                                                                              >  1:972000/1‑150 (MQ=255)
                                         ggAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGc                                                                                                             <  2:167459/150‑1 (MQ=255)
                                         ggAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGc                                                                                                             >  1:6471372/1‑150 (MQ=255)
                                           aaaCGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTc                                                                                                           <  2:7901074/150‑1 (MQ=255)
                                               gTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAg                                                                                                       >  2:6294758/1‑150 (MQ=255)
                                                  gAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGgcg                                                                                                    <  2:1294887/150‑1 (MQ=255)
                                                   aaCTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGgcgc                                                                                                   <  2:5825707/150‑1 (MQ=255)
                                                    aCTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCt                                                                                                  <  2:7337056/150‑1 (MQ=255)
                                                      tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGa                                                                                                <  1:4787611/150‑1 (MQ=255)
                                                      tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCAGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGa                                                                                                >  2:2624059/1‑150 (MQ=255)
                                                       cGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGaa                                                                                               <  2:4558923/150‑1 (MQ=255)
                                                        gAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAg                                                                                              >  1:641221/1‑150 (MQ=255)
                                                        gAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCACTCGAAGGCGCTGAAg                                                                                              >  1:6963899/1‑150 (MQ=255)
                                                         aaGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGgcgatcttca                                                                                             <  2:3285046/150‑8 (MQ=255)
                                                         aaGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGGTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGa                                                                                             >  1:3658380/1‑150 (MQ=255)
                                                          aGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGAt                                                                                            <  2:2512802/150‑1 (MQ=255)
                                                          aGCTTTCATCACTGTGTACACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGAt                                                                                            <  2:6093618/150‑1 (MQ=255)
                                                             tttCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGc                                                                                         >  1:581767/1‑150 (MQ=255)
                                                                catcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagcca                                                                                                                     <  2:7197426/119‑1 (MQ=255)
                                                                 atcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCg                                                                                     >  1:5243593/1‑150 (MQ=255)
                                                                  tcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGa                                                                                    >  2:5382612/1‑150 (MQ=255)
                                                                     cTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGatat                                                                                 >  1:6298373/1‑150 (MQ=255)
                                                                      tgtgTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATc                                                                                >  1:8130292/1‑150 (MQ=255)
                                                                          ttCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTcc                                                                            >  2:1782333/1‑150 (MQ=255)
                                                                               gcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggt                                                                       >  2:2960262/1‑150 (MQ=255)
                                                                                 acgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggtgg                                                                     <  1:4452202/150‑1 (MQ=255)
                                                                                     tCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAAt                                                                 <  2:4636675/150‑1 (MQ=255)
                                                                                     tCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAAt                                                                 >  1:7209758/1‑150 (MQ=255)
                                                                                       aaaCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCa                                                               >  1:434020/1‑150 (MQ=255)
                                                                                       aaaCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCa                                                               >  2:2011640/1‑150 (MQ=255)
                                                                                         aCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAAc                                                             <  1:4191700/150‑1 (MQ=255)
                                                                                            cTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                          >  1:1107250/1‑150 (MQ=255)
                                                                                            cTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                          <  1:3465785/150‑1 (MQ=255)
                                                                                            cTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                          <  2:152808/150‑1 (MQ=255)
                                                                                                tGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGAc                                                      >  2:4139256/1‑150 (MQ=255)
                                                                                                 gATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACa                                                     <  2:6547509/150‑1 (MQ=255)
                                                                                                         ttttACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGa                                             >  2:4748509/1‑150 (MQ=255)
                                                                                                           ttACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATc                                           <  1:1611408/150‑1 (MQ=255)
                                                                                                           ttACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATc                                           <  2:3710100/150‑1 (MQ=255)
                                                                                                             aCGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgt                                         <  2:6767186/150‑1 (MQ=255)
                                                                                                                  aGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCACTCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGt                                    <  2:1470180/150‑1 (MQ=255)
                                                                                                                   gTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTa                                   <  1:7324581/150‑1 (MQ=255)
                                                                                                                   gTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTa                                   <  1:4107774/150‑1 (MQ=255)
                                                                                                                    tAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGACGATCGCCTCGATCTCGTCCTGGGTGGTAATCAACGGCGGACCGAAGGCGATCGGGTCGTAg                                  >  2:6179046/1‑150 (MQ=255)
                                                                                                                      gCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGAt                                >  2:3882228/1‑150 (MQ=255)
                                                                                                                       ccATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATg                               >  2:7992169/1‑150 (MQ=255)
                                                                                                                       ccATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATg                               <  2:4019189/150‑1 (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:1753819/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)
                                                                                                                                   gCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCCGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGAGCGTGTCGTAGATGGCGCGACTGATc                   >  2:3559054/1‑150 (MQ=255)
                                                                                                                                    ccTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCa                  <  2:5201297/150‑1 (MQ=255)
                                                                                                                                     gcgcgCTCACAGCACCTGCTGTGAGCACGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                          <  1:3514882/105‑1 (MQ=39)
                                                                                                                                      tttGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGc                <  2:5530522/150‑1 (MQ=255)
                                                                                                                                      tttGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGc                <  2:6617620/150‑1 (MQ=255)
                                                                                                                                           tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATg           >  1:7788893/1‑150 (MQ=255)
                                                                                                                                           tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATg           >  2:8183492/1‑150 (MQ=255)
                                                                                                                                           tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATg           >  1:3601271/1‑150 (MQ=255)
                                                                                                                                           tCACAGCACCTGCTGTGAGCACGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATg           >  1:1702178/1‑150 (MQ=255)
                                                                                                                                            cacaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGc          <  1:1965959/150‑1 (MQ=255)
                                                                                                                                             acaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCt         >  2:2816406/1‑150 (MQ=255)
                                                                                                                                             acaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCt         >  2:8140519/1‑150 (MQ=255)
                                                                                                                                             acaGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCt         >  2:3407473/1‑150 (MQ=255)
                                                                                                                                              ccgcagGGGCTGTGAGCGCGCCCACGCCGCCGCGTCAGCCCGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGCTATCGTCCTGGGTGGTAATCACCGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGATTGATCAGCCCATGCtt        <  2:8649027/148‑1 (MQ=255)
                                                                                                                                              cagcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCtt        <  1:7753037/150‑1 (MQ=255)
                                                                                                                                                gcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgg                                                           >  1:8514592/1‑97 (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:1172560/150‑1 (MQ=255)
                                                                                                                                                    gTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGAGc  >  2:8265623/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACTTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGAGC  >  NC_002947/3194123‑3194420

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: