Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 961,717 C→T 9.3% G234D (GGC→GAC)  selD ← selenide, water dikinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947961,7170CT 9.3% 777.3 / 50.0 268G234D (GGC→GAC) selDselenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base C (103/140);  new base T (13/12);  total (116/152)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 273 total aligned reads displayed.

TCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGT  >  NC_002947/961568‑961865
                                                                                                                                                     |                                                                                                                                                    
tCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTACTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGc                                                                                                                                                      >  1:677336/1‑150 (MQ=255)
 cTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGcc                                                                                                                                                     <  1:3404456/150‑1 (MQ=255)
   tgtgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGa                                                                                                                                                   <  1:3794721/150‑1 (MQ=255)
   tgtgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACGAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGa                                                                                                                                                   <  2:4830197/150‑1 (MQ=255)
    gtgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGGCCACCTGGTCGAg                                                                                                                                                  <  1:5093315/150‑1 (MQ=255)
     tgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGc                                                                                                                                                 >  2:2965315/1‑150 (MQ=255)
        ccGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCTCTGCCTTCGGCCACCTGGTCGAGCTGGCCGaaggca                                                                                                                                              <  1:9005822/150‑6 (MQ=255)
         cGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCTGCTCGGCCACCTGGTCGAGCTGGc                                                                                                                                             <  2:6814322/150‑1 (MQ=255)
              cTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAg                                                                                                                                        >  1:4008153/1‑150 (MQ=255)
               tGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCTGCTCGGCCACCTGGTCGAGCTGGCCGAAgg                                                                                                                                       <  2:7954800/150‑2 (MQ=39)
                gTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGcc                                                                                                                                      >  1:3025300/1‑150 (MQ=255)
                 tCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCg                                                                                                                                     >  1:8379701/1‑150 (MQ=255)
                  cGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCgg                                                                                                                                    <  1:7240437/150‑1 (MQ=255)
                  cGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCgg                                                                                                                                    <  2:5976380/150‑1 (MQ=255)
                   gTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGaa                                                                                                                                         <  1:1385844/144‑1 (MQ=255)
                   gTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGt                                                                                                                                   >  2:2316164/1‑150 (MQ=255)
                      gTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGAc                                                                                                                                >  2:805752/1‑150 (MQ=255)
                         gcgcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTc                                                                                                                             <  1:1574906/150‑1 (MQ=255)
                         gcgcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTc                                                                                                                             >  2:7497793/1‑150 (MQ=255)
                            cAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGt                                                                                                                          >  1:2958986/1‑150 (MQ=255)
                                gTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                                      <  1:805752/150‑1 (MQ=255)
                                gTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAACTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                                      <  2:6458527/150‑1 (MQ=255)
                                 ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATgg                                                                                                                     <  2:5246776/150‑1 (MQ=255)
                                 ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATgg                                                                                                                     <  1:5515071/150‑1 (MQ=255)
                                  tCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGc                                                                                                                    <  1:6810337/150‑1 (MQ=255)
                                     aCCGGGGATGCAGCCTTCGGCCAGGTCATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCACGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCtt                                                                                                                 >  1:945139/1‑150 (MQ=255)
                                      ccGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTg                                                                                                                >  1:5702100/1‑150 (MQ=255)
                                       cGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGa                                                                                                               <  1:7666232/150‑1 (MQ=255)
                                        ggggATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAc                                                                                                              >  1:8541260/1‑150 (MQ=255)
                                           gATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAcccc                                                                                                           <  1:8637082/150‑1 (MQ=255)
                                           gATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAcccc                                                                                                           >  2:8700819/1‑150 (MQ=255)
                                             tGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGt                                                                                                         >  2:7259671/1‑150 (MQ=255)
                                             tGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGt                                                                                                         >  1:6179568/1‑150 (MQ=255)
                                                aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCa                                                                                                      <  1:2544402/150‑1 (MQ=255)
                                                aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTAGACCCCGACCa                                                                                                      >  2:6088223/1‑150 (MQ=255)
                                                    ttCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGttt                                                                                                  <  2:9021483/150‑1 (MQ=255)
                                                    ttCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGttt                                                                                                  <  1:7266292/150‑1 (MQ=255)
                                                    ttCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGttt                                                                                                  <  1:4995241/150‑1 (MQ=255)
                                                      cGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGc                                                                                                <  2:341303/150‑1 (MQ=255)
                                                      cGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGc                                                                                                <  2:7828095/150‑1 (MQ=255)
                                                        gCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCg                                                                                              >  1:4949209/1‑150 (MQ=255)
                                                          cAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGaa                                                                                            >  2:4512500/1‑150 (MQ=255)
                                                          cAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGaa                                                                                            <  1:3087758/150‑1 (MQ=255)
                                                           aGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCTGCTCGGCCACCTGGTCGAGCTGGCCGaa                                                                                                                                         >  2:1560929/1‑104 (MQ=37)
                                                           aGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAg                                                                                           >  2:2752753/1‑150 (MQ=255)
                                                            ggTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGc                                                                                          >  1:5312604/1‑150 (MQ=255)
                                                                 tGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTg                                                                                     >  1:4616461/1‑150 (MQ=255)
                                                                   gTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                                   >  2:3577319/1‑150 (MQ=255)
                                                                      gACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCggg                                                                                >  1:3294054/1‑150 (MQ=255)
                                                                       acacTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGc                                                                               <  2:333426/150‑1 (MQ=255)
                                                                         acTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGTCGTCGGTCAAGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgt                                                                             >  1:1089079/1‑150 (MQ=255)
                                                                         acTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgt                                                                             <  2:1873664/150‑1 (MQ=255)
                                                                          cTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgtg                                                                            <  2:4233335/150‑1 (MQ=255)
                                                                           tAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgtgt                                                                           >  1:2302354/1‑150 (MQ=255)
                                                                             ggcaggcGTGGCTCGGCGGCGTAGTTCAGGTGCGCTGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCAGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTg                                                                         <  2:3729471/150‑1 (MQ=255)
                                                                               caggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCTGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAg                                                                       <  1:5904389/150‑1 (MQ=255)
                                                                               caggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAg                                                                       >  1:1710312/1‑150 (MQ=255)
                                                                                 ggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAggg                                                                     <  1:8562925/150‑1 (MQ=255)
                                                                                  gcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGt                                                                    <  2:6116283/150‑1 (MQ=255)
                                                                                    gTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGt                                                                    <  1:6566299/148‑1 (MQ=255)
                                                                                    gTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGc                                                                  <  2:3759255/150‑1 (MQ=255)
                                                                                    gTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCGGCCGGGCGTGTTGAGGGTGc                                                                  >  1:7109764/1‑150 (MQ=255)
                                                                                     tGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAgg                                                                      <  1:2135505/145‑1 (MQ=255)
                                                                                              ggcgTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTc                                                        <  1:3499828/150‑1 (MQ=255)
                                                                                                 gTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTcgcg                                                     <  2:8210674/150‑1 (MQ=255)
                                                                                                 gTAGTTCAGGTGCGCGGAGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTcgcg                                                     >  1:7392466/1‑150 (MQ=255)
                                                                                                  tAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGa                                                    <  1:1037315/150‑1 (MQ=255)
                                                                                                    gTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGc                                                  <  1:3891696/150‑1 (MQ=255)
                                                                                                      tCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCa                                                <  2:1032863/150‑1 (MQ=255)
                                                                                                       cAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAg                                               <  1:1002792/150‑1 (MQ=255)
                                                                                                             cgcgGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTg                                         <  2:7903851/150‑1 (MQ=255)
                                                                                                             cgcgGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTg                                         <  2:3680263/150‑1 (MQ=255)
                                                                                                              gcgGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTgg                                        <  1:3238626/150‑1 (MQ=255)
                                                                                                                  tGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCt                                    <  2:2997830/150‑1 (MQ=255)
                                                                                                                  tGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCt                                    <  2:8528951/150‑1 (MQ=255)
                                                                                                                  tGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCt                                    <  2:8964373/150‑1 (MQ=255)
                                                                                                                  tGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCt                                    <  2:2517311/150‑1 (MQ=255)
                                                                                                                   gAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTg                                   <  2:7192064/150‑1 (MQ=255)
                                                                                                                      gCCGCTGCCTTCGGCCAGCACGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCAAGGCCATGACCCCGGCCAGTTTGCCGAAGCGGCTGCCCGGCGTGTTTAGGGTGCACAACCACTCGCGAGCCAGGCCTTGGTCCTGcac                                >  1:4110891/1‑150 (MQ=255)
                                                                                                                         gCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCCGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACg                             >  2:411114/1‑150 (MQ=255)
                                                                                                                         accggcTTCGGCCTGCTCGGCCACCTGGTCGAGCTGGCCGAAGGCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACg                             >  2:4643870/6‑150 (MQ=37)
                                                                                                                               ttCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCtt                       >  1:1288679/1‑150 (MQ=255)
                                                                                                                               ttCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCtt                       >  2:1457630/1‑150 (MQ=255)
                                                                                                                                tCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTg                      <  2:6691474/150‑1 (MQ=255)
                                                                                                                                 cGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTgg                     >  1:8770179/1‑150 (MQ=255)
                                                                                                                                  ggCCTGCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGc                    >  1:4867332/1‑150 (MQ=255)
                                                                                                                                  ggCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGc                    <  1:8841630/150‑1 (MQ=255)
                                                                                                                                   gCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGcc                   >  1:7914462/1‑150 (MQ=255)
                                                                                                                                    cctgCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCct                  >  2:2576182/4‑150 (MQ=255)
                                                                                                                                    cctgCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCct                  >  1:8317175/4‑150 (MQ=255)
                                                                                                                                      tgCTCGGCCACCTGGTCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCcttc                >  2:6698794/2‑150 (MQ=255)
                                                                                                                                           gACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTc           <  2:4357390/150‑1 (MQ=255)
                                                                                                                                            aCCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCg          >  1:7821361/1‑150 (MQ=255)
                                                                                                                                            aCCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCg          >  1:8262418/1‑150 (MQ=255)
                                                                                                                                             ccacctGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCgg         >  1:2202314/6‑150 (MQ=255)
                                                                                                                                             ccAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCgg         >  2:5596687/1‑150 (MQ=255)
                                                                                                                                               aGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCa       >  1:814807/1‑150 (MQ=255)
                                                                                                                                                 gagtccGAGCAGGCCGAAGCCGGTGACGTCGGACATGGCCTTGACCCCGTCCAGTTAGCCGAAGCGGCTGCCGGGCGTGTTGAGGGAGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCGTCTTTTCGGCAgt     >  2:2187299/5‑150 (MQ=255)
                                                                                                                                                 gTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgt     >  2:2467604/1‑150 (MQ=255)
                                                                                                                                                   ggCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgtgg   <  2:1303576/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGT  >  NC_002947/961568‑961865

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑