Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 231,378 C→A 100% R26S (CGC→AGC)  yafS → putative S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273231,3780CA96.2% 1011.9 / 22.5 333R26S (CGC→AGC) yafSputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base A (151/168);  total (157/176)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 333 total aligned reads displayed.

GTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCC  >  CP009273/231231‑231523
                                                                                                                                                   |                                                                                                                                                 
 tAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCa                                                                                                                                                >  1:2753589/1‑150 (MQ=255)
    ttGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAgg                                                                                                                                             >  1:1700583/1‑150 (MQ=255)
     tGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGgc                                                                                                                                            <  1:159216/150‑1 (MQ=255)
      gTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGgcg                                                                                                                                           >  1:4553950/1‑150 (MQ=255)
       tcatcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGgcgc                                                                                                                                          >  2:1030563/1‑150 (MQ=255)
       tcatcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGgcgc                                                                                                                                          >  2:4208977/1‑150 (MQ=255)
             aaGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTTCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAgc                                                                                                                                    <  2:2265219/150‑1 (MQ=255)
             aaGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAgc                                                                                                                                    >  2:3242767/1‑150 (MQ=255)
                 cGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCa                                                                                                                                <  1:179759/150‑1 (MQ=255)
                  gggAATACTGTTAAGATTCATAGATTACCTCTCAGTGTTAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCACTCAAACAGTCGTGGCTAATGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTATATAGCAAGGCGCTGGAGCGCCAg                                                                                                                               >  1:4455919/1‑150 (MQ=255)
                  gggAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAg                                                                                                                               >  2:257024/1‑150 (MQ=255)
                  gggAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAg                                                                                                                               >  2:3320958/1‑150 (MQ=255)
                       tACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCaa                                                                                                                          >  2:3170919/1‑150 (MQ=255)
                          tGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAAccc                                                                                                                       <  2:2894038/150‑1 (MQ=255)
                           gTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAAccca                                                                                                                      <  2:214971/150‑2 (MQ=255)
                           gTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTACCTCAAACTGTCTTGGCTCATTATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCg                                                                                                                      >  1:2819679/1‑150 (MQ=255)
                             tAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTg                                                                                                                    >  2:1677769/1‑150 (MQ=255)
                             tAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTg                                                                                                                    >  2:1572712/1‑150 (MQ=255)
                              aaGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTgg                                                                                                                   <  1:3056770/150‑1 (MQ=255)
                              aaGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTgg                                                                                                                   >  1:3112083/1‑150 (MQ=255)
                               aGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGt                                                                                                                  <  1:1520149/150‑1 (MQ=255)
                                 aTTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTc                                                                                                                <  2:29578/150‑1 (MQ=255)
                                  ttCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa                                                                                                               <  2:2817656/150‑1 (MQ=255)
                                  ttCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa                                                                                                               >  1:5015888/1‑150 (MQ=255)
                                  ttCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa                                                                                                               >  1:773064/1‑150 (MQ=255)
                                  ttCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa                                                                                                               <  1:308910/150‑1 (MQ=255)
                                     aTAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTa                                                                                                            <  1:3398905/150‑1 (MQ=255)
                                       aGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaa                                                                                                          >  1:3774712/1‑150 (MQ=255)
                                       aGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaa                                                                                                          >  1:4087076/1‑150 (MQ=255)
                                        gATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaaa                                                                                                         >  1:4771952/1‑150 (MQ=255)
                                        gATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaaa                                                                                                         >  1:2884348/1‑150 (MQ=255)
                                          ttACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATg                                                                                                       >  2:2085618/1‑150 (MQ=255)
                                                ctcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                 <  2:1506534/150‑1 (MQ=255)
                                                ctcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                 >  2:4122193/1‑150 (MQ=255)
                                                ctcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                 <  2:1780846/150‑1 (MQ=255)
                                                ctcAGTGTGAAACGGAAGGTTGTGATGAAACAGGCAATAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGAGATTTGCCCTGTGGAAAGCTATATAGCAAGGCGCTGGAGCTCCAGCTCAACCCGTGGTTCACTAAAATGGCTGGt                                                                                                 >  1:2954575/1‑150 (MQ=255)
                                                  cAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGttt                                                                                               >  2:4785954/1‑150 (MQ=255)
                                                    gtgtGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTc                                                                                             <  2:1203677/150‑1 (MQ=255)
                                                    gtgtGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTc                                                                                             >  1:1867126/1‑150 (MQ=255)
                                                     tgtgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCa                                                                                            >  2:2205585/1‑150 (MQ=255)
                                                       tgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATc                                                                                          >  2:3066148/1‑150 (MQ=255)
                                                        gAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCt                                                                                         >  1:2403088/1‑150 (MQ=255)
                                                        gAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCt                                                                                         >  1:200971/1‑150 (MQ=255)
                                                          aaCGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGc                                                                                       <  2:3724598/150‑1 (MQ=255)
                                                          aaCGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGc                                                                                       >  1:479640/1‑150 (MQ=255)
                                                           aCGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCt                                                                                      <  1:654828/150‑1 (MQ=255)
                                                           aCGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCt                                                                                      <  2:2374910/150‑1 (MQ=255)
                                                            gggACGGTTGTGCTGAAACCGGCAAGAGTCCCTCACACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTTGGGAAAGATTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTACAATGTATGGTTTTCATCTGCtt                                                                                     <  2:5103182/149‑1 (MQ=255)
                                                            cGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCtt                                                                                     >  2:728724/1‑150 (MQ=255)
                                                            cGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCtt                                                                                     <  1:4328370/150‑1 (MQ=255)
                                                             ggAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTa                                                                                    <  2:5042368/150‑1 (MQ=255)
                                                                 ggTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAt                                                                                >  2:588852/1‑150 (MQ=255)
                                                                  gTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                               <  1:4479566/150‑1 (MQ=255)
                                                                    tgtgATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTgg                                                                             <  1:2513608/150‑1 (MQ=255)
                                                                    tgtgATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTgg                                                                             >  2:2663384/1‑150 (MQ=255)
                                                                    tgtgATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTgg                                                                             <  2:471391/150‑1 (MQ=255)
                                                                    tgtgATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTgg                                                                             <  2:4378805/150‑1 (MQ=255)
                                                                       gatgaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCaa                                                                          >  1:3795815/1‑150 (MQ=255)
                                                                            aaCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGTGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTaa                                                                     <  1:3564865/150‑1 (MQ=255)
                                                                            aaCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTaa                                                                     <  2:519644/150‑1 (MQ=255)
                                                                            aaCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTaa                                                                     >  1:4016584/1‑150 (MQ=255)
                                                                             aCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAg                                                                    <  2:3541876/150‑1 (MQ=255)
                                                                               cGGCGAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgcg                                                                  <  1:70012/150‑1 (MQ=255)
                                                                               cGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgcg                                                                  <  1:2173177/150‑1 (MQ=255)
                                                                                  cAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTAAGGTTTTCAACTGCTTAAGATTGGCAAATTAAGCGCAg                                                               >  1:4172976/1‑150 (MQ=255)
                                                                                   aaGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTCTGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGa                                                              >  2:1396503/1‑150 (MQ=255)
                                                                                   aaGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGa                                                              >  1:2941436/1‑150 (MQ=255)
                                                                                     gagTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGaaa                                                            >  1:2337231/1‑150 (MQ=255)
                                                                                     gagTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGaaa                                                            <  1:3147861/150‑1 (MQ=255)
                                                                                              aaaCTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCg                                                   >  1:5102879/1‑150 (MQ=255)
                                                                                               aaCTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGa                                                  <  2:4829705/150‑1 (MQ=255)
                                                                                                   gTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCg                                              <  1:3462321/150‑1 (MQ=255)
                                                                                                   gTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCg                                              >  1:938547/1‑150 (MQ=255)
                                                                                                     cGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTg                                            <  2:3361134/150‑1 (MQ=255)
                                                                                                      gTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgc                                           >  1:368848/1‑150 (MQ=255)
                                                                                                      gTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgc                                           <  2:2807483/150‑1 (MQ=255)
                                                                                                              tGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtc                                   <  2:4450634/150‑1 (MQ=255)
                                                                                                              tGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtc                                   >  1:673831/1‑150 (MQ=255)
                                                                                                                aTTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtctc                                 <  2:4480470/150‑1 (MQ=255)
                                                                                                                 ttGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtca                                >  1:1403323/1‑150 (MQ=255)
                                                                                                                    cTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcatca                             >  2:89160/1‑150 (MQ=255)
                                                                                                                    cTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcatca                             >  2:4273629/1‑150 (MQ=255)
                                                                                                                      ggggCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAg                           >  1:2908259/1‑150 (MQ=255)
                                                                                                                       gggCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATAGCGAAGCGTGCGCGGTTTCTCATCAAGt                          >  2:1840016/1‑150 (MQ=255)
                                                                                                                         gCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGa                        <  1:1030563/150‑1 (MQ=255)
                                                                                                                           gATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAAt                      <  1:2228267/150‑1 (MQ=255)
                                                                                                                           gATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAAt                      >  2:5068353/1‑150 (MQ=255)
                                                                                                                           gATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAAt                      <  1:5000714/150‑1 (MQ=255)
                                                                                                                           gATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAAt                      >  1:741107/1‑150 (MQ=255)
                                                                                                                            aTTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATg                     >  1:405107/1‑150 (MQ=255)
                                                                                                                               tGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGttt                  >  2:3448888/1‑150 (MQ=255)
                                                                                                                                gCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTc                 <  2:638521/150‑1 (MQ=255)
                                                                                                                                   cTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTgc              <  1:4578146/150‑1 (MQ=255)
                                                                                                                                     ggggAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTgcgc            >  2:3854970/1‑150 (MQ=255)
                                                                                                                                      gggAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCa           <  1:3250552/150‑1 (MQ=255)
                                                                                                                                      gggAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCa           >  1:317053/1‑150 (MQ=255)
                                                                                                                                          aaGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGa       <  2:3885208/150‑1 (MQ=255)
                                                                                                                                             cTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCTGTTTCTCATCAAGTGAATGTTTCTGCGCATGGAATg    <  2:1935973/150‑1 (MQ=255)
                                                                                                                                             cTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATg    >  1:4212786/1‑150 (MQ=255)
                                                                                                                                              tttATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGc   >  2:4624170/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCC  >  CP009273/231231‑231523

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: