Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,335,001 Δ1 bp 6.4% coding (604/1704 nt) PP_2932 ← amidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,335,0010A. 6.4% 212.4 / 9.5 78coding (604/1704 nt)PP_2932amidase family protein
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base . (2/3);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

CCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGG  >  NC_002947/3334855‑3335146
                                                                                                                                                   |                                                                                                                                                 
ccACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGGCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgg                                                                                                                                                <  1:1209487/150‑1 (MQ=255)
 cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                               >  2:47739/1‑150 (MQ=255)
 cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTCGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                               >  2:1664425/1‑150 (MQ=255)
 cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGACCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                               >  2:1718552/1‑150 (MQ=255)
   caGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGgt                                                                                                                                             <  2:304603/150‑1 (MQ=255)
      gTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTa                                                                                                                                          <  2:1022163/150‑1 (MQ=255)
        cgccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgattt                                                                                                                                        <  1:1018682/150‑5 (MQ=255)
         gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGGGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                       >  1:1091240/1‑150 (MQ=255)
         gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                       <  1:690538/150‑1 (MQ=255)
         gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                       <  2:764070/150‑1 (MQ=255)
         gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                       >  2:1667511/1‑150 (MQ=255)
            gcgcTTGTCAGAGCCATCGCCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGca                                                                                                                                    <  2:900646/150‑1 (MQ=255)
              gcTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGcaca                                                                                                                                  <  1:89025/150‑1 (MQ=255)
                 tGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGc                                                                                                                               <  2:591227/150‑1 (MQ=255)
                     ggCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGCGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAAGCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgtt                                                                                                                           <  2:353309/150‑1 (MQ=255)
                      gCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgttg                                                                                                                          <  2:1066924/150‑1 (MQ=255)
                       ggTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgttgt                                                                                                                         <  2:1388720/149‑1 (MQ=255)
                         tCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CCTGCCGGAGTGTAGGCGCACGTGTCGtttttg                                                                                                                       <  2:1553502/150‑1 (MQ=255)
                            tCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGAc                                                                                                                    >  2:1615467/1‑150 (MQ=255)
                                     accacGTCGAGGATTTCAAGGAGGGCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCa                                                                                                           >  2:928851/1‑150 (MQ=255)
                                      ccacGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgat                                                                                                                                          >  1:190907/1‑116 (MQ=255)
                                      ccacGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgat                                                                                                                                          <  2:190907/118‑3 (MQ=255)
                                      ccacGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCAc                                                                                                          >  1:1620127/1‑150 (MQ=255)
                                           tCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCg                                                                                                     >  1:1680613/1‑150 (MQ=255)
                                             gAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCt                                                                                                   >  2:743077/1‑150 (MQ=255)
                                             gAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCt                                                                                                   >  1:1352471/1‑150 (MQ=255)
                                              aGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTg                                                                                                  >  2:426086/1‑150 (MQ=255)
                                              aGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTg                                                                                                  >  2:1086039/1‑150 (MQ=255)
                                               ggATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCTCCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTgg                                                                                                 <  2:1306923/150‑1 (MQ=255)
                                                   ttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTGGGACGCCGGGCCACGGCCGCTGGACCa                                                                                             >  1:1067026/1‑150 (MQ=255)
                                                      aaGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                          >  1:1045503/1‑150 (MQ=255)
                                                           ggTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTttctt                                                                                     <  2:1219426/150‑1 (MQ=255)
                                                            gTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTttcttc                                                                                    >  2:1610128/1‑150 (MQ=255)
                                                              cGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCgg                                                                                  >  1:1092668/1‑150 (MQ=255)
                                                              cGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCgg                                                                                  <  2:1352471/150‑1 (MQ=255)
                                                                 ccATGGTGTGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggcca                                                                               <  1:676346/150‑1 (MQ=255)
                                                                 ccATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggcca                                                                               <  1:1609411/150‑1 (MQ=255)
                                                                 ccATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggcca                                                                               <  2:1092668/150‑1 (MQ=255)
                                                                         cGGGCATAGGGCAACACCACGTCCATGGTCGGCATTAGCGAACAGTTAACAAACTTCGAAATCACCCCG‑CGCGACGGGTTTTAGACGCACAGGCCGTTGTTGGACGCCGGGCGACGGCCGCTGAACCATGTGTCTTCGTCCAGGCCAAAg                                                                       >  2:267260/1‑150 (MQ=255)
                                                                           ggCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGgc                                                                     >  1:372236/1‑150 (MQ=255)
                                                                              aTAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCgg                                                                                                               <  2:172110/105‑1 (MQ=255)
                                                                              aTAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCgg                                                                                                               >  1:172110/1‑105 (MQ=255)
                                                                                aGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGa                                                                <  1:743077/150‑1 (MQ=255)
                                                                                    caccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                               >  2:1750311/1‑67 (MQ=255)
                                                                                    caccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                               <  1:1750311/67‑1 (MQ=255)
                                                                                    caccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCCCAGGCCGTTGTTGGACGCCGGGCCACGGCCGTCGGTCCAGGTTTCTTCGGCCATGCCAAATTCGCTGAAGCt                                                            <  2:1163676/150‑1 (MQ=255)
                                                                                    caccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCt                                                            <  2:1045503/150‑1 (MQ=255)
                                                                                       caccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggc                                                         >  1:375249/1‑150 (MQ=255)
                                                                                       caccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggc                                                         >  1:795998/1‑150 (MQ=255)
                                                                                             gTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggt                                                   <  2:1770651/150‑1 (MQ=255)
                                                                                                cATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAc                                                                                               >  2:1300550/1‑103 (MQ=255)
                                                                                                cATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAc                                                                                               <  1:1300550/103‑1 (MQ=255)
                                                                                                cATGGTCGGCGTTAGCGGCAAGTAGCCACGCACCGAAATCCCCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggc                                                <  2:1841728/150‑1 (MQ=255)
                                                                                                 aTGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggc                                                                                 >  1:1288279/1‑116 (MQ=255)
                                                                                                 aTGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggc                                                                                 <  2:1288279/116‑1 (MQ=255)
                                                                                                 aTGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggcg                                               <  2:491196/150‑1 (MQ=255)
                                                                                                       ggCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑AGCCACGGGGTGAAAGGGCATAGGCCGTTGTTGGCAGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCg                                         <  1:426086/150‑1 (MQ=255)
                                                                                                           ttAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCAc                                     <  2:795069/150‑1 (MQ=255)
                                                                                                                         ccACGCACCGAAATCACCCCG‑CGCGCCGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGTCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGccc                       >  1:1907475/1‑150 (MQ=255)
                                                                                                                           aCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCg                                                                                   >  2:575418/1‑88 (MQ=255)
                                                                                                                           aCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCg                                                                                   <  1:575418/88‑1 (MQ=255)
                                                                                                                              cACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGcc                                                                                                      >  1:194942/1‑66 (MQ=255)
                                                                                                                              cACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGcc                                                                                                      <  2:194942/66‑1 (MQ=255)
                                                                                                                               aCCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCATGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCg                 >  1:27210/1‑150 (MQ=255)
                                                                                                                               aCCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACAAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCg                 >  2:1403212/1‑150 (MQ=255)
                                                                                                                                ccGAAATCACCCCG‑CGCGACGGGGTGTAGNCGCACAAGCCGTTGNTGGANGCCGGGCCACGGCCGCTGGCCCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGa                >  1:290464/1‑150 (MQ=255)
                                                                                                                                 cGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGaa               <  1:453591/150‑1 (MQ=255)
                                                                                                                                  gcctacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCg                                   <  2:1693816/124‑1 (MQ=255)
                                                                                                                                  gcctacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCg                                   >  1:1693816/6‑129 (MQ=255)
                                                                                                                                        aCCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGcgcg        <  2:709544/150‑1 (MQ=255)
                                                                                                                                          cccG‑CGCGACGGGGTGTAGGTGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGACAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt      >  1:226140/1‑150 (MQ=255)
                                                                                                                                          cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGTACGAGCCCGAGGCGAAGGGCGCGGt      >  2:1648094/1‑150 (MQ=255)
                                                                                                                                          cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt      <  2:1620127/150‑1 (MQ=255)
                                                                                                                                          cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt      >  1:681661/1‑150 (MQ=255)
                                                                                                                                          cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt      >  1:845743/1‑150 (MQ=255)
                                                                                                                                          cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAACCCGAGGCGAAGGGCGCGGt      >  2:438638/1‑150 (MQ=255)
                                                                                                                                            cg‑cgcgACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCa    <  1:928851/150‑1 (MQ=255)
                                                                                                                                               cgcgACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAgg  <  1:1718552/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGG  >  NC_002947/3334855‑3335146

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