Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,355,781 C→A 7.7% L102L (CTG→CTT pnp ← polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,355,7810CA 7.7% 536.3 / 22.3 164L102L (CTG→CTTpnppolyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base C (67/84);  new base A (7/6);  total (74/90)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

GCGCCGATCGGGCCTTCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGC  >  NC_002947/5355635‑5355929
                                                                                                                                                  |                                                                                                                                                    
  gcCGATCGGGCCTTCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGcgg                                                                                                                                                 >  1:2154611/1‑150 (MQ=255)
  gcCGATCGGGCCTTCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGcgg                                                                                                                                                 >  2:1066600/1‑150 (MQ=255)
   ccGATCGGGCCTTCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGcgga                                                                                                                                                <  1:1856853/150‑1 (MQ=255)
         gggCCTTCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATc                                                                                                                                          <  1:2121923/150‑1 (MQ=255)
            ccTTCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGa                                                                                                                                       >  2:2316541/1‑150 (MQ=255)
            ccTTCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGa                                                                                                                                       >  1:4022278/1‑150 (MQ=255)
            ccTTCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGa                                                                                                                                       <  1:1940595/150‑1 (MQ=255)
              ttCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACg                                                                                                                                     >  2:2885130/1‑150 (MQ=255)
              ttCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACg                                                                                                                                     <  2:2055112/150‑1 (MQ=255)
                cGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGt                                                                                                                                   <  2:1955639/150‑1 (MQ=255)
                 gAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTc                                                                                                                                  >  2:2737410/1‑150 (MQ=255)
                 gAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTc                                                                                                                                  <  1:3808925/150‑1 (MQ=255)
                   aCGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGa                                                                                                                                >  2:3936762/1‑150 (MQ=255)
                      gAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCa                                                                                                                             <  1:3796224/150‑1 (MQ=255)
                      gAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCAATCa                                                                                                                             <  2:2146811/150‑1 (MQ=255)
                         tGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGTTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGgc                                                                                                                          >  2:2584887/1‑150 (MQ=255)
                            cGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGcgcg                                                                                                                       >  1:2349912/1‑150 (MQ=255)
                              gAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAg                                                                                                                     <  2:424369/150‑1 (MQ=255)
                               aaaTGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAgg                                                                                                                    <  2:3867296/150‑1 (MQ=255)
                                aaTGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGt                                                                                                                   <  2:1381781/150‑1 (MQ=255)
                                    gCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGc                                                                                                               <  2:4105608/150‑1 (MQ=255)
                                    gCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGc                                                                                                               <  1:1550307/150‑1 (MQ=255)
                                    gCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGc                                                                                                               <  1:1924718/150‑1 (MQ=255)
                                    gCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGc                                                                                                               <  2:1534254/150‑1 (MQ=255)
                                     ccAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCg                                                                                                              >  2:40444/1‑150 (MQ=255)
                                              ccGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTtct                                                                                                     >  1:1634218/1‑150 (MQ=255)
                                                 aGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCgg                                                                                                  >  1:1136853/1‑150 (MQ=255)
                                                   gTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGaa                                                                                                <  1:1017350/150‑1 (MQ=255)
                                                   gTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGaa                                                                                                <  2:212854/150‑1 (MQ=255)
                                                   gTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGaa                                                                                                <  1:182079/150‑1 (MQ=255)
                                                   gTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCGTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGaa                                                                                                <  2:4504017/150‑1 (MQ=255)
                                                    tACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAg                                                                                               >  2:1882390/1‑150 (MQ=255)
                                                        gATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACg                                                                                           <  1:1989436/150‑1 (MQ=255)
                                                         atcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACgg                                                                                          >  1:3693730/1‑150 (MQ=255)
                                                          tcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGc                                                                                         <  2:4895407/150‑1 (MQ=255)
                                                          tcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGc                                                                                         >  1:1417326/1‑150 (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:1434834/1‑150 (MQ=255)
                                                                       gTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaa                                                                            >  1:1434568/1‑150 (MQ=255)
                                                                          cGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaagaa                                                                         <  2:1647934/150‑1 (MQ=255)
                                                                          cGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaagaa                                                                         <  2:1417326/150‑1 (MQ=255)
                                                                          cGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaagaa                                                                         <  1:2709866/150‑1 (MQ=255)
                                                                                  tcttctTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTg                                                                 >  1:3034116/1‑150 (MQ=255)
                                                                                    ttcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTggg                                                               >  2:572996/1‑150 (MQ=255)
                                                                                    ttcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTggg                                                               >  1:738126/1‑150 (MQ=255)
                                                                                     tcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGa                                                              <  1:3718774/150‑1 (MQ=255)
                                                                                     tcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGa                                                              <  2:3219146/150‑1 (MQ=255)
                                                                                      cttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGAt                                                             <  1:4205614/150‑1 (MQ=255)
                                                                                         gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                          <  2:799797/150‑1 (MQ=255)
                                                                                         gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                          <  2:462357/150‑1 (MQ=255)
                                                                                         gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                          <  1:3459609/150‑1 (MQ=255)
                                                                                         gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                          >  1:4211198/1‑150 (MQ=255)
                                                                                          cTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTg                                                         >  1:4642375/1‑150 (MQ=255)
                                                                                          cTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTg                                                         >  1:4128899/1‑150 (MQ=255)
                                                                                               ggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgc                                                    >  2:1056148/1‑150 (MQ=255)
                                                                                               ggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgc                                                    >  1:360802/1‑150 (MQ=255)
                                                                                                gAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgcg                                                   >  1:4934414/1‑150 (MQ=255)
                                                                                                  aaCCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgcggc                                                 >  1:3962242/1‑150 (MQ=255)
                                                                                                     cACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTa                                              >  1:1322581/1‑150 (MQ=255)
                                                                                                     cACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTa                                              >  2:50796/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)
                                                                                                               acgacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTcc                                    <  2:2062010/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                               <  1:2971208/143‑1 (MQ=255)
                                                                                                                    ctgcgtgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATa                               <  2:3807330/143‑1 (MQ=255)
                                                                                                                    cTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTa                               <  1:1730050/150‑1 (MQ=255)
                                                                                                                    cTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTa                               <  1:3976441/150‑1 (MQ=255)
                                                                                                                    cTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTa                               <  1:481119/150‑1 (MQ=255)
                                                                                                                     tGCACTTCGTTCATGAAGCCATCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGACAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGACACCTGGGATCTTGCCAGCGGCGAAGGTCTTTTCCTGGTAg                              >  2:4548324/1‑150 (MQ=255)
                                                                                                                        gtgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGttg                           <  1:498736/147‑3 (MQ=255)
                                                                                                                         tgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGttga                          <  1:1710888/148‑4 (MQ=255)
                                                                                                                         tgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCGGCGTAAGTCTTTTCCTGATAGttga                          >  2:1582952/3‑147 (MQ=255)
                                                                                                                         cTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGaa                          <  2:3455844/150‑1 (MQ=255)
                                                                                                                          ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGc                                                                                         >  1:1148241/1‑86 (MQ=255)
                                                                                                                          ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                         <  2:1238668/150‑1 (MQ=255)
                                                                                                                          ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                         >  2:1301089/1‑150 (MQ=255)
                                                                                                                          tgcGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                         >  1:2754749/3‑150 (MQ=255)
                                                                                                                               tCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACa                    >  1:1865252/1‑150 (MQ=255)
                                                                                                                                 aTGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGc                  >  1:4892508/1‑150 (MQ=255)
                                                                                                                                 aTGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGc                  >  2:541295/1‑150 (MQ=255)
                                                                                                                                 aTGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTACTGGTAGTGAACCGACAGc                  >  2:1429694/1‑150 (MQ=255)
                                                                                                                                   gAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCgg                >  2:3214770/1‑150 (MQ=255)
                                                                                                                                     aGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGga              <  1:3348745/150‑1 (MQ=255)
                                                                                                                                      gCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCAc                                                                                   >  2:3378489/1‑80 (MQ=255)
                                                                                                                                      gCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaa             <  2:3693730/150‑1 (MQ=255)
                                                                                                                                         ttCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaagaa          >  1:408538/1‑150 (MQ=255)
                                                                                                                                         ttCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaagaa          >  2:4775142/1‑150 (MQ=255)
                                                                                                                                         atcaacgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGgaagaa          >  1:2691825/7‑150 (MQ=255)
                                                                                                                                           caacgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGaaat        <  2:2932017/146‑2 (MQ=255)
                                                                                                                                            cGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc       <  2:3342728/150‑1 (MQ=255)
                                                                                                                                               gAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCtt    <  1:446266/150‑1 (MQ=255)
                                                                                                                                               gAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAATCCtg    >  1:1609826/1‑149 (MQ=255)
                                                                                                                                                 aCAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGc  <  1:3742750/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGCCGATCGGGCCTTCGAACGGAATGCCGGAAATGGCCAGGGCAGCCGAGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGC  >  NC_002947/5355635‑5355929

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: