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,3780CA97.7% 738.7 / 5.9 235R26S (CGC→AGC) yafSputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (143/86);  total (146/89)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 237 total aligned reads displayed.

ATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGC  >  CP009273/231234‑231522
                                                                                                                                                |                                                                                                                                                
aTTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCa                                                                                                                                               <  1:2424054/148‑1 (MQ=255)
 ttGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCaa                                                                                                                                              >  2:1412558/1‑148 (MQ=255)
  tGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAg                                                                                                                                             >  2:3273835/1‑148 (MQ=255)
       tcaAAGGCGGGACTACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCt                                                                                                                                        >  2:3073714/1‑148 (MQ=255)
       tcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCt                                                                                                                                        >  1:3383575/1‑148 (MQ=255)
       tcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCt                                                                                                                                        >  1:2368833/1‑148 (MQ=255)
       tcaAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCt                                                                                                                                        >  1:1774113/1‑148 (MQ=255)
        caAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTg                                                                                                                                       >  1:2254320/1‑148 (MQ=255)
         aaaGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTgg                                                                                                                                      >  1:3225691/1‑148 (MQ=255)
           aGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAg                                                                                                                                    >  2:1093348/1‑148 (MQ=255)
            ggCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAgc                                                                                                                                   >  1:2284308/1‑148 (MQ=255)
               gggAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGcc                                                                                                                                <  1:2355154/148‑1 (MQ=255)
               gggAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGcc                                                                                                                                <  1:2355089/148‑1 (MQ=255)
                    tACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTc                                                                                                                           >  2:3176729/1‑148 (MQ=255)
                     aCTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCa                                                                                                                          >  2:1199614/1‑148 (MQ=255)
                     aCTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCa                                                                                                                          >  2:257823/1‑148 (MQ=255)
                        gTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAAcc                                                                                                                       >  2:3182386/1‑148 (MQ=255)
                          tAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCg                                                                                                                     <  2:406285/148‑1 (MQ=255)
                          tAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCg                                                                                                                     <  2:3134522/148‑1 (MQ=255)
                             gATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTgg                                                                                                                  >  2:1416864/1‑148 (MQ=255)
                             gATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTgg                                                                                                                  >  2:326737/1‑148 (MQ=255)
                               ttCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGtt                                                                                                                >  1:358562/1‑148 (MQ=255)
                                 cATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa                                                                                                              <  1:3073714/148‑1 (MQ=255)
                                 cATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa                                                                                                              <  2:2207826/148‑1 (MQ=255)
                                  aTAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCAc                                                                                                             <  2:163028/148‑1 (MQ=255)
                                    aGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTa                                                                                                           >  2:727048/1‑148 (MQ=255)
                                    aGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTa                                                                                                           >  1:332736/1‑148 (MQ=255)
                                    aGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCGCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTa                                                                                                           >  1:1035621/1‑148 (MQ=255)
                                     gATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaa                                                                                                          >  2:2680/1‑148 (MQ=255)
                                      aTTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaa                                                                                                         >  2:2620154/1‑148 (MQ=255)
                                      aTTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaa                                                                                                         >  2:1329959/1‑148 (MQ=255)
                                      aTTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaa                                                                                                         >  1:2837711/1‑148 (MQ=255)
                                       ttACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaaa                                                                                                        <  1:1586679/148‑1 (MQ=255)
                                             ctcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATg                                                                                                  >  1:84966/1‑148 (MQ=255)
                                              tcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATgg                                                                                                 <  1:2120749/148‑1 (MQ=255)
                                               cAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                <  1:938603/148‑1 (MQ=255)
                                               cAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                <  1:2002404/148‑1 (MQ=255)
                                               cAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGt                                                                                                >  2:2298670/1‑148 (MQ=255)
                                                    tgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCa                                                                                           <  2:1277700/148‑1 (MQ=255)
                                                    tgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCa                                                                                           <  1:2434467/148‑1 (MQ=255)
                                                    tgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCa                                                                                           <  2:2987746/148‑1 (MQ=255)
                                                    tgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCa                                                                                           >  1:1730319/1‑148 (MQ=255)
                                                     gAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCAt                                                                                          <  1:2517386/148‑1 (MQ=255)
                                                     gAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCAt                                                                                          <  2:1359200/148‑1 (MQ=255)
                                                           gAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCtt                                                                                    >  2:350094/1‑148 (MQ=255)
                                                            aaGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTa                                                                                   <  1:2602240/148‑1 (MQ=255)
                                                            aaGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTa                                                                                   >  1:1201375/1‑148 (MQ=255)
                                                               gTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGa                                                                                <  1:2856219/148‑1 (MQ=255)
                                                               gTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGa                                                                                <  2:623361/148‑1 (MQ=255)
                                                                 tgtgATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                              <  2:3253147/148‑1 (MQ=255)
                                                                   tgatgaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTgg                                                                            >  2:922927/1‑148 (MQ=255)
                                                                     atgaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCa                                                                          <  1:1466072/148‑1 (MQ=255)
                                                                      tgaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTTCTTAAGATTGGCaa                                                                         >  2:21276/1‑148 (MQ=255)
                                                                       gaAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAt                                                                        >  1:1307469/1‑148 (MQ=255)
                                                                          aCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTa                                                                     >  1:905812/1‑148 (MQ=255)
                                                                          aCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTa                                                                     <  1:2647231/148‑1 (MQ=255)
                                                                          aCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTa                                                                     >  1:234716/1‑148 (MQ=255)
                                                                            cGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAg                                                                   >  2:3149594/1‑148 (MQ=255)
                                                                             ggCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgc                                                                  <  2:2325682/148‑1 (MQ=255)
                                                                               cAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgcgc                                                                <  2:1237785/148‑1 (MQ=255)
                                                                               cAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgcgc                                                                >  1:1382241/1‑148 (MQ=255)
                                                                                 agagTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAg                                                              >  2:3192215/1‑148 (MQ=255)
                                                                                 agagTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAg                                                              >  2:1306800/1‑148 (MQ=255)
                                                                                  gagTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGa                                                             >  2:449115/1‑148 (MQ=255)
                                                                                   agTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTTGCAATTTAAGCGCAGaa                                                            >  1:511204/1‑148 (MQ=255)
                                                                                   agTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGaa                                                            <  1:1666785/148‑1 (MQ=255)
                                                                                           aaaCTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTg                                                    <  1:2179613/148‑1 (MQ=255)
                                                                                             aCTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCg                                                  >  1:2653662/1‑148 (MQ=255)
                                                                                               tGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGaa                                                >  2:2580782/1‑148 (MQ=255)
                                                                                               tGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCAGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAAAGTATTGTTTTCATATGCTTAAGATTTGCAATTTAAGCGCAGAAATCAATTTCGaa                                                >  1:1023714/1‑148 (MQ=255)
                                                                                                gTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAg                                               <  1:1873657/148‑1 (MQ=255)
                                                                                                 tCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGc                                              >  2:2074197/1‑148 (MQ=255)
                                                                                                      gCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgcg                                         >  2:2859179/1‑148 (MQ=255)
                                                                                                       cTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgcgc                                        >  2:2101552/1‑148 (MQ=255)
                                                                                                        tCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGcgcg                                       <  1:257823/148‑1 (MQ=255)
                                                                                                         ccTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                                      >  2:2925984/1‑148 (MQ=255)
                                                                                                         ccTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                                      >  2:798190/1‑148 (MQ=255)
                                                                                                         ccTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                                      >  2:1223018/1‑148 (MQ=255)
                                                                                                         ccTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                                      <  2:1495809/148‑1 (MQ=255)
                                                                                                           tGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGtt                                    >  1:1149415/1‑148 (MQ=255)
                                                                                                            gATTGCTGGGGCTATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCAAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGttt                                   <  1:3306970/148‑1 (MQ=255)
                                                                                                            gATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGttt                                   >  1:707929/1‑148 (MQ=255)
                                                                                                              ttGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtct                                 >  1:2972047/1‑148 (MQ=255)
                                                                                                              ttGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtct                                 >  1:3145647/1‑148 (MQ=255)
                                                                                                              ttGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtct                                 >  1:2139929/1‑148 (MQ=255)
                                                                                                              ttGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtct                                 <  1:2088741/148‑1 (MQ=255)
                                                                                                                 cTGGGGCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcat                              <  2:2715800/148‑1 (MQ=255)
                                                                                                                    gggCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCCTGGTTCACTAAAATGTATGGTTTTAATCTGCTTAAGATTGGCAATTTAAGCGCAGAAAACAATTTCGAAGCGTGCGCGGTTTCTCATCaa                           >  1:182751/1‑148 (MQ=255)
                                                                                                                     ggCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAg                          <  2:1812675/148‑1 (MQ=255)
                                                                                                                     ggCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAg                          <  2:3432980/148‑1 (MQ=255)
                                                                                                                      gCGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGt                         <  2:3404666/148‑1 (MQ=255)
                                                                                                                       cGATTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTg                        >  1:2055387/1‑148 (MQ=255)
                                                                                                                         aTTTGCCCTGGGGAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGaa                      >  1:568997/1‑148 (MQ=255)
                                                                                                                                    ggAAAGCTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTgcgc           >  2:649873/1‑148 (MQ=255)
                                                                                                                                          cTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGaa     >  1:572667/1‑148 (MQ=255)
                                                                                                                                          cTTTATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGaa     >  1:1001774/1‑148 (MQ=255)
                                                                                                                                           tttATAGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAAt    <  1:350094/148‑1 (MQ=255)
                                                                                                                                            ttATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATg   <  1:303219/148‑1 (MQ=255)
                                                                                                                                             tataGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGc  <  2:234716/148‑1 (MQ=255)
                                                                                                                                             tataGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGc  <  2:2534795/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGC  >  CP009273/231234‑231522

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: