Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,355,823 C→T 9.0% E88E (GAG→GAA pnp ← polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,355,8230CT 9.0% 483.3 / 27.0 166E88E (GAG→GAApnppolyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base C (73/78);  new base T (8/7);  total (81/85)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

GGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGT  >  NC_002947/5355675‑5355972
                                                                                                                                                    |                                                                                                                                                     
gggCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGtct                                                                                                                                                      <  1:37694/150‑1 (MQ=255)
 ggCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGtctc                                                                                                                                                     <  1:3568571/150‑1 (MQ=255)
 ggCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGtctc                                                                                                                                                     <  2:1400840/150‑1 (MQ=255)
  gCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTctct                                                                                                                                                    <  1:83393/150‑1 (MQ=255)
    aGCCGAGGTACAGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTAGGTCATGAAGCCATCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTctctga                                                                                                                                                  >  1:998349/1‑148 (MQ=255)
      ccGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTtct                                                                                                                                                >  1:1634218/1‑150 (MQ=255)
       cGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTtctc                                                                                                                                               <  1:401376/150‑1 (MQ=255)
         aGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCgg                                                                                                                                             >  1:1136853/1‑150 (MQ=255)
           gTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGaa                                                                                                                                           <  1:1017350/150‑1 (MQ=255)
           gTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGaa                                                                                                                                           <  2:212854/150‑1 (MQ=255)
            tACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAg                                                                                                                                          >  2:1882390/1‑150 (MQ=255)
              ccGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGa                                                                                                                                        >  2:4602898/1‑150 (MQ=255)
                 atcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACgg                                                                                                                                     >  1:3693730/1‑150 (MQ=255)
                  tcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGc                                                                                                                                    <  2:4895407/150‑1 (MQ=255)
                  tcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGc                                                                                                                                    <  2:3867582/150‑1 (MQ=255)
                   catcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGcc                                                                                                                                   <  1:1081854/150‑1 (MQ=255)
                       gCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCa                                                                                                                               >  1:3128345/1‑150 (MQ=255)
                       gCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCa                                                                                                                               >  1:3128301/1‑150 (MQ=255)
                               gTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaa                                                                                                                       <  1:1618250/150‑1 (MQ=255)
                               gTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaa                                                                                                                       >  1:1434834/1‑150 (MQ=255)
                               gTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaa                                                                                                                       >  1:1434568/1‑150 (MQ=255)
                                  cGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaagaa                                                                                                                    <  2:1647934/150‑1 (MQ=255)
                                  cGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaagaa                                                                                                                    <  2:1909683/150‑1 (MQ=255)
                                  cGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaagaa                                                                                                                    <  2:2831812/150‑1 (MQ=255)
                                  cGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaagaa                                                                                                                    <  1:2709866/150‑1 (MQ=255)
                                          tcttctTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTg                                                                                                            >  1:3034116/1‑150 (MQ=255)
                                            ttcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTggg                                                                                                          >  1:738126/1‑150 (MQ=255)
                                             tcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGa                                                                                                         <  2:3219146/150‑1 (MQ=255)
                                                 gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                                                                     <  2:799797/150‑1 (MQ=255)
                                                 gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                                                                     <  1:3459609/150‑1 (MQ=255)
                                                 gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                                                                     >  1:4211198/1‑150 (MQ=255)
                                                 gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                                                                     <  2:3500925/150‑1 (MQ=255)
                                                     gtggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCa                                                                                                 <  1:1961384/150‑1 (MQ=255)
                                                       ggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgc                                                                                               >  1:360802/1‑150 (MQ=255)
                                                       ggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgc                                                                                               >  2:1056148/1‑150 (MQ=255)
                                                        gAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgcg                                                                                              >  1:4934414/1‑150 (MQ=255)
                                                          aaCCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgcggc                                                                                            >  1:3962242/1‑150 (MQ=255)
                                                             cACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTa                                                                                         >  2:50796/1‑150 (MQ=255)
                                                             cACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTa                                                                                         >  1:1322581/1‑150 (MQ=255)
                                                              aCGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAg                                                                                        >  2:4020249/1‑150 (MQ=255)
                                                                 gTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTc                                                                                     >  2:4497782/1‑150 (MQ=255)
                                                                      gacgacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTc                                                                                <  1:3521915/150‑1 (MQ=255)
                                                                      gacgacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTc                                                                                <  2:2200317/150‑1 (MQ=255)
                                                                       acgacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTcc                                                                               <  2:2062010/150‑1 (MQ=255)
                                                                         gacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTg                                                                             <  1:1334028/150‑1 (MQ=255)
                                                                         gacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTg                                                                             <  1:4187523/150‑1 (MQ=255)
                                                                         gacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTg                                                                             >  1:905530/1‑150 (MQ=255)
                                                                          acCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTgg                                                                            <  2:2271758/150‑1 (MQ=255)
                                                                          acCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTgg                                                                            <  1:4610957/150‑1 (MQ=255)
                                                                            ctgcgtgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATa                                                                          <  2:4068849/143‑1 (MQ=255)
                                                                            cTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTa                                                                          <  1:481119/150‑1 (MQ=255)
                                                                            cTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTa                                                                          <  1:1730050/150‑1 (MQ=255)
                                                                                gtgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGttg                                                                      <  1:498736/147‑3 (MQ=255)
                                                                                 tgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCGGCGTAAGTCTTTTCCTGATAGttga                                                                     >  2:1582952/3‑147 (MQ=255)
                                                                                 cTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGaa                                                                     <  2:3455844/150‑1 (MQ=255)
                                                                                  ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGc                                                                                                                                    <  2:1148241/86‑1 (MQ=255)
                                                                                  ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                    <  2:1238668/150‑1 (MQ=255)
                                                                                  ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                    >  2:1301089/1‑150 (MQ=255)
                                                                                  tgcGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                    >  1:2754749/3‑150 (MQ=255)
                                                                                     gTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGtct                                                                                                                                                      >  1:4248626/1‑65 (MQ=255)
                                                                                         aTGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGc                                                             >  2:541295/1‑150 (MQ=255)
                                                                                           gAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCgg                                                           >  2:3214770/1‑150 (MQ=255)
                                                                                             aGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGga                                                         <  1:3348745/150‑1 (MQ=255)
                                                                                              gCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaa                                                        <  2:3693730/150‑1 (MQ=255)
                                                                                                 ttCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaagaa                                                     >  2:4775142/1‑150 (MQ=255)
                                                                                                 ttCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaagaa                                                     >  2:1425216/1‑150 (MQ=255)
                                                                                                 atcaacgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGgaagaa                                                     >  1:4199174/7‑150 (MQ=255)
                                                                                                    cGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc                                                  <  2:3342728/150‑1 (MQ=255)
                                                                                                    cGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc                                                  <  2:4019775/150‑1 (MQ=255)
                                                                                                    aacgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGaaatc                                                  >  1:1860273/4‑148 (MQ=255)
                                                                                                     acgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAAtcc                                                 >  2:1573333/3‑150 (MQ=255)
                                                                                                         aaaGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAATCCtgac                                             <  1:413123/148‑4 (MQ=255)
                                                                                                         aCAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGc                                             <  1:3740146/150‑1 (MQ=255)
                                                                                                           aGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTATCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCt                                           >  2:2932558/1‑150 (MQ=255)
                                                                                                           aGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGCCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCGTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCt                                           <  1:4497803/150‑1 (MQ=255)
                                                                                                             cggacggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTgg                                         >  1:4676008/1‑150 (MQ=255)
                                                                                                               gacggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGACGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGAt                                       <  2:2055718/150‑1 (MQ=255)
                                                                                                                acggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATc                                      >  2:4543540/1‑150 (MQ=255)
                                                                                                                  ggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAg                                    >  2:4388475/1‑150 (MQ=255)
                                                                                                                   gaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGc                                   <  1:2885130/150‑1 (MQ=255)
                                                                                                                    aTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGcc                                  >  2:2426543/1‑150 (MQ=255)
                                                                                                                        gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGtt                              <  1:216626/150‑1 (MQ=255)
                                                                                                                        gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGtt                              >  1:3719564/1‑150 (MQ=255)
                                                                                                                        gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGttt                             <  2:336283/150‑1 (MQ=255)
                                                                                                                             tCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGgcg                         >  1:806686/1‑150 (MQ=255)
                                                                                                                               gATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGcc                       >  1:2613649/1‑150 (MQ=255)
                                                                                                                               gATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGcc                       >  1:170731/1‑150 (MQ=255)
                                                                                                                                      cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                <  1:2316541/150‑1 (MQ=255)
                                                                                                                                      cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                <  1:2316468/150‑1 (MQ=255)
                                                                                                                                      cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                <  1:630909/150‑1 (MQ=255)
                                                                                                                                      cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                <  1:3831709/150‑1 (MQ=255)
                                                                                                                                      cgcgAGGTCAGCGTCTCGTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                <  1:2737410/150‑1 (MQ=255)
                                                                                                                                         gAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTc             >  1:1024814/1‑150 (MQ=255)
                                                                                                                                            gTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCAcc          >  1:4116667/1‑150 (MQ=255)
                                                                                                                                             tCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCa         <  1:1882390/150‑1 (MQ=255)
                                                                                                                                             tCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCa         <  1:40444/150‑1 (MQ=255)
                                                                                                                                             tCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCa         >  2:4481867/1‑150 (MQ=255)
                                                                                                                                               aGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGc       >  1:1985135/1‑150 (MQ=255)
                                                                                                                                                gCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCa      >  2:4326571/1‑150 (MQ=255)
                                                                                                                                                    ctctTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGt  <  1:2584887/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGT  >  NC_002947/5355675‑5355972

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: