Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,355,903 G→T 5.1% H62N (CAC→AAC)  pnp ← polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,355,9030GT 5.1% 589.3 / 14.8 175H62N (CAC→AAC) pnppolyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (97/69);  new base T (7/2);  total (104/71)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 176 total aligned reads displayed.

CTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACGGTCGATT  >  NC_002947/5355756‑5356052
                                                                                                                                                   |                                                                                                                                                     
 ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                                                                                                    >  1:2754549/1‑150 (MQ=255)
 ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                                                                                                    >  1:4942241/1‑150 (MQ=255)
      tCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACa                                                                                                                                               >  1:1865252/1‑150 (MQ=255)
        aTGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGc                                                                                                                                             >  2:541292/1‑150 (MQ=255)
          gAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCgg                                                                                                                                           >  2:3214770/1‑150 (MQ=255)
          gAAGCCTTCAGGTAACAGCGGACGGATCTGACTGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCTGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCgg                                                                                                                                           <  2:3143473/150‑1 (MQ=255)
            aGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGga                                                                                                                                         <  1:3348745/150‑1 (MQ=255)
             gCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaa                                                                                                                                        <  2:3693730/150‑1 (MQ=255)
                ttCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaagaa                                                                                                                                     >  1:408538/1‑150 (MQ=255)
                ttCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaagaa                                                                                                                                     >  1:1121114/1‑150 (MQ=255)
                ttCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaagaa                                                                                                                                     >  2:1425216/1‑150 (MQ=255)
                atcaacgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGgaagaa                                                                                                                                     >  1:2691825/7‑150 (MQ=255)
                  caacgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGaaat                                                                                                                                   <  2:2932017/146‑2 (MQ=255)
                   cGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc                                                                                                                                  <  2:3343106/150‑1 (MQ=255)
                   cGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc                                                                                                                                  <  2:4019775/150‑1 (MQ=255)
                   cGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc                                                                                                                                  <  2:1865252/150‑1 (MQ=255)
                      gAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCtt                                                                                                                               <  1:446266/150‑1 (MQ=255)
                        aaaGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAATCCtgac                                                                                                                             <  1:413123/148‑4 (MQ=255)
                        aCAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGc                                                                                                                             <  1:3742750/150‑1 (MQ=255)
                             ggacggaTCGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGAt                                                                                                                       >  1:966140/1‑150 (MQ=255)
                                 ggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAg                                                                                                                    >  2:4388422/1‑150 (MQ=255)
                                 ggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAg                                                                                                                    >  2:3914586/1‑150 (MQ=255)
                                 ggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAg                                                                                                                    >  2:4388522/1‑150 (MQ=255)
                                  gaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGc                                                                                                                   <  1:2885130/150‑1 (MQ=255)
                                       gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGtt                                                                                                              >  1:3719564/1‑150 (MQ=255)
                                       gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGtt                                                                                                              <  1:216626/150‑1 (MQ=255)
                                           gTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGgc                                                                                                          <  1:3600246/150‑1 (MQ=255)
                                            tCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGgcg                                                                                                         >  1:806879/1‑150 (MQ=255)
                                              gATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGcc                                                                                                       >  1:2613649/1‑150 (MQ=255)
                                              gATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGcc                                                                                                       >  1:170731/1‑150 (MQ=255)
                                                  aGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTAcc                                                                                                   <  1:4235511/150‑1 (MQ=255)
                                                   ggCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACca                                                                                                  <  2:1136853/150‑1 (MQ=255)
                                                     cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                                                                                                <  1:630909/150‑1 (MQ=255)
                                                     cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                                                                                                <  1:2316541/150‑1 (MQ=255)
                                                     cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                                                                                                <  1:3831709/150‑1 (MQ=255)
                                                        gAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTc                                                                                             >  1:1024814/1‑150 (MQ=255)
                                                         aGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCa                                                                                            >  2:2668600/1‑150 (MQ=255)
                                                            tCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCa                                                                                         >  2:4481867/1‑150 (MQ=255)
                                                            tCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCa                                                                                         >  2:1904591/1‑150 (MQ=255)
                                                              aGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGc                                                                                       >  1:1985135/1‑150 (MQ=255)
                                                                 gTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACg                                                                                    <  1:3936762/150‑1 (MQ=255)
                                                                  tctcTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACgg                                                                                   <  2:2349912/150‑1 (MQ=255)
                                                                   ctctTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGt                                                                                  <  1:79626/150‑1 (MQ=255)
                                                                   ctctTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGt                                                                                  <  1:2584887/150‑1 (MQ=255)
                                                                        tctcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAcgt                                                                             >  1:1239912/1‑150 (MQ=255)
                                                                         ctcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAcgtc                                                                            <  2:4934414/150‑1 (MQ=255)
                                                                          tcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAcgtcg                                                                           <  2:682768/150‑1 (MQ=255)
                                                                            ggAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGtt                                                                         <  2:2797584/150‑1 (MQ=255)
                                                                            ggAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGtt                                                                         >  2:1369930/1‑150 (MQ=255)
                                                                              aaGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGt                                                                       <  2:1250222/150‑1 (MQ=255)
                                                                                      ccTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTa                                                               >  1:1259906/1‑150 (MQ=255)
                                                                                       cTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTaa                                                              >  1:3798737/1‑150 (MQ=255)
                                                                                        ttCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAAc                                                             >  2:2449200/1‑150 (MQ=255)
                                                                                        ttCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAAc                                                             >  2:3385953/1‑150 (MQ=255)
                                                                                        ttCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAAc                                                             >  1:4362871/1‑150 (MQ=255)
                                                                                             gCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCa                                                        >  2:3436717/1‑150 (MQ=255)
                                                                                              cTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCAc                                                       <  1:1056148/150‑1 (MQ=255)
                                                                                              cTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCAc                                                       >  2:2118683/1‑150 (MQ=255)
                                                                                               ttGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACg                                                      >  1:4871237/1‑150 (MQ=255)
                                                                                               ttGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACg                                                      <  1:4602898/150‑1 (MQ=255)
                                                                                               ttGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACg                                                      >  2:1902853/1‑150 (MQ=255)
                                                                                                   agaagCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGgcgc                                                  >  2:4955526/1‑150 (MQ=255)
                                                                                                     aagCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCg                                                >  2:3322637/1‑150 (MQ=255)
                                                                                                          aCCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGc                                           <  2:738126/150‑1 (MQ=255)
                                                                                                          aCCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGc                                           <  2:4642375/150‑1 (MQ=255)
                                                                                                          aCCTGGGATCGTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGc                                           <  2:4128899/150‑1 (MQ=255)
                                                                                                            cTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCt                                         <  2:3128301/150‑1 (MQ=255)
                                                                                                            cTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTAGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCt                                         >  1:2179081/1‑150 (MQ=255)
                                                                                                               ggATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGAc                                      <  2:1434568/150‑1 (MQ=255)
                                                                                                               ggATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGAc                                      <  2:1434834/150‑1 (MQ=255)
                                                                                                                 aTCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACgg                                    >  2:1617772/1‑150 (MQ=255)
                                                                                                                    ttGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCa                                 >  1:3366925/1‑150 (MQ=255)
                                                                                                                     tGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCaa                                >  2:4128881/1‑150 (MQ=255)
                                                                                                                     tGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCaa                                >  1:1370051/1‑150 (MQ=255)
                                                                                                                     tGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAACCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCaa                                >  2:2564914/1‑150 (MQ=255)
                                                                                                                        cAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATAc                             >  1:1418676/1‑150 (MQ=255)
                                                                                                                           cggcgTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGc                          <  2:3034116/150‑1 (MQ=255)
                                                                                                                            ggcgTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGcc                         <  1:4497782/150‑1 (MQ=255)
                                                                                                                               gTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGt                      >  2:5120753/1‑150 (MQ=255)
                                                                                                                                tAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGtt                     >  2:1021588/1‑150 (MQ=255)
                                                                                                                                 aGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGttt                    >  1:4682833/1‑150 (MQ=255)
                                                                                                                                 aGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGttt                    <  1:4548324/150‑1 (MQ=255)
                                                                                                                                  ggTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTc                   >  1:1878394/1‑150 (MQ=255)
                                                                                                                                  ggTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTc                   >  2:4251462/1‑150 (MQ=255)
                                                                                                                                    tCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCga                 <  2:4682833/150‑1 (MQ=255)
                                                                                                                                    tCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCga                 <  2:1322581/150‑1 (MQ=255)
                                                                                                                                    tCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCga                 <  2:4877640/150‑1 (MQ=255)
                                                                                                                                      ttttCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCgaga               <  2:2754549/150‑1 (MQ=255)
                                                                                                                                       tttCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGagag              <  1:1904591/150‑1 (MQ=255)
                                                                                                                                        ttCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGt             >  1:3804401/1‑150 (MQ=255)
                                                                                                                                           cTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAACCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGAc          >  1:403461/1‑150 (MQ=255)
                                                                                                                                            tGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACg         >  2:4034133/1‑150 (MQ=255)
                                                                                                                                            tGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACg         >  2:4033972/1‑150 (MQ=255)
                                                                                                                                            tGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACg         >  2:3628013/1‑150 (MQ=255)
                                                                                                                                            tGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACg         >  2:3268875/1‑150 (MQ=255)
                                                                                                                                            tGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACg         >  2:4792362/1‑150 (MQ=255)
                                                                                                                                            tGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACg         >  2:1510946/1‑150 (MQ=255)
                                                                                                                                            tGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACg         >  1:3005202/1‑150 (MQ=255)
                                                                                                                                                  tGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACGGTCGAt   <  2:905530/150‑1 (MQ=255)
                                                                                                                                                   gAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACGGTCGAtt  <  1:50796/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCGGTTGCCTGACGGGCAATACGGCCCGTTTCGAGAGTGACGGTCGATT  >  NC_002947/5355756‑5356052

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: