Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,916,830 C→A 100% D71Y (GAT→TAT)  rsmG ← 16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,916,8300CA99.6% 1021.4 / ‑4.2 298D71Y (GAT→TAT) rsmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (133/163);  total (135/163)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 298 total aligned reads displayed.

TGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTG  >  CP009273/3916684‑3916974
                                                                                                                                                  |                                                                                                                                                
tGGTCCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGa                                                                                                                                                <  2:4696425/149‑1 (MQ=255)
   ttCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATaa                                                                                                                                             >  1:4661307/1‑149 (MQ=255)
   ttCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATaa                                                                                                                                             >  1:2658434/1‑149 (MQ=255)
    tCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATaaa                                                                                                                                            <  2:226075/149‑1 (MQ=255)
     cAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAAc                                                                                                                                           <  2:4189695/149‑1 (MQ=255)
      aaTATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAAcc                                                                                                                                          <  2:1953108/149‑1 (MQ=255)
       atatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCg                                                                                                                                         <  2:53272/149‑1 (MQ=255)
        tatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGt                                                                                                                                        <  2:4729534/149‑1 (MQ=255)
        tatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGt                                                                                                                                        <  2:3016254/149‑1 (MQ=255)
         atTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGtt                                                                                                                                       >  2:3631520/1‑149 (MQ=255)
         atTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGtt                                                                                                                                       >  1:242904/1‑149 (MQ=255)
          ttCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTc                                                                                                                                      >  2:3959750/1‑149 (MQ=255)
           tctcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCa                                                                                                                                     >  2:860476/1‑149 (MQ=255)
            ctcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCAc                                                                                                                                    <  1:2945244/149‑1 (MQ=255)
            ctcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCAc                                                                                                                                    >  1:1678208/1‑149 (MQ=255)
             tcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCAcc                                                                                                                                   <  2:3102545/149‑1 (MQ=255)
             tcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCAcc                                                                                                                                   >  1:954627/1‑149 (MQ=255)
              ccAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCt                                                                                                                                  >  1:896843/1‑149 (MQ=255)
               cAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCtt                                                                                                                                 <  2:2658420/149‑1 (MQ=255)
                aGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTg                                                                                                                                >  2:3738709/1‑149 (MQ=255)
                aGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTg                                                                                                                                <  2:1748176/149‑1 (MQ=255)
                 gTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGc                                                                                                                               <  2:2156155/149‑1 (MQ=255)
                 gTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGc                                                                                                                               >  2:4072168/1‑149 (MQ=255)
                      aGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATa                                                                                                                          <  1:4749332/149‑1 (MQ=255)
                      aGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATa                                                                                                                          <  2:3496523/149‑1 (MQ=255)
                       gCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATAc                                                                                                                         <  1:3050155/149‑1 (MQ=255)
                       gCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATAc                                                                                                                         >  2:2649921/1‑149 (MQ=255)
                          cATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGt                                                                                                                      <  2:4115238/149‑1 (MQ=255)
                            tGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGc                                                                                                                    >  2:4134122/1‑149 (MQ=255)
                            tGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGc                                                                                                                    >  2:3048083/1‑149 (MQ=255)
                            tGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGc                                                                                                                    <  1:277232/149‑1 (MQ=255)
                                    cTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACaa                                                                                                            <  1:4277489/149‑1 (MQ=255)
                                       aCGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGc                                                                                                         <  2:3407806/149‑1 (MQ=255)
                                        cGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCt                                                                                                        <  2:2974514/149‑1 (MQ=255)
                                          aaGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCCCAATGCTAt                                                                                                      >  1:4104955/1‑149 (MQ=255)
                                          aaGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTAt                                                                                                      >  1:3252232/1‑149 (MQ=255)
                                           aGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc                                                                                                     <  1:3419216/149‑1 (MQ=255)
                                               aaCGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCgaga                                                                                                 >  1:1840931/1‑149 (MQ=255)
                                                   cacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAatat                                                                                             >  2:3651579/1‑149 (MQ=255)
                                                   cacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAatat                                                                                             >  1:2352998/1‑149 (MQ=255)
                                                    acgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATg                                                                                            >  2:2771209/1‑149 (MQ=255)
                                                     cgcgTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATgg                                                                                           <  2:4373395/149‑1 (MQ=255)
                                                      gcgTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGc                                                                                          <  1:2257316/149‑1 (MQ=255)
                                                         tttACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTa                                                                                       >  2:2419163/1‑149 (MQ=255)
                                                          ttACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAc                                                                                      <  1:3797274/149‑1 (MQ=255)
                                                           tACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAcc                                                                                     <  1:3185967/149‑1 (MQ=255)
                                                           tACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAcc                                                                                     <  2:4275610/149‑1 (MQ=255)
                                                           tACCAAGACTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAcc                                                                                     <  2:1372203/149‑1 (MQ=255)
                                                            aCCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCa                                                                                    <  1:219908/149‑1 (MQ=255)
                                                                aGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAt                                                                                >  1:4192338/1‑149 (MQ=255)
                                                                  gCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAtct                                                                              >  2:3610155/1‑149 (MQ=255)
                                                                  gCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAtct                                                                              >  1:4787376/1‑149 (MQ=255)
                                                                  gCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAtct                                                                              <  1:817898/149‑1 (MQ=255)
                                                                       ccAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTa                                                                         <  1:300931/149‑1 (MQ=255)
                                                                       ccAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACAAGCATCTCATTa                                                                         >  1:1122605/1‑149 (MQ=255)
                                                                          aCAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGa                                                                      <  1:3959750/149‑1 (MQ=255)
                                                                           cAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGAt                                                                     <  2:1868620/149‑1 (MQ=255)
                                                                           cAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGAt                                                                     <  2:2563873/149‑1 (MQ=255)
                                                                            agagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATc                                                                    >  2:3295926/1‑149 (MQ=255)
                                                                            agagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATc                                                                    >  2:342831/1‑149 (MQ=255)
                                                                             gagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcg                                                                   >  2:1445173/1‑149 (MQ=255)
                                                                             gagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcg                                                                   >  1:3820627/1‑149 (MQ=255)
                                                                             gagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcg                                                                   >  1:1280105/1‑149 (MQ=255)
                                                                              agTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcgc                                                                  <  1:1072195/149‑1 (MQ=255)
                                                                               gTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcgcg                                                                 <  2:1595857/149‑1 (MQ=255)
                                                                                      gggCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGaa                                                          >  1:739729/1‑149 (MQ=255)
                                                                                      gggCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGaa                                                          <  1:3777491/149‑1 (MQ=255)
                                                                                        gCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGt                                                        <  1:1342887/149‑1 (MQ=255)
                                                                                        gCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGt                                                        >  1:2130886/1‑149 (MQ=255)
                                                                                         cTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTc                                                       >  1:1501666/1‑149 (MQ=255)
                                                                                          ttCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCa                                                      >  1:63942/1‑149 (MQ=255)
                                                                                              ggACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGtt                                                  >  1:1968281/1‑149 (MQ=255)
                                                                                                aCGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGt                                                <  1:3695438/149‑1 (MQ=255)
                                                                                                 cGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTa                                               <  1:4562235/149‑1 (MQ=255)
                                                                                                 cGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTa                                               <  2:612490/149‑1 (MQ=255)
                                                                                                    aCGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGc                                            <  1:4226664/149‑1 (MQ=255)
                                                                                                     cGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCt                                           <  2:3837529/149‑1 (MQ=255)
                                                                                                     cGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCt                                           >  2:1020143/1‑149 (MQ=255)
                                                                                                      gATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCtt                                          <  2:3865651/149‑1 (MQ=255)
                                                                                                        tAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCTCGGACCGAAGTCAGGTTGTACGCTTTg                                        >  1:4322206/1‑149 (MQ=255)
                                                                                                        tAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTg                                        <  2:879536/149‑1 (MQ=255)
                                                                                                          gagagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGtt                                      >  2:1218772/1‑149 (MQ=255)
                                                                                                           agagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTc                                     <  2:122537/149‑1 (MQ=255)
                                                                                                             agTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCa                                   >  1:2998638/1‑149 (MQ=255)
                                                                                                                   aTGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATg                             <  2:1568597/149‑1 (MQ=255)
                                                                                                                        tGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgcagca                        >  1:4000960/1‑149 (MQ=255)
                                                                                                                          gCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCata                      <  1:596618/149‑1 (MQ=255)
                                                                                                                           cAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCatat                     <  1:1891230/149‑1 (MQ=255)
                                                                                                                            aGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATAtt                    >  1:2536363/1‑149 (MQ=255)
                                                                                                                            aGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATAtt                    <  1:1626394/149‑1 (MQ=255)
                                                                                                                             gTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTc                   <  1:3738709/149‑1 (MQ=255)
                                                                                                                              tCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCa                  <  2:2039681/149‑1 (MQ=255)
                                                                                                                               ccTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCAc                 >  1:865279/1‑149 (MQ=255)
                                                                                                                                 tGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGt               >  2:1104323/1‑149 (MQ=255)
                                                                                                                                  ggTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTa              >  2:4148541/1‑149 (MQ=255)
                                                                                                                                    tCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAgg            >  1:3116816/1‑149 (MQ=255)
                                                                                                                                    tCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAgg            >  1:3612388/1‑149 (MQ=255)
                                                                                                                                         tGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCCTTTATGCAGCATATTCACGTAGGCAATa       >  2:3181154/1‑149 (MQ=255)
                                                                                                                                            cGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGc    <  1:4618596/149‑1 (MQ=255)
                                                                                                                                            cGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGc    <  1:2568921/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTG  >  CP009273/3916684‑3916974

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: