Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,755,404 A→G 6.9% H388H (CAT→CAC thiI ← tRNA sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,755,4040AG 6.9% 855.2 / 33.9 275H388H (CAT→CACthiItRNA sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (135/120);  new base G (12/7);  total (147/128)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

Only 100 of 287 total aligned reads displayed.

AGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGG  >  NC_002947/5755256‑5755552
                                                                                                                                                    |                                                                                                                                                    
  aaCGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGa                                                                                                                                                   <  2:908843/150‑1 (MQ=255)
    cGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATc                                                                                                                                                 <  2:3065874/150‑1 (MQ=255)
         tGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGATga                                                                                                                                            <  2:427517/150‑2 (MQ=255)
         tGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGATga                                                                                                                                            <  2:2387160/150‑2 (MQ=255)
         tGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATgg                                                                                                                                            >  1:4923842/1‑150 (MQ=255)
         tGGCCTGGACGTCGAGACCGTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCGGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCATGTGATCGATga                                                                                                                                            <  2:5218821/150‑2 (MQ=255)
            ccTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGaaa                                                                                                                                         <  2:4004918/150‑1 (MQ=255)
            ccTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGaaa                                                                                                                                         <  2:1893421/150‑1 (MQ=255)
             cTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAAt                                                                                                                                        >  1:1605611/1‑150 (MQ=255)
             cTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGAGGGAAAt                                                                                                                                        >  2:393004/1‑150 (MQ=255)
               ggACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGa                                                                                                                                      >  1:2966362/1‑150 (MQ=255)
                   ttCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGACATGAACtt                                                                                                                                  >  2:5419868/1‑150 (MQ=255)
                    tCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTg                                                                                                                                 <  1:4441783/150‑1 (MQ=255)
                    ccGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTg                                                                                                                                 <  1:1073204/149‑1 (MQ=255)
                      gATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGgc                                                                                                                               >  1:5527286/1‑150 (MQ=255)
                      gATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGgc                                                                                                                               >  2:5334719/1‑150 (MQ=255)
                      gATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGgc                                                                                                                               <  2:267141/150‑1 (MQ=255)
                                    ccAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgc                                                                                                                 >  1:3475796/1‑150 (MQ=255)
                                    ccAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgc                                                                                                                 <  2:3960813/150‑1 (MQ=255)
                                     cAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgcg                                                                                                                >  2:2032852/1‑150 (MQ=255)
                                      agagGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgcgc                                                                                                               >  1:1867516/1‑150 (MQ=255)
                                        agGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTc                                                                                                             <  1:5384768/150‑1 (MQ=255)
                                         ggTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgccaa                                                                                                            <  1:1229361/150‑4 (MQ=255)
                                         ggTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTAGCGCGCTCTTGAgcgcgccaa                                                                                                            <  1:5277172/150‑4 (MQ=255)
                                          gTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGa                                                                                                                                      <  1:3624523/123‑1 (MQ=255)
                                          gTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGa                                                                                                                                      >  2:3624523/1‑123 (MQ=255)
                                          gTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCa                                                                                                           >  1:2319/1‑150 (MQ=255)
                                          gTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCa                                                                                                           >  2:1303074/1‑150 (MQ=255)
                                          gTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCa                                                                                                           >  1:2160115/1‑150 (MQ=255)
                                            tCGTCTTCCTGGGCATGGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgccaagtt                                                                                                         <  1:7758/150‑7 (MQ=255)
                                            tCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgccaagtt                                                                                                         <  1:487871/150‑7 (MQ=255)
                                               tCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATgg                                                                                                                                            >  2:1937365/1‑112 (MQ=255)
                                                cTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACgg                                                                                                     <  2:2841524/150‑1 (MQ=255)
                                                  tCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGcc                                                                                                   <  2:1973921/150‑1 (MQ=255)
                                                  tCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGcc                                                                                                   <  2:4604367/150‑1 (MQ=255)
                                                   ccTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCa                                                                                                  <  1:3878322/150‑1 (MQ=255)
                                                     tGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATg                                                                                                >  1:3496879/1‑150 (MQ=255)
                                                     tGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATg                                                                                                >  1:5602464/1‑150 (MQ=255)
                                                      gggCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGt                                                                                               <  2:1007388/150‑1 (MQ=255)
                                                       ggCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTc                                                                                              <  1:4351218/150‑1 (MQ=255)
                                                          aTCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGaa                                                                                           >  2:2155840/1‑150 (MQ=255)
                                                          aTCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGaa                                                                                           >  2:4786906/1‑150 (MQ=255)
                                                          aTCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGaa                                                                                           <  1:2405794/150‑1 (MQ=255)
                                                              gggTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGc                                                                                       <  1:2402916/150‑1 (MQ=255)
                                                                   gCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttc                                                                                  <  2:294360/150‑1 (MQ=255)
                                                                   gCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttc                                                                                  <  1:469203/150‑1 (MQ=255)
                                                                     gAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttctt                                                                                >  2:1139158/1‑150 (MQ=255)
                                                                     gAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttctt                                                                                >  2:546511/1‑150 (MQ=255)
                                                                      aaTGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttcttc                                                                               <  2:5032892/150‑1 (MQ=255)
                                                                      aaTGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttcttc                                                                               <  2:2525419/150‑1 (MQ=255)
                                                                          tCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGc                                                                           <  1:3724917/150‑1 (MQ=255)
                                                                             aTGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTcc                                                                        >  2:24227/1‑150 (MQ=255)
                                                                                 cGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGc                                                                    <  2:3500146/150‑1 (MQ=255)
                                                                                  gATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGc                                                                                                    <  1:3844291/117‑1 (MQ=255)
                                                                                  gATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGc                                                                                                    >  2:3844291/1‑117 (MQ=255)
                                                                                  gATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCg                                                                   >  1:220137/1‑150 (MQ=255)
                                                                                    tCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCTATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGt                                                                 <  2:1520111/150‑1 (MQ=255)
                                                                                      tggcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGa                                                               >  2:5496004/1‑150 (MQ=255)
                                                                                       ggcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGAc                                                              >  1:209630/1‑150 (MQ=255)
                                                                                           tggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCtt                                                          >  1:3220232/1‑150 (MQ=255)
                                                                                              cAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGc                                                       >  1:4781573/1‑150 (MQ=255)
                                                                                                gcgcCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCAt                                                     >  1:5080957/1‑150 (MQ=255)
                                                                                                gcgcCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCAt                                                     <  2:3410414/150‑1 (MQ=255)
                                                                                                   ccTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCAgggg                                                  >  2:5080160/1‑150 (MQ=255)
                                                                                                    cTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggt                                                 <  2:5113100/150‑1 (MQ=255)
                                                                                                    cTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggt                                                 <  2:4680896/150‑1 (MQ=255)
                                                                                                    cTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggt                                                 <  1:1962166/150‑1 (MQ=255)
                                                                                                     tCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggtg                                                <  2:3319424/150‑1 (MQ=255)
                                                                                                     tCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGggtg                                                >  2:2821116/1‑150 (MQ=255)
                                                                                                        gCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTc                                             <  2:5296059/150‑1 (MQ=255)
                                                                                                          cACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCgg                                           <  1:5217867/150‑1 (MQ=255)
                                                                                                           aCTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCggg                                          >  1:685551/1‑150 (MQ=255)
                                                                                                           aCTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCggg                                          >  1:14072/1‑150 (MQ=255)
                                                                                                              tCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGttt                                       >  1:3680767/1‑150 (MQ=255)
                                                                                                                cTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAc                                     >  2:2552704/1‑150 (MQ=255)
                                                                                                                 tCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAcc                                    >  1:2782078/1‑150 (MQ=255)
                                                                                                                 tCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAcc                                    >  1:1023674/1‑150 (MQ=255)
                                                                                                                 tCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAcc                                    >  2:4112734/1‑150 (MQ=255)
                                                                                                                 tCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAcc                                    >  2:322127/1‑150 (MQ=255)
                                                                                                                   gATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCga                                  <  2:1634791/150‑1 (MQ=255)
                                                                                                                   gATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCga                                  <  2:4925307/150‑1 (MQ=255)
                                                                                                                         gATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCAc                            >  1:1726197/1‑150 (MQ=255)
                                                                                                                         gATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCAc                            >  1:3966610/1‑150 (MQ=255)
                                                                                                                          aTGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCAcg                           <  2:4698223/150‑1 (MQ=255)
                                                                                                                           tGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCAcgc                          >  2:3880697/1‑150 (MQ=255)
                                                                                                                                ttGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAg                     >  1:959262/1‑150 (MQ=255)
                                                                                                                                    ccAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATggg                                                  >  2:1851307/1‑117 (MQ=255)
                                                                                                                                     cAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTc                >  1:3354684/1‑150 (MQ=255)
                                                                                                                                       gCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCg              >  1:1294387/1‑150 (MQ=255)
                                                                                                                                          cATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCa           <  2:2966362/150‑1 (MQ=255)
                                                                                                                                          cATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCa           <  2:2288837/150‑1 (MQ=255)
                                                                                                                                           aTCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt          >  1:1702667/1‑150 (MQ=255)
                                                                                                                                           aTCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt          <  2:61371/150‑1 (MQ=255)
                                                                                                                                           aTCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt          >  1:1524252/1‑150 (MQ=255)
                                                                                                                                               aTCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGc      <  1:361585/150‑1 (MQ=255)
                                                                                                                                               aTCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGc      >  2:4372522/1‑150 (MQ=255)
                                                                                                                                                  acaTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTg   <  2:727614/150‑1 (MQ=255)
                                                                                                                                                   caTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTgg  <  2:3701402/150‑1 (MQ=255)
                                                                                                                                                   caTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTgg  >  2:1775100/1‑150 (MQ=255)
                                                                                                                                                   caTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTgg  >  1:5121194/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGAACGGCATGGCCTGGACTTCGATACCTTCCAGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGG  >  NC_002947/5755256‑5755552

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: