Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,557,688 C→A 16.8% D752Y (GAT→TAT)  glgP ← glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,557,6880CA16.8% 763.5 / 92.0 262D752Y (GAT→TAT) glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (110/108);  new base A (29/15);  total (139/123)
Fisher's exact test for biased strand distribution p-value = 6.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 263 total aligned reads displayed.

GATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAG  >  CP009273/3557542‑3557835
                                                                                                                                                  |                                                                                                                                                   
gATAGTACGGTCAGAAGAGAAGTAGCCCCTATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCg                                                                                                                                                    <  2:2049661/148‑1 (MQ=255)
   aGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCg                                                                                                                                                 >  2:1143925/1‑148 (MQ=255)
     tACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGaa                                                                                                                                               >  1:880129/1‑148 (MQ=255)
     tACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGaa                                                                                                                                               <  1:2762526/148‑1 (MQ=255)
       cGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGt                                                                                                                                             <  1:689939/148‑1 (MQ=255)
       cGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGt                                                                                                                                             >  1:3387993/1‑148 (MQ=255)
         gTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTg                                                                                                                                           >  2:1407561/1‑148 (MQ=255)
           cAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGATGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAAACCGCCAGTACCTGGTAGTGATCGCCGAAGTTGAt                                                                                                                                         <  2:2789664/148‑1 (MQ=255)
           cAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGAt                                                                                                                                         <  1:1266131/148‑1 (MQ=255)
            agaagagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATc                                                                                                                                        <  1:2472048/148‑1 (MQ=255)
              aagagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAg                                                                                                                                      >  1:1019268/1‑148 (MQ=255)
              aagagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAg                                                                                                                                      <  1:712858/148‑1 (MQ=255)
              aagagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAg                                                                                                                                      <  2:1214256/148‑1 (MQ=255)
                gagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCg                                                                                                                                    >  1:1717622/1‑148 (MQ=255)
                gagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCg                                                                                                                                    >  2:2813509/1‑148 (MQ=255)
                 agaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGa                                                                                                                                   >  2:1090633/1‑148 (MQ=255)
                 agaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGa                                                                                                                                   >  2:3518856/1‑148 (MQ=255)
                   aagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAAt                                                                                                                                 <  2:2974503/148‑1 (MQ=255)
                     ggaGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCa                                                                                                                               <  2:2367046/146‑1 (MQ=255)
                         cccATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCa                                                                                                                           >  2:3037241/1‑148 (MQ=255)
                              atTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATcg                                                                                                                      >  2:3115859/1‑148 (MQ=255)
                              atTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATcg                                                                                                                      >  1:149135/1‑148 (MQ=255)
                                tGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATcgcg                                                                                                                    >  2:2841844/1‑148 (MQ=255)
                                tGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATcgcg                                                                                                                    <  2:240410/148‑1 (MQ=255)
                                 ggCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGa                                                                                                                   <  1:2036522/148‑1 (MQ=255)
                                   cAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATa                                                                                                                 <  1:3117512/148‑1 (MQ=255)
                                   cAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATa                                                                                                                 >  2:639078/1‑148 (MQ=255)
                                   cAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATa                                                                                                                 >  2:142911/1‑148 (MQ=255)
                                     aTGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGc                                                                                                               >  2:1555298/1‑148 (MQ=255)
                                      tGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCg                                                                                                              >  1:1698370/1‑148 (MQ=255)
                                       gTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGa                                                                                                             <  1:1198795/148‑1 (MQ=255)
                                        ttCAGCATCGCTTTTGCGGTCCACTCTTCATGAAGCTGGTAGTTTTAATCGACTTTATCCTTACAATCTACATAGCTGCGATAATCCGCCAGTACGTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGAc                                                                                                            <  1:3239308/148‑1 (MQ=255)
                                         tCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGAcc                                                                                                           <  2:1736921/148‑1 (MQ=255)
                                               tCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATc                                                                                                     <  2:3187787/148‑1 (MQ=255)
                                               tCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATc                                                                                                     >  1:2326385/1‑148 (MQ=255)
                                                cGCTTTGGCGGTCCACACTTCCTGAAGCTGGTAGATTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCt                                                                                                    <  1:3388497/148‑1 (MQ=255)
                                                 gCTTTTGCGGTCCACTGTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCtt                                                                                                   <  1:3664287/148‑1 (MQ=255)
                                                 gCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCtt                                                                                                   <  2:2561093/148‑1 (MQ=255)
                                                   ttttGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTcc                                                                                                 >  2:3597134/1‑148 (MQ=255)
                                                     ttGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCgg                                                                                               >  1:1375116/1‑148 (MQ=255)
                                                           tCCACTCTTCCTGATGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGaa                                                                                         <  1:2434760/148‑1 (MQ=255)
                                                           tCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGaa                                                                                         <  2:1071594/148‑1 (MQ=255)
                                                            ccACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCTATAGCGACCCGGATCTTCCGGACTGAAt                                                                                        <  1:1831551/148‑1 (MQ=255)
                                                             cACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTCTCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATa                                                                                       <  1:69879/148‑1 (MQ=255)
                                                               ctctTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATaca                                                                                     >  1:3419533/1‑148 (MQ=255)
                                                                   tCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGc                                                                                 >  2:972310/1‑148 (MQ=255)
                                                                   tCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGc                                                                                 >  2:3122071/1‑148 (MQ=255)
                                                                       gAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGcc                                                                             >  2:1486871/1‑148 (MQ=255)
                                                                       gAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGcc                                                                             >  1:3163941/1‑148 (MQ=255)
                                                                         aGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGa                                                                           >  1:532564/1‑148 (MQ=255)
                                                                         aGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGa                                                                           <  2:1411652/148‑1 (MQ=255)
                                                                         aGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGa                                                                           <  2:1592734/148‑1 (MQ=255)
                                                                          gCTGGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGAt                                                                          <  2:2307005/148‑1 (MQ=255)
                                                                          gCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGAt                                                                          <  2:703484/148‑1 (MQ=255)
                                                                          gCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGAt                                                                          <  2:3411798/148‑1 (MQ=255)
                                                                           cTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGAtt                                                                         >  2:2702346/1‑148 (MQ=255)
                                                                                 gagTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTc                                                                   <  1:1906754/148‑1 (MQ=255)
                                                                                  agTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCa                                                                  >  2:377171/1‑148 (MQ=255)
                                                                                   gTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAg                                                                 <  1:1243874/148‑1 (MQ=255)
                                                                                        tCGACTTTATCCTGACAATCGCCATAGCTGCGAGAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCt                                                            <  2:54250/148‑1 (MQ=255)
                                                                                        tCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCt                                                            <  2:52130/148‑1 (MQ=255)
                                                                                         cgtcTTTTTCCTGACTTTCGTCATAGGTGCTATAATCCGCCAGTACCTGTTAGTTATCGCCGATGTTTTTCAGCGAATAAACCAGATCGCGTTAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTg                                                           <  1:2951039/145‑1 (MQ=255)
                                                                                         cGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTg                                                           >  1:34755/1‑148 (MQ=255)
                                                                                         cGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTg                                                           <  2:1698370/148‑1 (MQ=255)
                                                                                         cGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTg                                                           >  2:3231619/1‑148 (MQ=255)
                                                                                         cGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTg                                                           >  1:2198822/1‑148 (MQ=255)
                                                                                          gACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGa                                                          <  2:2326385/148‑1 (MQ=255)
                                                                                            cTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATg                                                        >  2:3269384/1‑148 (MQ=255)
                                                                                             tttATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGc                                                       <  1:1879302/148‑1 (MQ=255)
                                                                                              ttATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCa                                                      >  1:2683133/1‑148 (MQ=255)
                                                                                                 tCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGct                                                   >  1:1205143/1‑148 (MQ=255)
                                                                                                 tCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGct                                                   <  1:3791220/148‑1 (MQ=255)
                                                                                                  ccTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctc                                                  <  1:1143880/148‑1 (MQ=255)
                                                                                                   cTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctcc                                                 <  1:3637367/148‑1 (MQ=255)
                                                                                                    tGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctcct                                                >  1:353653/1‑148 (MQ=255)
                                                                                                      aCAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCa                                              <  2:743088/148‑1 (MQ=255)
                                                                                                       cAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCAt                                             <  2:1817688/148‑1 (MQ=255)
                                                                                                         aTCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCt                                           <  1:2385043/148‑1 (MQ=255)
                                                                                                             aCATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTtct                                       >  2:2252406/1‑148 (MQ=255)
                                                                                                               aTAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCg                                     >  2:3027869/1‑148 (MQ=255)
                                                                                                                tAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCgt                                    <  1:174741/148‑1 (MQ=255)
                                                                                                                tAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCgt                                    >  2:1653412/1‑148 (MQ=255)
                                                                                                                tAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCgt                                    <  2:3680241/148‑1 (MQ=255)
                                                                                                                tAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTGTCTCgt                                    >  1:2021540/1‑148 (MQ=255)
                                                                                                                  gCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCgtag                                  <  2:3387993/148‑1 (MQ=255)
                                                                                                                  gCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCgtag                                  <  2:3388849/148‑1 (MQ=255)
                                                                                                                     gCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTAt                               >  1:1530713/1‑148 (MQ=255)
                                                                                                                     gCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTAt                               <  1:2732593/148‑1 (MQ=255)
                                                                                                                        aTAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCa                            >  2:2025529/1‑148 (MQ=255)
                                                                                                                         tAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAc                           >  2:3491338/1‑148 (MQ=255)
                                                                                                                         tAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAc                           >  2:2920184/1‑148 (MQ=255)
                                                                                                                           aTCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAcgc                         >  1:60001/1‑148 (MQ=255)
                                                                                                                                 cAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTg                   >  1:84504/1‑148 (MQ=255)
                                                                                                                                  aGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGt                  >  1:1572564/1‑148 (MQ=255)
                                                                                                                                      ccTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGcc              <  2:2389650/148‑1 (MQ=255)
                                                                                                                                         ggTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTg           >  1:3015766/1‑148 (MQ=255)
                                                                                                                                         ggTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTg           <  2:1526432/148‑1 (MQ=255)
                                                                                                                                            aGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGTATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTgacg        <  2:1928334/148‑1 (MQ=255)
                                                                                                                                              tGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTgacgac      >  1:2643330/1‑148 (MQ=255)
                                                                                                                                                 tCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCa   >  2:1891623/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAG  >  CP009273/3557542‑3557835

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