Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,559,196 T→G 10.9% Q249P (CAG→CCG)  glgP ← glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,559,1960TG10.9% 557.7 / 4.1 246Q249P (CAG→CCG) glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (78/108);  new base G (32/2);  total (136/110)
Fisher's exact test for biased strand distribution p-value = 6.62e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

TTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACC  >  CP009273/3559047‑3559344
                                                                                                                                                     |                                                                                                                                                    
ttAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCt                                                                                                                                                      <  1:397923/150‑1 (MQ=255)
  aaaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTgg                                                                                                                                                    <  2:2241202/150‑1 (MQ=255)
   aaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGt                                                                                                                                                   <  1:886405/150‑1 (MQ=255)
   aaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTACACTGCCGCGAAGTAGTCACCCAGGt                                                                                                                                                   >  2:1666605/1‑150 (MQ=255)
    aaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGtt                                                                                                                                                  <  1:113344/150‑1 (MQ=255)
    aaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACAGCCGCGAAGTAGTAACCCTGGtt                                                                                                                                                  >  1:1563406/1‑150 (MQ=255)
     aTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTATTCACCCGGGTTg                                                                                                                                                 >  1:3464434/1‑150 (MQ=255)
      tGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGa                                                                                                                                                >  2:3852392/1‑150 (MQ=255)
      tGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGGTGa                                                                                                                                                >  2:853347/1‑150 (MQ=255)
      tGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACTGCCGCGTAGTAGTGACCCTGGTTgt                                                                                                                                                >  2:2392471/1‑149 (MQ=255)
          cTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGTTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAAttt                                                                                                                                            <  2:4433727/150‑1 (MQ=255)
             aaTGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTAcc                                                                                                                                         <  1:3647314/150‑1 (MQ=255)
             aaTGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTAcc                                                                                                                                         >  2:922836/1‑150 (MQ=255)
               tGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGa                                                                                                                                       >  2:1094010/1‑150 (MQ=255)
                 gTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAgg                                                                                                                                     <  1:2847902/150‑1 (MQ=255)
                    gAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTa                                                                                                                                  >  2:3327243/1‑150 (MQ=255)
                    gAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTa                                                                                                                                  >  2:3327105/1‑150 (MQ=255)
                    gAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGGTa                                                                                                                                  >  1:2549054/1‑150 (MQ=255)
                      ggAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAAt                                                                                                                                <  1:192038/150‑1 (MQ=255)
                       gAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTTGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAAtt                                                                                                                               <  2:419063/150‑1 (MQ=255)
                           ccAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTACCGAGGTTAATTTCGc                                                                                                                           >  1:4242544/1‑150 (MQ=255)
                            cAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTATTACACAGCCCCGAACTAGTAACCCAGGTTTAATTTAACGAGGTTAATTTAGCt                                                                                                                          >  1:475042/1‑150 (MQ=255)
                            cAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTACAAAGCCCCGAAGTAGTCACCCAGGTTGTATTTAACGAGGTTTATTTAGCt                                                                                                                          >  1:134311/1‑150 (MQ=255)
                             aGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTCTTACACAGCCGCGAAGTAGTCACCCGGGGTGAATTTACCGAGGTTTATTTCGCTa                                                                                                                         >  2:3575523/1‑150 (MQ=255)
                               gAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACt                                                                                                                       <  1:1045176/150‑1 (MQ=255)
                               gAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTAGCTaat                                                                                                                       >  2:485843/1‑148 (MQ=255)
                                 aGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTgg                                                                                                                     <  1:2927645/150‑1 (MQ=255)
                                 aGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTgg                                                                                                                     >  2:2560712/1‑150 (MQ=255)
                                  gTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCATCCACAGCCGCGAAGTAGTCACCCAGGGTGAATTTACCGAGGTTTATTTCGCTACAGGc                                                                                                                    >  1:3885679/1‑150 (MQ=255)
                                   tATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTTACGAGGTTAATTTCGCTACTGGcg                                                                                                                   >  2:1595751/1‑149 (MQ=255)
                                    aTTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCtt                                                                                                                  <  2:4476440/150‑1 (MQ=255)
                                      tCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTgc                                                                                                                <  1:4149240/150‑1 (MQ=255)
                                          gACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCCCAGCCGCGCAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCAc                                                                                                            >  2:2597851/1‑150 (MQ=255)
                                          gACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTCTTCCACAGCCGCGAAGTAGTCACCCCGGGTGAATTTACCGAGGTTAAATTCGCTACTGGCGTGCGCAc                                                                                                            >  1:1683734/1‑150 (MQ=255)
                                           aCGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACt                                                                                                           <  1:3477195/150‑1 (MQ=255)
                                            cGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGACGTAGTCACCCTGGTTGAATTTACCGAGGTTACTGTCGCTACTGGCGTGCGCCCTc                                                                                                          >  2:2632698/1‑150 (MQ=255)
                                            cGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTc                                                                                                          <  2:2522914/150‑1 (MQ=255)
                                                   gcAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTATTCCACAGGCCCGAACTAGTAACCCAGGTTTAATTTTAACAGGTTAATTTAGCAACTGGCATGAGCACACCACaaaa                                                                                                   >  1:1105690/1‑149 (MQ=255)
                                                     aGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGc                                                                                                 >  2:2679578/1‑150 (MQ=255)
                                                     aGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGc                                                                                                 >  1:4606114/1‑150 (MQ=255)
                                                      gCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCa                                                                                                <  1:2098193/150‑1 (MQ=255)
                                                       cTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTTAATTTACCGAGGTTAATTTAGCTACTGGCTTGCGCACTCCACAAACGCAg                                                                                               >  2:1509056/1‑150 (MQ=255)
                                                       cTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGACGCGAAGTAGTAACCCAGGTTGAATTTAACCAGGTTAATTTAGCTAATGGCTTTCGCACTCCACAAACACAg                                                                                               >  1:3956345/1‑150 (MQ=255)
                                                        tCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTAACCCTGGTTGAATTTAACGAGGTTAATTTCGCTTCGGGCTTGCGCACTCCACAAACGCAGc                                                                                              >  2:173345/1‑150 (MQ=255)
                                                          aCGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTAGTACACAGCCCCGAAATATTGACCCAGGTTGTAATTAACCAGGGTTATTTCGCGAAAGGCTTGCGCACACCACAAAccccaccc                                                                                            >  1:2851703/1‑143 (MQ=255)
                                                           cGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCgtg                                                                                           <  1:1069127/150‑1 (MQ=255)
                                                           cGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCgtg                                                                                           <  1:1186627/150‑1 (MQ=255)
                                                            gCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCAGGTTGAATTTACCGAGGTTAAGTTCGCTACAGGCTTGCGCACTCCACCAACGCcgcgggg                                                                                          >  2:3787393/1‑147 (MQ=255)
                                                            gCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTAACGAGGGTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCgtgt                                                                                          >  1:1947299/1‑150 (MQ=255)
                                                            gCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTACACAGCCGCGAAGTAGTCACCCTGGTTGAATTTAACGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCgtgt                                                                                          >  2:4018502/1‑150 (MQ=255)
                                                             cccGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCCCGAAGTAGTCACCCAGGGTGAATTTACCGATGTTAATTTCGCTACTGGCTTGCGCACACCCCAAACGCAGCGGGtt                                                                                         >  1:683602/1‑150 (MQ=255)
                                                                 gAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGtggtgggc                                                                                     >  2:4671845/1‑148 (MQ=255)
                                                                  aGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTcg                                                                                    <  1:3828878/150‑1 (MQ=255)
                                                                  aGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACAGCCCCGAACTAGTCACCCCGGTTTAAATTAACGAGGTTAATTTAGCTAAAGGCTTGCGCACACCACAAAACCACCGGGTTGGTcg                                                                                    >  2:1991432/1‑150 (MQ=255)
                                                                   gTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAAAGCAGCGTGTTGGTcgc                                                                                   >  2:4092761/1‑150 (MQ=255)
                                                                         ggAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCGACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgt                                                                             <  2:1685416/150‑1 (MQ=255)
                                                                           aGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgtgt                                                                           <  2:3541464/150‑1 (MQ=255)
                                                                           aGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgtgt                                                                           <  1:3831552/150‑1 (MQ=255)
                                                                           aGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATAATCCACAGCCGCGAACTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTAGCTACTGGCTTGCGCACACCACAAACCCACCGTGTTGGTCGCGGCGgtgt                                                                           >  1:2993836/1‑150 (MQ=255)
                                                                            gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTCTTACACAGCCGCGAAGTTGTCCCCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACCAACGCAGCGGGGTGGTCGCGTCGGTGTc                                                                          >  2:216787/1‑150 (MQ=255)
                                                                            gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCCCCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAAAGCAGCGGGTTGGTCGCGTCGGTGTc                                                                          >  2:617131/1‑150 (MQ=255)
                                                                             tCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCg                                                                         <  1:3360694/150‑1 (MQ=255)
                                                                              cATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACAGCCGCGAAGTAGTCACCCGGGTTGAATTTAACGAGGTTAATTTCGCTAAAGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGt                                                                        >  2:2897531/1‑150 (MQ=255)
                                                                               aTCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTa                                                                       <  2:531655/150‑1 (MQ=255)
                                                                                   ggATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCa                                                                   <  1:3202958/150‑1 (MQ=255)
                                                                                        cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGAt                                                              <  1:3683911/150‑1 (MQ=255)
                                                                                        cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACAGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGAt                                                              >  2:3437731/1‑150 (MQ=255)
                                                                                         aGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGAtt                                                             <  1:1502357/150‑1 (MQ=255)
                                                                                          gTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTa                                                            <  2:1386646/150‑1 (MQ=255)
                                                                                          gTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTa                                                            <  1:4263741/150‑1 (MQ=255)
                                                                                          gTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTa                                                            <  1:2494347/150‑1 (MQ=255)
                                                                                           tACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTAt                                                           <  1:2392471/150‑1 (MQ=255)
                                                                                           tACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTAt                                                           <  1:4342065/150‑1 (MQ=255)
                                                                                                 agaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATc                                                     <  2:3615892/150‑1 (MQ=255)
                                                                                                  gaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCg                                                    <  1:3327148/150‑1 (MQ=255)
                                                                                                  gaTACGTTCTCGGAGTGGTTTTTATATTACACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTAATGGCTTGCGCACACCACAAACGCCGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCg                                                    >  2:4575355/1‑150 (MQ=255)
                                                                                                    tACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTa                                                  <  1:599179/150‑1 (MQ=255)
                                                                                                      cGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAg                                                <  1:987134/150‑1 (MQ=255)
                                                                                                        ttCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAAATTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAAAGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGAACGTAAGCg                                              >  2:1127759/1‑150 (MQ=255)
                                                                                                         tctcGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGa                                             <  2:1105740/150‑1 (MQ=255)
                                                                                                           tcGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCGGGTTTAATTTACCGAGGTTAATTTCGCTACAGGCTTGCGCACACCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCAGATCGTAAGCGACt                                           >  2:4746205/1‑150 (MQ=255)
                                                                                                                gtgtTTTTTATCTTCCACTGCCGCGATGTAGTCACCCTGGTTTAATTTACCGTGGTTATTTTCGCTAGTGGGTTGCGCACTCAACAAACGCAGCTTGTTGGTCGCGTGGGTTTGGTAACCAGGGATTATCTTATCTTAAGCGACTCCCAg                                      <  2:1334921/150‑1 (MQ=255)
                                                                                                                  ggTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGTTTTATCTGATCGTAAGCGACTCCCAGaa                                    <  2:3478437/150‑1 (MQ=255)
                                                                                                                    tttttATTTTCCAGTGCCGCTATGTTGTCACCGTGGTTGATTTTACCGTGGTTATTTTCGCTACGGGCTTGCGCACGCCACAAACGCTGGGTGTTGGTCGCGTGGGTGTCGTAACCTGGTATTATCTGATCGTAAGCGACTCCCAGAAtc                                  <  2:4456396/150‑1 (MQ=255)
                                                                                                                     ttttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtct                                 >  2:1017847/1‑150 (MQ=255)
                                                                                                                     ttttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtct                                 >  2:4172868/1‑150 (MQ=255)
                                                                                                                       ttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATctct                               >  2:2950440/1‑150 (MQ=255)
                                                                                                                             tCCACTGCCGCGAAGTTGTAACCCTGGTTGAATTTACCGTGGTTATTTTCGCTACTGGCTTGCGCACTCCACAAACGCCGCGTGTTGGTCGCGTCGGTGTCGTACCCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGtt                         <  2:1237595/150‑1 (MQ=255)
                                                                                                                              ccACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGttt                        >  1:1989240/1‑150 (MQ=255)
                                                                                                                                        gAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAgcg              >  1:225202/1‑150 (MQ=255)
                                                                                                                                        gAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAgcg              >  1:224388/1‑150 (MQ=255)
                                                                                                                                        gAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAgcg              >  2:4636484/1‑150 (MQ=255)
                                                                                                                                          agtagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGcgcg            >  1:3055888/1‑150 (MQ=255)
                                                                                                                                            ttgtCACCGTGGTTTATTTTACCGTGGTTATTTTCGCTACGGGCTTGCGCACCCCAAAACGGCTGGGTGTTGGTGGCGTCGGTGTCTAAACCTGGTATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtt          <  1:386453/148‑1 (MQ=255)
                                                                                                                                             agtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACGGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCCGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGttt         <  2:2354678/150‑1 (MQ=255)
                                                                                                                                              gtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtttt        >  2:3966887/1‑150 (MQ=255)
                                                                                                                                              gtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtttt        >  2:3966771/1‑150 (MQ=255)
                                                                                                                                                 aCCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGttttttt     >  2:2707153/1‑150 (MQ=255)
                                                                                                                                                 aCCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGttttttt     >  1:803167/1‑150 (MQ=255)
                                                                                                                                                    cTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTAcc  >  1:4164210/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACC  >  CP009273/3559047‑3559344

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑