Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,726,527 A→T 12.2% L560Q (CTG→CAG)  plsB ← glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,726,5270AT12.2% 278.4 / 25.1 98L560Q (CTG→CAG) plsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (42/44);  new base T (8/4);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 108 total aligned reads displayed.

CAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCAGGCGG  >  NC_002947/1726379‑1726676
                                                                                                                                                    |                                                                                                                                                     
caGCGCCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcag                                                                                                                                                      <  1:2264200/150‑1 (MQ=255)
  gCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcagca                                                                                                                                                    >  1:2226711/1‑150 (MQ=255)
  gCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcagca                                                                                                                                                    >  1:1794474/1‑150 (MQ=255)
    aCCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGAGAGTAGGGCACTTGGCGCtgctgc                                                                                                                                                  >  1:1356556/1‑150 (MQ=255)
      cgcgTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAgcgc                                                                                                                                                >  1:1950431/1‑150 (MQ=255)
         gTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAg                                                                                                                                             <  1:1176937/150‑1 (MQ=255)
         gTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAg                                                                                                                                             >  1:918537/1‑150 (MQ=255)
         gTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAg                                                                                                                                             <  2:212468/150‑1 (MQ=255)
         gTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAg                                                                                                                                             <  2:1950431/150‑1 (MQ=255)
         gTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGAGAGTAGGGCACTTGGCGCTGCTGCGCCaa                                                                                                                                             <  2:1356556/150‑2 (MQ=255)
            gtctgcATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCTGCTGCGCCAAGtg                                                                                                                                          <  2:727881/145‑2 (MQ=255)
             ccTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTAc                                                                                                                                         >  2:1877435/1‑150 (MQ=255)
              cTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACAGGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACa                                                                                                                                        >  1:1976645/1‑150 (MQ=255)
               ttCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAg                                                                                                                                       <  1:1046764/150‑1 (MQ=255)
                tCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAgg                                                                                                                                      <  2:619976/150‑1 (MQ=255)
                 cATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGt                                                                                                                                     <  2:2389476/150‑1 (MQ=255)
                    ccAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGa                                                                                                                                  <  2:338541/150‑1 (MQ=255)
                    ccAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGa                                                                                                                                  <  1:1877435/150‑1 (MQ=255)
                    ccAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGa                                                                                                                                  <  1:393039/150‑1 (MQ=255)
                    ccAGGTAGAGGATGCAGCCCTGGGCGTTCTTCTGCTCGGCCACCAGGTTCACGCTGCGCACATGTCCGTTAAGTGCCTGGCCGCCCCCCTGGGGCAGCGTGGTATGCGGGGTGTAGGGCACTTGGGGCAGCAGCGCCAGGTACAGGTCGa                                                                                                                                  <  2:738558/150‑1 (MQ=255)
                     cAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCGCCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAg                                                                                                                                 <  1:869568/150‑1 (MQ=255)
                          agagGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGTCAGGGTGGTGTGCGTCCTGTAGGGCGCGTTGCGTAGCAGCGCCAGGTACAGGTCGAGTAcgc                                                                                                                            <  2:276991/150‑1 (MQ=255)
                           gagGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCTGCTGCGCCaa                                                                                                                                             >  1:1598580/1‑131 (MQ=255)
                           gagGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCTGCTGCGCCaa                                                                                                                                             <  2:1598580/132‑2 (MQ=255)
                           gagGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTAcgcg                                                                                                                           <  1:732750/150‑1 (MQ=255)
                              gATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGt                                                                                                                        >  2:1088620/1‑150 (MQ=255)
                                 gCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAg                                                                                                                     <  2:1976645/150‑1 (MQ=255)
                                        aGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgcgcgc                                                                                                              <  2:239306/150‑1 (MQ=255)
                                          ggCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTc                                                                                                            <  2:1796810/150‑1 (MQ=255)
                                             gTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTc                                                                                                         <  1:1088620/150‑1 (MQ=255)
                                               ccTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCa                                                                                                       >  1:665028/1‑150 (MQ=255)
                                                cTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCgg                                                                                                                          <  2:481514/130‑1 (MQ=255)
                                                cTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCgg                                                                                                                          >  1:481514/1‑130 (MQ=255)
                                                 ttCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCCGTGCCTGCCCGTCGCCCTCGGGCAGGGCGGTGTGCGGCGCGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTGCAgg                                                                                                     >  1:1041744/1‑150 (MQ=255)
                                                 ttCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAgg                                                                                                     <  1:1522496/150‑1 (MQ=255)
                                                 ttCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAgg                                                                                                     <  2:2226711/150‑1 (MQ=255)
                                                       tCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAgg                                                                                               >  1:1732838/1‑150 (MQ=255)
                                                         ggCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAggcg                                                                                             >  2:502934/1‑150 (MQ=255)
                                                           ccaccaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAggcggc                                                                                           >  1:563743/1‑150 (MQ=255)
                                                               caGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtggt                                                                                       >  1:1673429/1‑150 (MQ=255)
                                                                  gTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGa                                                                                    >  1:828429/1‑150 (MQ=255)
                                                                   ttCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGAc                                                                                   >  1:640085/1‑150 (MQ=255)
                                                                       tgctgcGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGCGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGAcagca                                                                               >  2:1183257/1‑150 (MQ=255)
                                                                         ctgcGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGt                                                                             <  1:1183257/150‑1 (MQ=255)
                                                                           gcgcACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGc                                                                           <  2:665028/150‑1 (MQ=255)
                                                                             gcACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgc                                                                                                                  <  2:1027885/109‑1 (MQ=255)
                                                                             gcACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgc                                                                                                                  >  1:1027885/1‑109 (MQ=255)
                                                                               acaTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAgg                                                                       <  2:640085/150‑1 (MQ=255)
                                                                                caTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAgg                                                                                                     <  1:389360/119‑1 (MQ=255)
                                                                                caTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAgg                                                                                                     >  2:389360/1‑119 (MQ=255)
                                                                                 aTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGc                                                                     <  2:1794474/150‑1 (MQ=255)
                                                                                   gTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGcca                                                                   >  2:364105/1‑150 (MQ=255)
                                                                                       gATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAgg                                                                                                                    <  1:504987/97‑1 (MQ=255)
                                                                                       gATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAgg                                                                                                                    >  2:504987/1‑97 (MQ=255)
                                                                                        gtCAGGGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGCTCTGTGGGTAGGGCACTTGGCGGAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAgg                                                              <  2:2371002/149‑1 (MQ=255)
                                                                                        aTCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAgg                                                              >  1:639110/1‑150 (MQ=255)
                                                                                           aGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTg                                                           >  2:457256/1‑150 (MQ=255)
                                                                                            gTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGa                                                          >  2:1071478/1‑150 (MQ=255)
                                                                                              gCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcag                                                                                                                                                      <  1:124495/56‑1 (MQ=255)
                                                                                              gCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcag                                                                                                                                                      >  2:124495/1‑56 (MQ=255)
                                                                                              gCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGAcc                                                        >  2:1969135/1‑150 (MQ=255)
                                                                                                cTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCa                                                                                                       <  2:1020629/101‑1 (MQ=255)
                                                                                                cTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCa                                                                                                       >  1:1020629/1‑101 (MQ=255)
                                                                                                 tGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCggg                                                     <  2:828429/150‑1 (MQ=255)
                                                                                                  ggCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgc                                                                                                                  >  1:142851/1‑88 (MQ=255)
                                                                                                     cGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGACACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGa                                                 >  1:791750/1‑150 (MQ=255)
                                                                                                      gTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGAt                                                >  1:1312257/1‑150 (MQ=255)
                                                                                                         gCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggc                                             >  1:1037810/1‑150 (MQ=255)
                                                                                                           ccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCCGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcgg                                           >  1:2307456/1‑150 (MQ=255)
                                                                                                           ccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcgg                                           >  1:2192655/1‑150 (MQ=255)
                                                                                                           ccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcgg                                           >  1:704285/1‑150 (MQ=255)
                                                                                                           ccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATggcgg                                           >  1:2330658/1‑150 (MQ=255)
                                                                                                            cTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAg                                                               >  2:607622/1‑129 (MQ=255)
                                                                                                            cTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAg                                                               <  1:607622/129‑1 (MQ=255)
                                                                                                             tCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCggg                                                     >  2:482192/1‑138 (MQ=255)
                                                                                                               gggCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGgccgc                                       >  1:2326952/1‑150 (MQ=255)
                                                                                                                   aGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCg                                   >  1:66849/1‑150 (MQ=255)
                                                                                                                     ggtggtGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGtt                                 <  2:834595/150‑1 (MQ=255)
                                                                                                                       tggtgTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGa                               <  2:1037810/150‑1 (MQ=255)
                                                                                                                         gtgtGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAgg                             <  1:364105/150‑1 (MQ=255)
                                                                                                                           gtgCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCTGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTg                           >  1:1868535/1‑150 (MQ=255)
                                                                                                                           gtgCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTg                           <  2:2326952/150‑1 (MQ=255)
                                                                                                                                       ctgtacctGGCGCTGCTGCGCCAAGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGt               >  1:2274209/8‑150 (MQ=255)
                                                                                                                                        tgtacctGGCGCTGCTGCGCCAAGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAgg                             <  2:1292016/129‑1 (MQ=255)
                                                                                                                                        tgtacctGGCGCTGCTGCGCCAAGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAgg                             >  1:1292016/7‑135 (MQ=255)
                                                                                                                                         gCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGttt             <  2:66849/150‑1 (MQ=255)
                                                                                                                                          cACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGCTGATGGCGGCCGGCTCGTTGAGGTGGCGGGCCACGGTTtc            >  2:571033/1‑150 (MQ=255)
                                                                                                                                           aCTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACCCTGATTCCGGACCCAGCAGGTGGGCCGGGTTGATGGCGGCCGCCTCGTTGAGGGGGCGGGCCACGGTTtct           >  1:1852894/1‑150 (MQ=255)
                                                                                                                                            cctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCACCGCGGTGATGTGGCGGGCCACGGTTtctc          >  2:1508616/3‑150 (MQ=255)
                                                                                                                                            cTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTtctc          <  1:502934/150‑1 (MQ=255)
                                                                                                                                            cTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTtctc          <  2:704285/150‑1 (MQ=255)
                                                                                                                                             ttGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGa                               >  2:1910181/1‑128 (MQ=255)
                                                                                                                                             ttGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGa                               <  1:1910181/128‑1 (MQ=255)
                                                                                                                                             ctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCGGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTcc         >  1:147596/2‑150 (MQ=255)
                                                                                                                                             ctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTcc         >  1:2188564/2‑150 (MQ=255)
                                                                                                                                             ctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGCCCACGGTTTCTcc         >  1:101618/2‑150 (MQ=255)
                                                                                                                                             ctGGCGCTGCTGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTCCGGTCGCGCGCCACGGTTTCTcc         >  2:912114/2‑150 (MQ=255)
                                                                                                                                                  gcagcaGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCAggc    >  2:1918665/1‑150 (MQ=255)
                                                                                                                                                    agcagcGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGCCCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCAggcgg  <  1:807643/150‑1 (MQ=255)
                                                                                                                                                    agcagcGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCa       <  1:56772/145‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGACCGGGTTGATGGCGGCCGCCTCGTTGAGGTGGCGGGCCACGGTTTCTCCCAGGCGG  >  NC_002947/1726379‑1726676

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