Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,436,760 G→T 11.7% T506T (ACG→ACT mnmC → fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,436,7600GT11.7% 818.7 / 70.0 247T506T (ACG→ACTmnmCfused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base G (135/83);  new base T (11/18);  total (146/101)
Fisher's exact test for biased strand distribution p-value = 1.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 248 total aligned reads displayed.

GTGGTACTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCAAAAGAGGTT  >  CP009273/2436611‑2436907
                                                                                                                                                     |                                                                                                                                                   
gTGGTACTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACg                                                                                                                                                     <  1:666180/150‑1 (MQ=255)
 tGGTACTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGc                                                                                                                                                    >  1:5049390/1‑150 (MQ=255)
 tGGTACTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGc                                                                                                                                                    >  1:2121347/1‑150 (MQ=255)
     aCTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCcaca                                                                                                                                                >  1:3069059/1‑150 (MQ=255)
      cTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGACTTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACaa                                                                                                                                               >  2:569915/1‑150 (MQ=255)
      cTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACaa                                                                                                                                               >  1:3624861/1‑150 (MQ=255)
      cTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACaa                                                                                                                                               >  1:3181082/1‑150 (MQ=255)
      cTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACaa                                                                                                                                               <  1:268335/150‑1 (MQ=255)
       tGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACaaa                                                                                                                                              >  1:1083590/1‑150 (MQ=255)
        ngCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACaaaa                                                                                                                                             >  1:3376071/2‑150 (MQ=255)
          cGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATc                                                                                                                                           >  2:2536884/1‑150 (MQ=255)
           gAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTAGTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCTGGTGCTGTGCTATGACGGTTATCTCCCGCCCCACAATcc                                                                                                                                          >  2:4906342/1‑150 (MQ=255)
             aCGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGACTTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACACAATCCgg                                                                                                                                        >  2:3129046/1‑150 (MQ=255)
             aCGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCgg                                                                                                                                        >  2:1299711/1‑150 (MQ=255)
                ggCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGa                                                                                                                                     >  2:852781/1‑150 (MQ=255)
                 gCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGaa                                                                                                                                    <  1:4540446/150‑1 (MQ=255)
                  cATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAAt                                                                                                                                   >  1:2819519/1‑150 (MQ=255)
                   aTCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGACTc                                                                                                                                  >  2:4422153/1‑150 (MQ=255)
                   aTCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATc                                                                                                                                  >  1:75948/1‑150 (MQ=255)
                    tCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACTCCACAAAATCCGGCGAATCa                                                                                                                                 <  2:2250863/150‑1 (MQ=255)
                           aGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAAcatcat                                                                                                                          >  1:1488118/1‑150 (MQ=255)
                           aGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAAcatcat                                                                                                                          <  1:4883712/150‑1 (MQ=255)
                            gCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCAtt                                                                                                                         <  2:2972838/150‑1 (MQ=255)
                              cGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGt                                                                                                                       <  1:3775185/150‑1 (MQ=255)
                               gATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTa                                                                                                                      >  1:4729234/1‑150 (MQ=255)
                               gATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTa                                                                                                                      >  1:2144447/1‑150 (MQ=255)
                                    aGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGt                                                                                                                 <  2:3440637/150‑1 (MQ=255)
                                        aaaCGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCa                                                                                                             >  2:1263782/1‑150 (MQ=255)
                                              cGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTatc                                                                                                       <  1:4611176/150‑1 (MQ=255)
                                                 ctctCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTatcatc                                                                                                    <  2:1029234/150‑1 (MQ=255)
                                                   ctcCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATcgc                                                                                                  >  1:1152509/1‑150 (MQ=255)
                                                       cGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCa                                                                                              >  1:4499311/1‑150 (MQ=255)
                                                       cGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCAGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCa                                                                                              >  1:820059/1‑150 (MQ=255)
                                                         gtgtATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGc                                                                                            <  1:3980574/150‑1 (MQ=255)
                                                          tgtATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCg                                                                                           <  1:2070448/150‑1 (MQ=255)
                                                           gtATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGa                                                                                          >  2:3385294/1‑150 (MQ=255)
                                                            tATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGaa                                                                                         >  2:3766488/1‑150 (MQ=255)
                                                                cGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGACTCAACATCATTGTCTTGGTGCCAGTTATCATCGCGGCAGCGAAGATa                                                                                     >  1:2325989/1‑150 (MQ=255)
                                                                cGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATa                                                                                     <  1:4000689/150‑1 (MQ=255)
                                                                cGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATa                                                                                     >  2:1665529/1‑150 (MQ=255)
                                                                 gggtGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATAc                                                                                    <  1:1783919/147‑1 (MQ=255)
                                                                  gTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACg                                                                                   >  1:2536259/1‑150 (MQ=255)
                                                                    tGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGc                                                                                 >  2:3078563/1‑150 (MQ=255)
                                                                    tGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGc                                                                                 >  1:3205406/1‑150 (MQ=255)
                                                                                 aGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGAc                                                                    <  2:2121347/150‑1 (MQ=255)
                                                                                  gCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACg                                                                   >  2:377691/1‑150 (MQ=255)
                                                                                    cATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGAt                                                                 >  1:2801291/1‑150 (MQ=255)
                                                                                    cATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGAt                                                                 <  2:446394/150‑1 (MQ=255)
                                                                                     atatTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATc                                                                <  2:3513736/150‑1 (MQ=255)
                                                                                     atatTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATc                                                                >  2:2415703/1‑150 (MQ=255)
                                                                                     atatTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATc                                                                >  2:975404/1‑150 (MQ=255)
                                                                                      tatTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATca                                                               <  2:1412954/150‑1 (MQ=255)
                                                                                       atTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATcag                                                              >  2:5063582/1‑150 (MQ=255)
                                                                                           cGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATcagcagc                                                          <  1:1299711/150‑1 (MQ=255)
                                                                                                  gCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCg                                                   >  1:2677411/1‑150 (MQ=255)
                                                                                                   ccGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGc                                                  >  2:1138186/1‑150 (MQ=255)
                                                                                                    cGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGcc                                                 >  2:4274194/1‑150 (MQ=255)
                                                                                                    cGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGcc                                                 <  2:279633/150‑1 (MQ=255)
                                                                                                    cGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTCTTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGACTCGcc                                                 >  1:1011946/1‑150 (MQ=255)
                                                                                                        aTTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCg                                             >  2:3986610/1‑150 (MQ=255)
                                                                                                          tGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGt                                           <  2:1334058/150‑1 (MQ=255)
                                                                                                          tGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGt                                           <  2:2909838/150‑1 (MQ=255)
                                                                                                              agagCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGACTCGCCAGCGGttgat                                       >  2:4310562/1‑150 (MQ=255)
                                                                                                                 gCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGttgattga                                    <  2:1399498/150‑1 (MQ=255)
                                                                                                                   tGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGGTgagtgatt                                  >  1:3285972/1‑150 (MQ=255)
                                                                                                                    gAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGattgattg                                 >  2:4496862/1‑150 (MQ=255)
                                                                                                                    gAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGattgattg                                 >  1:4434040/1‑150 (MQ=255)
                                                                                                                    gAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGagtgattg                                 >  2:3486064/1‑150 (MQ=255)
                                                                                                                     aaGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGt                                <  2:318215/150‑1 (MQ=255)
                                                                                                                      aGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGtt                               >  1:4914991/1‑150 (MQ=255)
                                                                                                                      aGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGtt                               <  2:700470/150‑1 (MQ=255)
                                                                                                                        cAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTc                             <  2:820059/150‑1 (MQ=255)
                                                                                                                         aGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTcc                            >  2:256221/1‑150 (MQ=255)
                                                                                                                         aGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTcc                            >  1:4159385/1‑150 (MQ=255)
                                                                                                                         aGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTcc                            >  2:4547447/1‑150 (MQ=255)
                                                                                                                         aGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTcc                            >  2:1625508/1‑150 (MQ=255)
                                                                                                                          ggTGCTGTGCTATGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTccc                           <  1:4144387/150‑1 (MQ=255)
                                                                                                                           gtgctgtgctATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCg                          >  2:259770/1‑150 (MQ=255)
                                                                                                                            tgctgtgctATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGc                         <  2:4951887/150‑1 (MQ=255)
                                                                                                                               tgtgctATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAgg                      >  1:1385327/1‑150 (MQ=255)
                                                                                                                               tgtgctATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAgg                      >  1:4710842/1‑150 (MQ=255)
                                                                                                                                gtgctATGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGc                     <  1:3253634/150‑1 (MQ=255)
                                                                                                                                gtgctATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCCGCGGGGGAGTGATTGTTTCCCGCAGGc                     >  2:3285619/1‑150 (MQ=255)
                                                                                                                                   ctATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGcaca                  <  2:4080399/150‑1 (MQ=255)
                                                                                                                                    tATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAg                 >  2:4847315/1‑150 (MQ=255)
                                                                                                                                     aTGACGGTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGt                >  2:4722916/1‑150 (MQ=255)
                                                                                                                                     aTGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGAGGTTGATTGATTGTTTCCCGCAGGCACAGt                >  2:3556919/1‑150 (MQ=255)
                                                                                                                                      tGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTg               <  2:4684999/150‑1 (MQ=255)
                                                                                                                                        aCGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTggg             >  1:3651269/1‑150 (MQ=255)
                                                                                                                                         cGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGc            <  2:5128863/150‑1 (MQ=255)
                                                                                                                                          ggTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCTCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCa           >  2:4160316/1‑150 (MQ=255)
                                                                                                                                          ggTTATCTCACTCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCa           >  2:4095912/1‑150 (MQ=255)
                                                                                                                                          ggTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCa           >  1:1195028/1‑150 (MQ=255)
                                                                                                                                            ttATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCaaa         >  1:1680347/1‑150 (MQ=255)
                                                                                                                                             tATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCaaaa        >  2:3645828/1‑150 (MQ=255)
                                                                                                                                             tATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCaaaa        >  2:2765803/1‑150 (MQ=255)
                                                                                                                                              aTCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGGTGATTGATTGTTTCCCGCAGGCACAGTGGGCACAag       >  2:1503929/1‑150 (MQ=255)
                                                                                                                                                ctcACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCACAagag     >  1:476598/1‑150 (MQ=255)
                                                                                                                                                ctcACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCAAAagag     >  1:125468/1‑150 (MQ=255)
                                                                                                                                                   aCGCCACAACATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCAAAAGAGGtt  <  2:404521/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTGGTACTGGCGAACGGGCATCAAATCAGCCGATTCAGCCAAACGTCGACTCTCCCGGTGTATTCGGTTGCCGGGCAGGTCAGCCATATTCCGACAACGCCGGAATTGGCAGAGCTGAAGCAGGTGCTGTGCTATGACGGTTATCTCACGCCACAAAATCCGGCGAATCAACATCATTGTATTGGTGCCAGTTATCATCGCGGCAGCGAAGATACGGCGTACAGTGAGGACGATCAGCAGCAGAATCGCCAGCGGTTGATTGATTGTTTCCCGCAGGCACAGTGGGCAAAAGAGGTT  >  CP009273/2436611‑2436907

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