Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,003,597 G→T 100% D100Y (GAC→TAC)  rlmL → fused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,003,5970GT98.4% 598.5 / 1.1 185D100Y (GAC→TAC) rlmLfused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (93/89);  total (93/92)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 185 total aligned reads displayed.

CCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCG  >  CP009273/1003452‑1003741
                                                                                                                                                 |                                                                                                                                                
 cgccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTaca                                                                                                                                               <  1:3336723/148‑1 (MQ=255)
    cTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCa                                                                                                                                            >  2:1434115/1‑148 (MQ=255)
      ggCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATc                                                                                                                                          <  1:2978315/148‑1 (MQ=255)
      ggCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATc                                                                                                                                          >  2:2515545/1‑148 (MQ=255)
          tCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCa                                                                                                                                      >  1:3541868/1‑148 (MQ=255)
           cgcgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCaa                                                                                                                                     >  2:277677/1‑148 (MQ=255)
                attatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtc                                                                                                                                <  2:276994/148‑1 (MQ=255)
                attatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtc                                                                                                                                <  1:2132825/148‑1 (MQ=255)
                  tatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCATGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTAAACCATCCGCAAcagtcag                                                                                                                              >  2:2108011/1‑148 (MQ=255)
                    tGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTa                                                                                                                            >  2:1661529/1‑148 (MQ=255)
                       tGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACgg                                                                                                                         <  1:377865/148‑1 (MQ=255)
                       tGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACgg                                                                                                                         >  2:301905/1‑148 (MQ=255)
                       tGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACgg                                                                                                                         <  1:2302147/148‑1 (MQ=255)
                          cGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGgcgc                                                                                                                      >  1:723239/1‑148 (MQ=255)
                           gCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGcgcg                                                                                                                     <  1:401654/148‑1 (MQ=255)
                           gCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGcgcg                                                                                                                     >  1:2041401/1‑148 (MQ=255)
                           gCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATACGCAACAGTCAGTACGGcgcg                                                                                                                     >  1:2590234/1‑148 (MQ=255)
                               ggCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtga                                                                                                                 >  2:1067732/1‑148 (MQ=255)
                               ggCGAGTGTAAGGTATACAGCGATTTAGACCTCAATCTCGGTTTTCATGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAAAGTCAGTACGGCGCGAtga                                                                                                                 >  1:3119957/1‑148 (MQ=255)
                                  gAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaag                                                                                                              <  2:3064988/148‑1 (MQ=255)
                                   aGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagt                                                                                                             >  2:2624290/1‑148 (MQ=255)
                                      gtAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaa                                                                                                          >  2:57007/1‑148 (MQ=255)
                                      gtAAGGTTTACAGCGAGTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaa                                                                                                          >  2:2894893/1‑148 (MQ=255)
                                        aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaag                                                                                                        <  1:2533518/148‑1 (MQ=255)
                                        aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaag                                                                                                        <  1:3660459/148‑1 (MQ=255)
                                        aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaag                                                                                                        >  1:3407108/1‑148 (MQ=255)
                                        aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaag                                                                                                        <  2:1529235/148‑1 (MQ=255)
                                        aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaag                                                                                                        >  1:1507554/1‑148 (MQ=255)
                                         aGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGa                                                                                                       <  2:3688588/148‑1 (MQ=255)
                                         aGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGa                                                                                                       <  1:1922105/148‑1 (MQ=255)
                                           gTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCACCAGTCAGTACGGCGCGATGAAAGAGAAAGAcg                                                                                                     >  2:2672910/1‑148 (MQ=255)
                                            tttACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAcgc                                                                                                    >  1:3362812/1‑148 (MQ=255)
                                              tACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGa                                                                                                  <  2:948841/148‑1 (MQ=255)
                                                 aGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtcg                                                                                               <  1:788234/148‑1 (MQ=255)
                                                 aGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtcg                                                                                               <  1:2103035/148‑1 (MQ=255)
                                                  gCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtcgt                                                                                              >  1:2469565/1‑148 (MQ=255)
                                                     aTTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCTCTGTCCACTTCAGTGGTTTGAATTACACCTTCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGa                                                                                           <  2:3331774/148‑1 (MQ=255)
                                                     aTTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGa                                                                                           <  2:3330425/148‑1 (MQ=255)
                                                       ttAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATg                                                                                         <  2:352496/148‑1 (MQ=255)
                                                       ttAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATg                                                                                         <  2:1619979/148‑1 (MQ=255)
                                                          gACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCtt                                                                                      <  1:1127074/148‑1 (MQ=255)
                                                          gACCTATATCTCGGTGTTCAGGCGATAAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTTGTTTGAATTACACCATCCGCAACAGTCAGTTCGGCGCGATGAAAGTGAATGACGCGATCGTCGATGCtt                                                                                      <  2:2862992/148‑1 (MQ=255)
                                                           aCCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCttt                                                                                     >  1:1258424/1‑148 (MQ=255)
                                                           aCCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCttt                                                                                     <  1:2393312/148‑1 (MQ=255)
                                                             ctctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCa                                                                                   >  2:2935893/1‑148 (MQ=255)
                                                              tctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAc                                                                                  >  2:3777921/1‑148 (MQ=255)
                                                               ctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAcg                                                                                 >  2:3253157/1‑148 (MQ=255)
                                                               ctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAcg                                                                                 <  2:1249911/148‑1 (MQ=255)
                                                               ctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAcg                                                                                 >  2:3252832/1‑148 (MQ=255)
                                                                 aTCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAcgcg                                                                               >  1:3756369/1‑148 (MQ=255)
                                                                   ctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGa                                                                             >  2:2511809/1‑148 (MQ=255)
                                                                     cGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaa                                                                           >  2:684593/1‑148 (MQ=255)
                                                                      ggTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaa                                                                          <  2:60336/148‑1 (MQ=255)
                                                                      ggTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaa                                                                          >  2:3278480/1‑148 (MQ=255)
                                                                           tCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGCCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTg                                                                     <  1:1023357/148‑1 (MQ=255)
                                                                              ggCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcg                                                                  >  2:1404292/1‑148 (MQ=255)
                                                                               gCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcgc                                                                 >  1:3657862/1‑148 (MQ=255)
                                                                                cGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcgcg                                                                <  1:3205934/148‑1 (MQ=255)
                                                                                  aTCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTc                                                              >  1:1405107/1‑148 (MQ=255)
                                                                                  aTCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTc                                                              <  2:53753/148‑1 (MQ=255)
                                                                                    cAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCa                                                            >  2:1875437/1‑148 (MQ=255)
                                                                                        tGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATg                                                        <  2:3541868/148‑1 (MQ=255)
                                                                                          gACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGtt                                                      >  1:2286114/1‑148 (MQ=255)
                                                                                           aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTg                                                     <  2:2277211/148‑1 (MQ=255)
                                                                                           aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTg                                                     >  1:491061/1‑148 (MQ=255)
                                                                                           aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTg                                                     <  1:3202858/148‑1 (MQ=255)
                                                                                            cAGAGATGTTTACTCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGa                                                    <  2:3269994/148‑1 (MQ=255)
                                                                                            cAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGa                                                    <  2:2552738/148‑1 (MQ=255)
                                                                                            cAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGa                                                    <  1:1707956/148‑1 (MQ=255)
                                                                                            cAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGa                                                    <  1:1707642/148‑1 (MQ=255)
                                                                                             agagATGTTTAATCCTGGCGCGACCGTCGCTGTCCACTTCAGTGGTTTGCATTACCCCATCCGCAACAGTCAGTAGGTCGAGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGAt                                                   <  2:1940364/148‑1 (MQ=255)
                                                                                                gaTGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGACAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATcgc                                                >  1:1960976/1‑148 (MQ=255)
                                                                                                  tGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGa                                              <  2:1709797/148‑1 (MQ=255)
                                                                                                      gaagccTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATgcg                                          <  2:196226/144‑1 (MQ=255)
                                                                                                         tCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCg                                       >  1:3063035/1‑148 (MQ=255)
                                                                                                          ccTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCgg                                      <  1:1024882/148‑1 (MQ=255)
                                                                                                          ccTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCgg                                      <  1:2931711/148‑1 (MQ=255)
                                                                                                           cTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGa                                     >  2:51485/1‑148 (MQ=255)
                                                                                                           cTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGa                                     >  2:51338/1‑148 (MQ=255)
                                                                                                            tGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGat                                    >  1:1631603/1‑148 (MQ=255)
                                                                                                              gcgcGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGatat                                  >  2:1368156/1‑148 (MQ=255)
                                                                                                              gcgcGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGatat                                  <  1:277677/148‑1 (MQ=255)
                                                                                                                   aCCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCcgcg                             <  2:824606/148‑1 (MQ=255)
                                                                                                                   aCCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCcgcg                             >  1:3310488/1‑148 (MQ=255)
                                                                                                                      ttCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTa                          <  2:3148075/148‑1 (MQ=255)
                                                                                                                      ttCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTa                          <  2:1201187/148‑1 (MQ=255)
                                                                                                                      ttCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTa                          <  1:2309026/148‑1 (MQ=255)
                                                                                                                       tCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTaa                         <  2:1800780/148‑1 (MQ=255)
                                                                                                                         gCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACg                       <  1:121504/148‑1 (MQ=255)
                                                                                                                         gCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACg                       <  1:2894893/148‑1 (MQ=255)
                                                                                                                         gCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACg                       <  1:301905/148‑1 (MQ=255)
                                                                                                                              ccACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTgg                  >  1:2346550/1‑148 (MQ=255)
                                                                                                                                 cTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTg               >  2:2930317/1‑148 (MQ=255)
                                                                                                                                       tGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaaa         >  1:1571567/1‑148 (MQ=255)
                                                                                                                                       tGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaaa         >  2:864427/1‑148 (MQ=255)
                                                                                                                                          tttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa      <  1:595140/148‑1 (MQ=255)
                                                                                                                                          tttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa      >  2:1383172/1‑148 (MQ=255)
                                                                                                                                           ttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaaa     >  1:3102932/1‑148 (MQ=255)
                                                                                                                                              aaTTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCg  >  2:3563694/1‑148 (MQ=255)
                                                                                                                                              aaTTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCg  >  1:2411277/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCG  >  CP009273/1003452‑1003741

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: