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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,557,6880CA25.4% 530.5 / 137.5 232D752Y (GAT→TAT) glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (74/99);  new base A (25/34);  total (99/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 232 total aligned reads displayed.

GTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCA  >  CP009273/3557546‑3557834
                                                                                                                                              |                                                                                                                                                  
gTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGa                                                                                                                                               >  1:2194867/1‑148 (MQ=255)
 tACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGaa                                                                                                                                              >  2:2723097/1‑148 (MQ=255)
   cGGTCAGAAGAGAAGTAGCCCATCTTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACCTAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGt                                                                                                                                            <  1:1144762/148‑1 (MQ=255)
   cGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGt                                                                                                                                            >  2:2478306/1‑148 (MQ=255)
      tCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGa                                                                                                                                         <  2:2170157/148‑1 (MQ=255)
      tCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGa                                                                                                                                         >  1:1534669/1‑148 (MQ=255)
        agaagagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATc                                                                                                                                       >  1:1567322/1‑148 (MQ=255)
        agaagagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATc                                                                                                                                       <  2:1662219/148‑1 (MQ=255)
            gagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCg                                                                                                                                   >  2:1213750/1‑148 (MQ=255)
            gagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCg                                                                                                                                   >  1:2977025/1‑148 (MQ=255)
            gagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCg                                                                                                                                   >  2:2386137/1‑148 (MQ=255)
            gagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCg                                                                                                                                   <  2:2094200/148‑1 (MQ=255)
            gagaagTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGATGCGATAATCCGCCAGTACCAGGTATTGATCGCCGAAGGTGTTCAGCg                                                                                                                                   >  1:2230829/1‑148 (MQ=255)
                 gTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCa                                                                                                                              <  1:2010231/148‑1 (MQ=255)
                 gTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCa                                                                                                                              >  2:1151007/1‑148 (MQ=255)
                 gTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCa                                                                                                                              <  1:1588505/148‑1 (MQ=255)
                  tAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCaa                                                                                                                             <  1:905115/148‑1 (MQ=255)
                  tAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCaa                                                                                                                             <  1:905265/148‑1 (MQ=255)
                  tAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCaa                                                                                                                             >  2:869615/1‑148 (MQ=255)
                    gCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAAcc                                                                                                                           >  2:1297141/1‑148 (MQ=255)
                      ccATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAg                                                                                                                         <  1:2828816/148‑1 (MQ=255)
                       cATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGa                                                                                                                        <  1:305674/148‑1 (MQ=255)
                        atatTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGAt                                                                                                                       >  2:797603/1‑148 (MQ=255)
                        atatTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGAt                                                                                                                       <  1:2987526/148‑1 (MQ=255)
                           ttGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATcgc                                                                                                                    >  2:258616/1‑148 (MQ=255)
                            tGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATcgcg                                                                                                                   <  1:2120303/148‑1 (MQ=255)
                            tGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATcgcg                                                                                                                   >  1:1692317/1‑148 (MQ=255)
                             ggCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGa                                                                                                                  <  2:296929/148‑1 (MQ=255)
                              gCAATGTTCATCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGAt                                                                                                                 <  2:2080995/148‑1 (MQ=255)
                              gCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGAt                                                                                                                 >  2:140075/1‑148 (MQ=255)
                               cAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATa                                                                                                                <  1:258616/148‑1 (MQ=255)
                               cAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATa                                                                                                                <  1:2059094/148‑1 (MQ=255)
                                 aTGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGc                                                                                                              <  2:1817331/148‑1 (MQ=255)
                                   gTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGa                                                                                                            <  2:1163644/148‑1 (MQ=255)
                                    ttCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGAc                                                                                                           >  2:3522140/1‑148 (MQ=255)
                                    ttCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGAc                                                                                                           >  1:1151435/1‑148 (MQ=255)
                                     tCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGAcc                                                                                                          >  2:3560756/1‑148 (MQ=255)
                                      cAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGAccc                                                                                                         <  2:2947132/148‑1 (MQ=255)
                                      cAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGAccc                                                                                                         >  1:863022/1‑148 (MQ=255)
                                         cATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGa                                                                                                      >  2:273411/1‑148 (MQ=255)
                                         cATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGa                                                                                                      >  1:1187633/1‑148 (MQ=255)
                                         cATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGa                                                                                                      >  1:1188217/1‑148 (MQ=255)
                                          aTCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGAt                                                                                                     <  2:1348813/148‑1 (MQ=255)
                                           tCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATc                                                                                                    <  2:1245530/148‑1 (MQ=255)
                                                tttGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCg                                                                                               >  1:2307621/1‑148 (MQ=255)
                                                      gTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGa                                                                                         <  1:1149588/148‑1 (MQ=255)
                                                       tCCACTCTTCCTTATGCTGGTAGAGTTCATCTACTTTATCCTGAAAATCGACATAGCTGCTATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGaa                                                                                        <  1:2710130/148‑1 (MQ=255)
                                                       tCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCCGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGaa                                                                                        <  2:3304995/148‑1 (MQ=255)
                                                        ccACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAAt                                                                                       >  2:3146156/1‑148 (MQ=255)
                                                         cACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATa                                                                                      <  2:2771454/148‑1 (MQ=255)
                                                          aCTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATac                                                                                     >  1:1885859/1‑148 (MQ=255)
                                                          aCTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATac                                                                                     <  1:869615/148‑1 (MQ=255)
                                                           ctctTCCTTAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATCGCTGCGATAATCCGCCAGTACCTGGTAGGGATCGCCGAAGTTGATCAGCTAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATaca                                                                                    <  1:1024038/148‑1 (MQ=255)
                                                           ctctTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATaca                                                                                    <  2:880404/148‑1 (MQ=255)
                                                              ttCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCg                                                                                 <  2:3392474/148‑1 (MQ=255)
                                                               tCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGc                                                                                >  2:686030/1‑148 (MQ=255)
                                                                       cTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGAtt                                                                        <  1:2375261/148‑1 (MQ=255)
                                                                            agagTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGt                                                                   >  2:1455202/1‑148 (MQ=255)
                                                                            agagTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGt                                                                   >  1:3455148/1‑148 (MQ=255)
                                                                               gTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAg                                                                <  2:596336/148‑1 (MQ=255)
                                                                                ttCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGc                                                               <  2:1710694/148‑1 (MQ=255)
                                                                                   aTCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCAcc                                                            >  1:535832/1‑148 (MQ=255)
                                                                                         tttATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGc                                                      <  1:2741500/148‑1 (MQ=255)
                                                                                           tATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAg                                                    <  2:1579027/148‑1 (MQ=255)
                                                                                              ccTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctc                                                 <  2:1092153/148‑1 (MQ=255)
                                                                                              ccTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctc                                                 <  2:2649604/148‑1 (MQ=255)
                                                                                               cTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctcc                                                <  1:562358/148‑1 (MQ=255)
                                                                                                tGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctcct                                               <  2:750778/148‑1 (MQ=255)
                                                                                                 gACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctcctc                                              <  1:2077957/148‑1 (MQ=255)
                                                                                                 gACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctcctc                                              >  1:2440376/1‑148 (MQ=255)
                                                                                                 gACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGctcctc                                              <  2:460484/148‑1 (MQ=255)
                                                                                                  aCAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCa                                             >  2:1501934/1‑148 (MQ=255)
                                                                                                   cAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCAt                                            <  2:2389081/148‑1 (MQ=255)
                                                                                                    aaTCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATc                                           <  2:3533770/148‑1 (MQ=255)
                                                                                                      tCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCTAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCtt                                         <  1:2334596/148‑1 (MQ=255)
                                                                                                      tCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCtt                                         >  1:894570/1‑148 (MQ=255)
                                                                                                      tCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCtt                                         >  2:3076353/1‑148 (MQ=255)
                                                                                                          cATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCAGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTtctc                                     >  2:152687/1‑148 (MQ=255)
                                                                                                           aTAGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCg                                    <  1:1492544/148‑1 (MQ=255)
                                                                                                             aGCTGCGATAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCgta                                  <  2:41831/148‑1 (MQ=255)
                                                                                                              gCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCgtag                                 >  1:386820/1‑148 (MQ=255)
                                                                                                                tGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCgtagta                               >  2:1801599/1‑148 (MQ=255)
                                                                                                                   gATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCTAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTc                            <  1:328157/148‑1 (MQ=255)
                                                                                                                    aTAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCa                           <  1:3530279/148‑1 (MQ=255)
                                                                                                                    aTAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCa                           <  1:1269237/148‑1 (MQ=255)
                                                                                                                     tAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAc                          <  1:1746772/148‑1 (MQ=255)
                                                                                                                     tAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAc                          >  2:3271665/1‑148 (MQ=255)
                                                                                                                     tAATCCGCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAc                          >  1:2612622/1‑148 (MQ=255)
                                                                                                                         ccGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCgg                      <  2:1110365/148‑1 (MQ=255)
                                                                                                                           gCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAAACAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGtt                    <  1:1213599/148‑1 (MQ=255)
                                                                                                                           gCCAGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGtt                    <  1:394627/148‑1 (MQ=255)
                                                                                                                              aGTACCTGGTAGTGATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGt                 <  2:886591/148‑1 (MQ=255)
                                                                                                                               gTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATCGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTa                <  1:327060/148‑1 (MQ=255)
                                                                                                                               gTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTa                <  1:1985413/148‑1 (MQ=255)
                                                                                                                               gTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTa                >  1:594388/1‑148 (MQ=255)
                                                                                                                                 aCCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGc              >  1:690638/1‑148 (MQ=255)
                                                                                                                                   cTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGccc            >  1:3107131/1‑148 (MQ=255)
                                                                                                                                          tGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTgacgac     >  2:1846658/1‑148 (MQ=255)
                                                                                                                                           gATAGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGacgacg    <  2:1941316/148‑1 (MQ=255)
                                                                                                                                             tCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCa  >  2:3056743/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCATCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCA  >  CP009273/3557546‑3557834

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: