Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,248,553 A→G 5.3% S506S (AGT→AGC PP_3724 ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,248,5530AG 5.3% 294.0 / 3.8 94S506S (AGT→AGCPP_3724acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (42/47);  new base G (2/3);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

ACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGT  >  NC_002947/4248407‑4248701
                                                                                                                                                  |                                                                                                                                                    
aCCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTc                                                                                                                                                   >  2:1242212/1‑150 (MQ=255)
 ccAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCa                                                                                                                                                  >  2:847462/1‑150 (MQ=255)
 ccAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCa                                                                                                                                                  <  1:862564/150‑1 (MQ=255)
 ccAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCAACAGATTCGGTCTTGTTCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCa                                                                                                                                                  <  1:862057/150‑1 (MQ=255)
  cAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAg                                                                                                                                                 >  2:1130931/1‑150 (MQ=255)
    gTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCa                                                                                                                                               <  1:458444/150‑1 (MQ=255)
     tGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCAc                                                                                                                                              >  2:1511133/1‑150 (MQ=255)
       cAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCa                                                                                                                                               >  2:757535/1‑147 (MQ=255)
       cAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCa                                                                                                                                               <  1:757535/147‑1 (MQ=255)
       cAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGc                                                                                                                                            >  2:540432/1‑150 (MQ=255)
        aGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGcc                                                                                                                                           <  1:1242212/150‑1 (MQ=255)
         gCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGccc                                                                                                                                          <  1:1511133/150‑1 (MQ=255)
           ttCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTGCGGGCGTGCTGAGTACGCCCga                                                                                                                                        >  1:771921/1‑149 (MQ=255)
            tCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGCTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTATGCCACGGCGCGTGCGGGCGTGCTGAGTACGCCCgaa                                                                                                                                       <  1:1893453/150‑3 (MQ=255)
             ccAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTGCGGGCGTGCTGAGTACGCCCgaac                                                                                                                                      <  2:771921/150‑4 (MQ=255)
               aGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGc                                                                                                                                    >  2:537073/1‑150 (MQ=255)
               aGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGc                                                                                                                                    >  1:553676/1‑150 (MQ=255)
                gggCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCa                                                                                                                                   >  2:2120641/1‑150 (MQ=255)
                   cgcgGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCg                                                                                                                                >  1:1962911/1‑150 (MQ=255)
                    gcgGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCgg                                                                                                                               <  1:1130931/150‑1 (MQ=255)
                    gcgGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCgg                                                                                                                               <  1:640976/150‑1 (MQ=255)
                    gcgGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCgg                                                                                                                               <  1:1214469/150‑1 (MQ=255)
                          aTCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCCGGCACCGGCAGAcc                                                                                                                         >  1:536067/1‑150 (MQ=255)
                             gcgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGaa                                                                                                                      <  1:2045919/150‑1 (MQ=255)
                             gcgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACAGGCAGACCGaa                                                                                                                      >  1:959132/1‑150 (MQ=255)
                              cgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACGCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                     <  1:2004336/150‑1 (MQ=255)
                               gcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAgg                                                                                                                    >  2:1761705/1‑150 (MQ=255)
                                       aaCGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCa                                                                                                                                               >  1:950378/1‑115 (MQ=255)
                                       aaCGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCa                                                                                                                                               <  2:950378/115‑1 (MQ=255)
                                         cGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGc                                                                                                          <  2:1962911/150‑1 (MQ=255)
                                         cGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGc                                                                                                          <  2:182769/150‑1 (MQ=255)
                                         cGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGc                                                                                                          <  1:540432/150‑1 (MQ=255)
                                          gCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGc                                                                                                                              >  1:700847/1‑129 (MQ=255)
                                          gCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGc                                                                                                                              <  2:700847/129‑1 (MQ=255)
                                           ccACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCt                                                                                                        >  1:1180511/1‑150 (MQ=255)
                                                gCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTg                                                                                                   <  2:1180511/150‑1 (MQ=255)
                                                  aCCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGNACGCCCNCAGGCACCGGCANACCGAAGGCGCGCAGGGCCTGCTTGGa                                                                                                 <  1:2324086/150‑1 (MQ=255)
                                                  aCCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGa                                                                                                 <  1:2011528/150‑1 (MQ=255)
                                                    cAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCAc                                                                                                                                  >  1:1167165/1‑115 (MQ=255)
                                                    cAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCAc                                                                                                                                  <  2:1167165/115‑1 (MQ=255)
                                                    cAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGAtt                                                                                               <  2:553676/150‑1 (MQ=255)
                                                        ttCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAg                                                                                                             >  2:135858/1‑132 (MQ=255)
                                                        ttCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAg                                                                                                             <  1:135858/132‑1 (MQ=255)
                                                           ggTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGa                                                                                                 >  2:1691472/1‑141 (MQ=255)
                                                           ggTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGa                                                                                                 <  1:1691472/141‑1 (MQ=255)
                                                           ggTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTc                                                                                        <  1:2315789/150‑1 (MQ=255)
                                                            gTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCg                                                                                       >  2:41622/1‑150 (MQ=255)
                                                                 gtgCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCaa                                                                                  <  2:2163491/150‑1 (MQ=255)
                                                                  tgCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAAc                                                                                 <  1:2120641/150‑1 (MQ=255)
                                                                      gCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCtt                                                                                                    <  1:2139513/127‑1 (MQ=255)
                                                                      gCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCtt                                                                                                    >  2:2139513/1‑127 (MQ=255)
                                                                      gCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTc                                                                             <  1:847462/150‑1 (MQ=255)
                                                                       cAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTgg                                                                                                  >  2:1255597/1‑128 (MQ=255)
                                                                       cAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTgg                                                                                                  <  1:1255597/128‑1 (MQ=255)
                                                                          cTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGAc                                                                                                                          <  1:1749324/101‑1 (MQ=255)
                                                                          cTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGAc                                                                                                                          >  2:1749324/1‑101 (MQ=255)
                                                                          cTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGCCGTACTCAGCACGCCCGCAAGCACCCCCACCCCGCAGGCGGGCAGGGCCTACTTGGCTTGCCAGTCGTCCAACAGTCGGCt                                                                         >  2:255987/1‑150 (MQ=255)
                                                                           tGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTc                                                                                     <  2:358291/137‑1 (MQ=255)
                                                                           tGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTc                                                                                     >  1:358291/1‑137 (MQ=255)
                                                                                   gACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCCGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGt                                                                <  1:542554/150‑1 (MQ=255)
                                                                                      cGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAg                                                             >  2:115767/1‑150 (MQ=255)
                                                                                        ccTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATCGCCATTCGTCCAACCGTCGGCTTTCTCCCGTGAggg                                                           <  2:1785738/150‑1 (MQ=255)
                                                                                          ttCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGgcgcg                                                                                                                <  1:1703808/95‑1 (MQ=255)
                                                                                          ttCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGgcgcg                                                                                                                >  2:1703808/1‑95 (MQ=255)
                                                                                               tACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCgg                                                    <  1:380112/150‑1 (MQ=255)
                                                                                               tACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGCCACTGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCgg                                                    <  1:1449194/150‑1 (MQ=255)
                                                                                                  cAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCAGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGca                                                 >  1:78774/1‑150 (MQ=255)
                                                                                                     ggggTAACCGACGCGGGCGGGGTCGGCCACGGCGAGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGGGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAg                                              <  1:350771/150‑1 (MQ=255)
                                                                                                         tAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGCCGGCTTTCGCCCGTGAGGGCCt                                                        >  2:970077/1‑136 (MQ=255)
                                                                                                         tAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGCCGGCTTTCGCCCGTGAGGGCCt                                                        <  1:970077/136‑1 (MQ=255)
                                                                                                         tAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCAGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGt                                          >  1:127078/1‑150 (MQ=255)
                                                                                                            ccGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAcc                                       >  1:1025259/1‑150 (MQ=255)
                                                                                                            ccGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACAGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAcc                                       <  2:959132/150‑1 (MQ=255)
                                                                                                             cGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAcac                                      <  2:594386/150‑3 (MQ=255)
                                                                                                               aCGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGc                                    >  2:1658831/1‑150 (MQ=255)
                                                                                                                 gcgggcggCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGcac                                  >  1:1737255/1‑150 (MQ=255)
                                                                                                                    ggcggcGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCAGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGc                               <  2:78774/150‑1 (MQ=255)
                                                                                                                       ggcgTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTc                            >  1:820169/1‑150 (MQ=255)
                                                                                                                         cgTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGc                          >  2:1340692/1‑150 (MQ=255)
                                                                                                                             tGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTg                      >  2:1342382/1‑150 (MQ=255)
                                                                                                                              gCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGTCTGCTCTGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTgg                     <  1:41622/150‑1 (MQ=255)
                                                                                                                               ccACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGt                    <  1:1340692/150‑1 (MQ=255)
                                                                                                                               ccACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCCTTCGCCAGCAACGACCTCCATGCACAGCGGGACCCGCACAGCCTCGCCCTGGt                    >  2:22861/1‑150 (MQ=255)
                                                                                                                                 aCGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAccc                                      <  1:951589/130‑1 (MQ=255)
                                                                                                                                 aCGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAccc                                      >  2:951589/1‑130 (MQ=255)
                                                                                                                                 aCGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGACCTGGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTcc                  >  1:734490/1‑150 (MQ=255)
                                                                                                                                      gcgTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCaa             >  1:1595739/1‑150 (MQ=255)
                                                                                                                                       cgTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAAc            <  1:1761705/150‑1 (MQ=255)
                                                                                                                                         ttCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCg          >  2:1249539/1‑150 (MQ=255)
                                                                                                                                         tgcGGGCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCCCAGCCTCGCCCTGTTCCATCAAccc          >  2:851059/3‑149 (MQ=255)
                                                                                                                                           cGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCt        <  2:933508/150‑1 (MQ=255)
                                                                                                                                            gggCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTg       <  1:851059/150‑1 (MQ=255)
                                                                                                                                               cGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCg    <  2:1595739/150‑1 (MQ=255)
                                                                                                                                                 tACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGt  >  1:2057513/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGT  >  NC_002947/4248407‑4248701

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: