Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,257,737 A→C 5.3% V38G (GTG→GGG)  PP_4634 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,257,7370AC 5.3% 292.4 / 5.8 103V38G (GTG→GGG) PP_4634alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base A (46/50);  new base C (3/3);  total (49/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 103 total aligned reads displayed.

GGAGGCCTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAA  >  NC_002947/5257589‑5257885
                                                                                                                                                     |                                                                                                                                                    
ggaggCCTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCAc                                                                                                                                                     <  2:556211/150‑1 (MQ=255)
     ccTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCg                                                                                                                                                >  2:1903798/1‑150 (MQ=255)
      cTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCgg                                                                                                                                               >  2:1086156/1‑150 (MQ=255)
      cTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCgg                                                                                                                                               >  1:405096/1‑150 (MQ=255)
      cTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCgg                                                                                                                                               >  1:203724/1‑150 (MQ=255)
      cTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCgg                                                                                                                                               >  1:810007/1‑150 (MQ=255)
       tGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCggg                                                                                                                                              >  1:1621167/1‑150 (MQ=255)
        ggAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGgc                                                                                                                                             >  1:750521/1‑150 (MQ=255)
          aaaCCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGgcgc                                                                                                                                           <  1:1874331/150‑1 (MQ=255)
             ccACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTgg                                                                                                                                        <  1:1089854/150‑1 (MQ=255)
              cACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTggg                                                                                                                                       <  1:1477430/150‑1 (MQ=255)
              cACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTggg                                                                                                                                       <  1:294202/150‑1 (MQ=255)
              cACGCAGCGAAGCCGGCAGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTggg                                                                                                                                       >  2:168513/1‑150 (MQ=255)
               aCGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTgggg                                                                                                                                      >  2:745385/1‑150 (MQ=255)
                cGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCCCCAGCGCCCGCTGGTGa                                                                                                                                     <  2:478088/150‑1 (MQ=255)
                   aGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCa                                                                                                                                  <  1:1581946/150‑1 (MQ=255)
                     cGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCCCCAGCGCCCGCTGGTGATCatc                                                                                                                                <  2:2608/150‑3 (MQ=255)
                      gAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGt                                                                                                                               >  2:1728011/1‑150 (MQ=255)
                       aaGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGtt                                                                                                                              <  2:78265/150‑1 (MQ=255)
                       aaGCCGGCCGGCATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCAGAACAAAACTGGTTATTCGGACCACTGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGtt                                                                                                                              <  1:134432/150‑1 (MQ=255)
                         gccggccgGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTgg                                                                                                                            <  1:745385/150‑1 (MQ=255)
                         gccggccgGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTgg                                                                                                                            <  2:1761379/150‑1 (MQ=255)
                         gccggccgGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAAGTGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTgg                                                                                                                            <  2:457319/150‑1 (MQ=255)
                          ccggccggGATTCGCCGCTGCCACGGTAGTCGTAAGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGGTCACCAGCGGGCGATGGGGATCAGGTTgga                                                                                                                           <  2:27299/150‑2 (MQ=255)
                           cggccggGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGc                                                                                                                          <  1:1903798/150‑1 (MQ=255)
                             gccggGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGt                                                                                                                        <  1:933822/150‑1 (MQ=255)
                                gggATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCCCCAGCGCCCGCTg                                                                                                                                         >  2:1356048/1‑130 (MQ=255)
                                gggATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATTAGGTTGGGCGTGgc                                                                                                                     <  2:1504669/150‑1 (MQ=255)
                                    ttCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCca                                                                                                                 <  2:750521/150‑1 (MQ=255)
                                    ttCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCca                                                                                                                 <  1:1086156/150‑1 (MQ=255)
                                          gATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGTTTGGGCGTGGCGCCAGCAGaa                                                                                                           <  1:335672/150‑1 (MQ=255)
                                          gATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGaa                                                                                                           <  2:489119/150‑1 (MQ=255)
                                          gATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGaa                                                                                                           <  2:365703/150‑1 (MQ=255)
                                          gATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGAGGCAGCGTACCGAAGTCGTAGCCATGATGATCATCAGCGGGCGCTGGGGATCAGGTTTGGCGTGGCGCCAGCAGaa                                                                                                           <  1:1674244/150‑1 (MQ=255)
                                             gCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGAGGGCGAATAGGTAGTTGACGAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAgcc                                                                                                        >  1:712626/1‑150 (MQ=255)
                                              ccACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAgccg                                                                                                       <  2:238145/150‑1 (MQ=255)
                                               cACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAgccgc                                                                                                      <  1:1263773/150‑1 (MQ=255)
                                               cACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAgccgc                                                                                                      <  2:1881854/150‑1 (MQ=255)
                                               cACGGTAGTCGTAGGTCAGCACATCGCAGCCGAGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGAGGCAGCGTTGATGATCCCCAGCGGGCGCTGGGGTTCAGGTTGGGCGTGGCGCTAGCAGAAgccgc                                                                                                      >  2:395777/1‑150 (MQ=255)
                                                aCGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGc                                                                                                                                                 >  2:1823440/1‑106 (MQ=255)
                                                aCGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGc                                                                                                                                                 <  1:1823440/106‑1 (MQ=255)
                                                aCGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAgccgcc                                                                                                     >  2:5955/1‑150 (MQ=255)
                                                aCGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAgccgcc                                                                                                     >  1:139907/1‑150 (MQ=255)
                                                  ggTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGa                                                                                                   <  2:712626/150‑1 (MQ=255)
                                                    tAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAgg                                                                                                 <  1:1513336/150‑1 (MQ=255)
                                                      gTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGTGTGGCGCTAGCAGAAGCCGCCGAGGCt                                                                                               <  1:1825302/150‑1 (MQ=255)
                                                      gTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCt                                                                                               <  2:1621167/150‑1 (MQ=255)
                                                         gTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATc                                                                                                                                   >  2:728074/1‑111 (MQ=255)
                                                         gTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATc                                                                                                                                   <  1:728074/111‑1 (MQ=255)
                                                         gTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTa                                                                                            <  1:168513/150‑1 (MQ=255)
                                                             gTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCg                                                                                        >  2:1206277/1‑150 (MQ=255)
                                                              tCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGt                                                                                       >  2:617876/1‑150 (MQ=255)
                                                              tCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGt                                                                                       <  2:139907/150‑1 (MQ=255)
                                                               cAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTGGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTc                                                                                      >  2:1312443/1‑150 (MQ=255)
                                                                aGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCCCCAGCGCCCGCTGGTGATCa                                                                                                                                  <  1:848112/105‑1 (MQ=255)
                                                                aGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCCCCAGCGCCCGCTGGTGATCa                                                                                                                                  >  2:848112/1‑105 (MQ=255)
                                                                aGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCAGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGAGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCg                                                                                     >  1:733950/1‑150 (MQ=255)
                                                                 gCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCgg                                                                                    >  2:736346/1‑150 (MQ=255)
                                                                     aTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCcgc                                                                                >  2:1278279/1‑150 (MQ=255)
                                                                       cGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCgg                                                                              <  1:1728011/150‑1 (MQ=255)
                                                                       cGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCgg                                                                              <  1:1206277/150‑1 (MQ=255)
                                                                        gCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCggg                                                                             <  1:107747/150‑1 (MQ=255)
                                                                        gCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCggg                                                                             <  1:107711/150‑1 (MQ=255)
                                                                               tGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCAGCAGCGGGCGCTGGGGATCAGGTTGGCCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGa                                                                      >  2:744227/1‑150 (MQ=255)
                                                                                 ggCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGaaa                                                                    >  1:25890/1‑150 (MQ=255)
                                                                                    gAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATgg                                                                 <  1:1352068/150‑1 (MQ=255)
                                                                                     aaCAGGGAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATggg                                                                <  2:1526302/150‑1 (MQ=255)
                                                                                      aCAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATgggg                                                               <  2:711410/150‑1 (MQ=255)
                                                                                          gTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTg                                                                                              >  1:951224/1‑115 (MQ=255)
                                                                                          gTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGc                                                           >  2:701034/1‑150 (MQ=255)
                                                                                          gTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCAGGGAAACAAGTTGGGCGGGGCGCCTGCAGAGGTCTCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGc                                                           >  2:1159604/1‑150 (MQ=255)
                                                                                              gTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCt                                                                                               >  2:582927/1‑111 (MQ=255)
                                                                                              gTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCt                                                                                               <  1:582927/111‑1 (MQ=255)
                                                                                              gTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCggg                                                        >  2:1693001/1‑150 (MQ=255)
                                                                                                cGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGtt                                                      <  1:997960/150‑1 (MQ=255)
                                                                                                    aaaGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTc                                                                                      <  2:1605067/114‑1 (MQ=255)
                                                                                                    aaaGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTc                                                                                      >  1:1605067/1‑114 (MQ=255)
                                                                                                       gCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCGCCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAg                                               <  1:393756/150‑1 (MQ=255)
                                                                                                        cGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAgg                                              <  1:901031/150‑1 (MQ=255)
                                                                                                               aaTAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGCTGAGGCTACTcc                                       <  2:25890/150‑2 (MQ=255)
                                                                                                                     ggCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTgg                                                                                                                                        >  2:520289/1‑47 (MQ=255)
                                                                                                                     ggCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGa                                                >  2:387208/1‑135 (MQ=255)
                                                                                                                     ggCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGa                                                <  1:387208/135‑1 (MQ=255)
                                                                                                                     ggCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGgcgc                                 >  1:1546672/1‑150 (MQ=255)
                                                                                                                       cAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGAGCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCAt                               >  2:269002/1‑150 (MQ=255)
                                                                                                                           gTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTATATGAGGCTACTCATGGCGCATGTcc                           <  1:1483879/150‑1 (MQ=255)
                                                                                                                            tACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATAACGTATCGGGTTGGATGAGGCTACTCATGGCGCATGTCCt                          >  2:1894114/1‑150 (MQ=255)
                                                                                                                             aCCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGt                             <  1:1721052/146‑1 (MQ=255)
                                                                                                                             aCCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGt                             >  2:1721052/1‑146 (MQ=255)
                                                                                                                               cGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgtt                       >  1:650703/1‑150 (MQ=255)
                                                                                                                               cGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgtt                       >  2:1423913/1‑150 (MQ=255)
                                                                                                                                 aaGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgttgt                     >  2:915936/1‑150 (MQ=255)
                                                                                                                                      gCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCAt                               >  1:456755/1‑135 (MQ=255)
                                                                                                                                      gCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCAt                               <  2:456755/135‑1 (MQ=255)
                                                                                                                                      gCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTgc                >  1:1695439/1‑150 (MQ=255)
                                                                                                                                           gTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCaa           >  2:926675/1‑150 (MQ=255)
                                                                                                                                           gTTGATGATCACCAGCGGGAGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGGTGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCaa           >  2:885803/1‑150 (MQ=255)
                                                                                                                                               atgatCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACt       >  2:170200/1‑150 (MQ=255)
                                                                                                                                                   tCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCa   >  1:1752631/1‑150 (MQ=255)
                                                                                                                                                    cACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCaa  >  1:598400/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGAGGCCTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAA  >  NC_002947/5257589‑5257885

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