Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 456,298 T→A 100% L652Q (CTG→CAG)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273456,2980TA96.4% 944.8 / 17.7 278L652Q (CTG→CAG) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (8/2);  new base A (145/123);  total (153/125)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 278 total aligned reads displayed.

GACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATG  >  CP009273/456152‑456444
                                                                                                                                                  |                                                                                                                                                  
   ggAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGAt                                                                                                                                              >  1:127219/1‑150 (MQ=255)
   ggAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGAt                                                                                                                                              >  1:127266/1‑150 (MQ=255)
    gAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATc                                                                                                                                             >  2:4188576/1‑150 (MQ=255)
    gAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATc                                                                                                                                             >  2:2631442/1‑150 (MQ=255)
       gTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAAc                                                                                                                                          <  1:3157804/150‑1 (MQ=255)
        tGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAAcc                                                                                                                                         >  2:922064/1‑150 (MQ=255)
        tGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAAcc                                                                                                                                         <  1:172234/150‑1 (MQ=255)
          ggCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCt                                                                                                                                       >  2:666698/1‑150 (MQ=255)
          ggCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCt                                                                                                                                       >  2:4285873/1‑150 (MQ=255)
           gCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTg                                                                                                                                      <  1:1081298/150‑1 (MQ=255)
            ctgTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGa                                                                                                                                     <  2:1795890/148‑1 (MQ=255)
            cGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGa                                                                                                                                     <  1:1982087/150‑1 (MQ=255)
            cGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGa                                                                                                                                     <  1:209750/150‑1 (MQ=255)
            cGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGa                                                                                                                                     >  1:3880788/1‑150 (MQ=255)
             ggTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAACAACAGGGGATCAACCCTGAt                                                                                                                                    >  2:3628287/1‑150 (MQ=255)
             ggTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAt                                                                                                                                    >  2:480608/1‑150 (MQ=255)
             ggTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAt                                                                                                                                    >  1:3238674/1‑150 (MQ=255)
               tGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAttt                                                                                                                                  >  2:1979751/1‑150 (MQ=255)
               tGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAttt                                                                                                                                  <  2:802132/150‑1 (MQ=255)
               tGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAttt                                                                                                                                  >  2:3986357/1‑150 (MQ=255)
                gACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAtttt                                                                                                                                 <  2:671222/150‑1 (MQ=255)
                   ttGCTGACCATTGAAACAGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTAc                                                                                                                              <  2:2648687/150‑1 (MQ=255)
                       tGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGaaa                                                                                                                          >  2:1331273/1‑150 (MQ=255)
                          ccATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAAc                                                                                                                       >  2:2486723/1‑150 (MQ=255)
                          ccATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAAc                                                                                                                       <  1:3399750/150‑1 (MQ=255)
                           cATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACg                                                                                                                      <  1:1421643/150‑1 (MQ=255)
                           cATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACg                                                                                                                      >  1:3927946/1‑150 (MQ=255)
                            aTTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGt                                                                                                                     >  2:1090670/1‑150 (MQ=255)
                            aTTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGt                                                                                                                     >  1:2792097/1‑150 (MQ=255)
                             ttGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTg                                                                                                                    <  2:3994277/150‑1 (MQ=255)
                             ttGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTg                                                                                                                    <  2:1488097/150‑1 (MQ=255)
                             ttGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTg                                                                                                                    >  2:4463553/1‑150 (MQ=255)
                               gAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGAc                                                                                                                  <  1:3717637/150‑1 (MQ=255)
                               gAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGAc                                                                                                                  >  1:3354418/1‑150 (MQ=255)
                                aaaCCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACa                                                                                                                 >  2:3619598/1‑150 (MQ=255)
                                aaaCCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACa                                                                                                                 >  2:3619396/1‑150 (MQ=255)
                                    cGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAATCGTGACAtcca                                                                                                             >  2:3506783/1‑150 (MQ=255)
                                     gCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccac                                                                                                            >  1:3485141/1‑150 (MQ=255)
                                     gCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccac                                                                                                            >  2:728532/1‑150 (MQ=255)
                                        tgtgTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtc                                                                                                         >  2:1958368/1‑150 (MQ=255)
                                         gtgtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtcc                                                                                                        <  1:2764214/150‑1 (MQ=255)
                                         gtgtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtcc                                                                                                        <  1:2763628/150‑1 (MQ=255)
                                          tgtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtcca                                                                                                       >  1:1197177/1‑150 (MQ=255)
                                           gtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtccac                                                                                                      <  2:1446782/150‑1 (MQ=255)
                                            ttCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtccacg                                                                                                     <  2:1267091/150‑1 (MQ=255)
                                                       ggCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGt                                                                                          <  1:2819557/150‑1 (MQ=255)
                                                          aaaCTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCg                                                                                       >  2:2096787/1‑150 (MQ=255)
                                                          aaaCTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCg                                                                                       >  2:1783481/1‑150 (MQ=255)
                                                           aaCTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGa                                                                                      >  2:2189094/1‑150 (MQ=255)
                                                             cTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACg                                                                                    >  1:4244009/1‑150 (MQ=255)
                                                               gACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGAGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGcc                                                                                  <  1:3449648/150‑1 (MQ=255)
                                                                aCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCg                                                                                 <  2:3690721/150‑1 (MQ=255)
                                                                    ataCCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAg                                                                             >  1:829981/1‑150 (MQ=255)
                                                                        cGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATgg                                                                         <  1:2948587/150‑1 (MQ=255)
                                                                         ggTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGt                                                                        >  2:639585/1‑150 (MQ=255)
                                                                         ggTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGt                                                                        >  2:632065/1‑150 (MQ=255)
                                                                          gTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTc                                                                       <  2:2089203/150‑1 (MQ=255)
                                                                          gTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTc                                                                       >  1:354185/1‑150 (MQ=255)
                                                                          gTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTc                                                                       >  2:2559372/1‑150 (MQ=255)
                                                                           ttCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTcc                                                                      <  2:379306/150‑1 (MQ=255)
                                                                            tCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCg                                                                     >  1:552090/1‑150 (MQ=255)
                                                                            tCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCg                                                                     >  1:2339085/1‑150 (MQ=255)
                                                                             cGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGa                                                                    <  1:976419/150‑1 (MQ=255)
                                                                                tCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTg                                                                 <  2:3880665/150‑1 (MQ=255)
                                                                                 cGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGc                                                                <  1:96063/150‑1 (MQ=255)
                                                                                    cGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCgg                                                             >  2:1902434/1‑150 (MQ=255)
                                                                                     gAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGCAAGATGGTCCGAGTGCCGGt                                                            >  2:3052020/1‑150 (MQ=255)
                                                                                     gAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGt                                                            >  2:3051634/1‑150 (MQ=255)
                                                                                     gAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGt                                                            <  1:4355047/150‑1 (MQ=255)
                                                                                       aGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTAt                                                          >  1:4292945/1‑150 (MQ=255)
                                                                                       aGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTAt                                                          >  1:1368925/1‑150 (MQ=255)
                                                                                       aGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTAt                                                          <  2:2697328/150‑1 (MQ=255)
                                                                                       aGTGATGCAGGAGTCCATTCAGGCGGCGATAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTAt                                                          <  2:2697528/150‑1 (MQ=255)
                                                                                            tGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTa                                                     >  1:1503183/1‑150 (MQ=255)
                                                                                                ggAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAtgtg                                                 >  2:4446289/1‑150 (MQ=255)
                                                                                                   gTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAtctcaaa                                              >  1:1674876/1‑144 (MQ=255)
                                                                                                    tCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCAcc                                             >  2:3193641/1‑150 (MQ=255)
                                                                                                    tCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCAcc                                             <  2:605816/150‑1 (MQ=255)
                                                                                                                gcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCt                                 <  1:2180697/150‑1 (MQ=255)
                                                                                                                 cgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCtt                                <  2:3255451/150‑1 (MQ=255)
                                                                                                                 cgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCtt                                <  2:2297210/150‑1 (MQ=255)
                                                                                                                   ttAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGc                              >  1:1256887/1‑150 (MQ=255)
                                                                                                                    tAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGcc                             >  1:2646497/1‑150 (MQ=255)
                                                                                                                     aaCGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCt                            >  2:1692653/1‑150 (MQ=255)
                                                                                                                     aaCGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCt                            >  2:1691499/1‑150 (MQ=255)
                                                                                                                      aCGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTg                           >  2:3393403/1‑150 (MQ=255)
                                                                                                                        ggtggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGAc                         >  2:3401024/1‑150 (MQ=255)
                                                                                                                        ggtggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGAc                         >  2:1814351/1‑150 (MQ=255)
                                                                                                                           ggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCgg                      >  1:1019729/1‑150 (MQ=255)
                                                                                                                                gtgcgcgtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAAcc                 >  1:4119290/1‑150 (MQ=255)
                                                                                                                                   cgcgtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCgg              >  2:1373890/1‑150 (MQ=255)
                                                                                                                                     cgtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGtt            <  1:2024548/150‑1 (MQ=255)
                                                                                                                                      gtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTc           >  2:4105001/1‑150 (MQ=255)
                                                                                                                                       tgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCg          >  1:3597518/1‑150 (MQ=255)
                                                                                                                                         cgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTg        <  2:1539395/150‑1 (MQ=255)
                                                                                                                                          ggAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGc       >  2:2929999/1‑150 (MQ=255)
                                                                                                                                           gAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGcc      >  2:4536805/1‑150 (MQ=255)
                                                                                                                                           gAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGcc      <  2:2807034/150‑1 (MQ=255)
                                                                                                                                              aaaCAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGAt   >  2:3029892/1‑150 (MQ=255)
                                                                                                                                               aaCAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGAtg  >  2:3325908/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATG  >  CP009273/456152‑456444

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: