Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,024,827 G→T 91.5% D103Y (GAT→TAT)  pepQ → proline dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,024,8270GT91.5% 528.7 / 40.9 203D103Y (GAT→TAT) pepQproline dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (6/12);  new base T (94/91);  total (100/103)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

CCGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGGTTA  >  CP009273/4024680‑4024974
                                                                                                                                                   |                                                                                                                                                   
   cAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtg                                                                                                                                                  >  1:2092509/1‑148 (MQ=255)
   cAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtg                                                                                                                                                  >  1:1814968/1‑148 (MQ=255)
   cAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtg                                                                                                                                                  >  1:2861110/1‑148 (MQ=255)
      ttCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTgg                                                                                                                                               >  1:1278396/1‑148 (MQ=255)
      ttCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTgg                                                                                                                                               <  2:60399/148‑1 (MQ=255)
       tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                              <  2:2301272/148‑1 (MQ=255)
       tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACGGAATATGTGGa                                                                                                                                              >  2:2256841/1‑148 (MQ=255)
        cAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGaa                                                                                                                                             >  1:2691802/1‑148 (MQ=255)
          aaGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGt                                                                                                                                           >  2:494542/1‑148 (MQ=255)
          aaGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGt                                                                                                                                           <  2:3018919/148‑1 (MQ=255)
           aGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTg                                                                                                                                          <  1:772185/148‑1 (MQ=255)
            gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGa                                                                                                                                         >  1:1544157/1‑148 (MQ=255)
            gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGa                                                                                                                                         >  2:2482670/1‑148 (MQ=255)
             cGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGAt                                                                                                                                        <  1:3057584/148‑1 (MQ=255)
             cGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGAGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGAt                                                                                                                                        <  1:2495105/148‑1 (MQ=255)
                 ggTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATcgcg                                                                                                                                    >  2:1586998/1‑148 (MQ=255)
                 ggTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATcgcg                                                                                                                                    <  2:395058/148‑1 (MQ=255)
                  gTGCCGGGAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCgcgc                                                                                                                                   >  2:1728858/1‑148 (MQ=255)
                           aCTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAg                                                                                                                          <  1:2013786/148‑1 (MQ=255)
                             tCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGcc                                                                                                                        >  1:2985026/1‑148 (MQ=255)
                              cAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCg                                                                                                                       <  1:2568322/148‑1 (MQ=255)
                                   gCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGc                                                                                                                  <  1:1128777/148‑1 (MQ=255)
                                    ccAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCa                                                                                                                 <  1:99331/148‑1 (MQ=255)
                                    ccAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCa                                                                                                                 <  1:375413/148‑1 (MQ=255)
                                      aaaCTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCAtt                                                                                                               >  1:165357/1‑148 (MQ=255)
                                        aCTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTgg                                                                                                             >  1:1947966/1‑148 (MQ=255)
                                          tgctggttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTa                                                                                                           <  2:826276/148‑1 (MQ=255)
                                               gttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctg                                                                                                      <  2:2641508/148‑1 (MQ=255)
                                               gttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctg                                                                                                      <  2:3219572/148‑1 (MQ=255)
                                                   ctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgc                                                                                                  <  1:3272175/148‑1 (MQ=255)
                                                      gtggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctg                                                                                               <  2:2874578/148‑1 (MQ=255)
                                                       tggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctgc                                                                                              <  2:1271135/148‑1 (MQ=255)
                                                              cGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcgcgc                                                                                       >  1:1174992/1‑148 (MQ=255)
                                                                  aaCAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCa                                                                                   >  1:2993438/1‑148 (MQ=255)
                                                                  aaCAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCa                                                                                   <  2:2525833/148‑1 (MQ=255)
                                                                    cAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAat                                                                                 >  2:2371481/1‑148 (MQ=255)
                                                                      aGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAatat                                                                               <  1:442515/148‑1 (MQ=255)
                                                                        ccGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCg                                                                             >  2:543501/1‑148 (MQ=255)
                                                                        ccGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCg                                                                             >  2:1379778/1‑148 (MQ=255)
                                                                        ccGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCg                                                                             >  2:1380625/1‑148 (MQ=255)
                                                                          gAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGAGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGt                                                                           >  1:1107369/1‑148 (MQ=255)
                                                                          gAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGACCAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAACGCCGATGGCAATTGTAGCCTGATGCCGGCTGCGCGCGGCAATAACGGt                                                                           >  1:559813/1‑148 (MQ=255)
                                                                           aaaCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGtt                                                                          >  2:1717390/1‑148 (MQ=255)
                                                                           aaaCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGtt                                                                          >  1:884983/1‑148 (MQ=255)
                                                                           aaaCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGAGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGtt                                                                          >  1:2898793/1‑148 (MQ=255)
                                                                           aaaCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGtt                                                                          >  1:494589/1‑148 (MQ=255)
                                                                                   gTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGt                                                                  <  1:1565995/148‑1 (MQ=255)
                                                                                   gTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGt                                                                  >  1:1831575/1‑148 (MQ=255)
                                                                                    ttCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTc                                                                 >  1:264060/1‑148 (MQ=255)
                                                                                     tCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTcc                                                                >  2:383977/1‑148 (MQ=255)
                                                                                     tCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTcc                                                                >  2:1339181/1‑148 (MQ=255)
                                                                                          cTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGc                                                           <  2:1278396/148‑1 (MQ=255)
                                                                                                 ttGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACg                                                    >  2:907268/1‑148 (MQ=255)
                                                                                                 ttGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACg                                                    <  1:2219646/148‑1 (MQ=255)
                                                                                                 ttGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACg                                                    >  1:168828/1‑148 (MQ=255)
                                                                                                     ttACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTgcg                                                <  2:1156912/148‑1 (MQ=255)
                                                                                                      tACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTgcgc                                               <  2:462152/148‑1 (MQ=255)
                                                                                                          ggCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGAGCCGGAACGTGCGCTGCa                                           <  2:3208044/148‑1 (MQ=255)
                                                                                                           gCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCaa                                          <  2:2719843/148‑1 (MQ=255)
                                                                                                           gCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCaa                                          >  1:334045/1‑148 (MQ=255)
                                                                                                           gCACAACGTCGAACCGCTGCCGACCTCTATCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCaa                                          >  2:2141776/1‑148 (MQ=255)
                                                                                                              caACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGACCGTGCGCTGCAACTg                                       >  2:3283694/1‑148 (MQ=255)
                                                                                                              caACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTg                                       <  2:334868/148‑1 (MQ=255)
                                                                                                              caACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTg                                       <  2:1987418/148‑1 (MQ=255)
                                                                                                              caACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTg                                       <  2:1988559/148‑1 (MQ=255)
                                                                                                                aCGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTggg                                     >  1:1293374/1‑148 (MQ=255)
                                                                                                                 cGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGt                                    <  1:2061480/148‑1 (MQ=255)
                                                                                                                  gTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTa                                   <  2:678864/148‑1 (MQ=255)
                                                                                                                    cGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAtt                                 >  2:1222554/1‑148 (MQ=255)
                                                                                                                    cGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAtt                                 >  1:1715358/1‑148 (MQ=255)
                                                                                                                    cGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAtt                                 >  1:1849710/1‑148 (MQ=255)
                                                                                                                     gAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTg                                <  2:2589771/148‑1 (MQ=255)
                                                                                                                      aaCCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGAGCCGGAACGTGCGCTGCAACTGGGTATTGa                               >  2:1863134/1‑148 (MQ=255)
                                                                                                                       aCCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGAGCCGGAACGTGCGCTGCAACTGGGTATTGAg                              >  1:2455441/1‑148 (MQ=255)
                                                                                                                        ccGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAgg                             >  2:1894283/1‑148 (MQ=255)
                                                                                                                        ccGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAgg                             <  1:1379778/148‑1 (MQ=255)
                                                                                                                        ccGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAgg                             <  1:1380625/148‑1 (MQ=255)
                                                                                                                        ccGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAgg                             >  1:3083623/1‑148 (MQ=255)
                                                                                                                           cTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCa                          <  2:1544201/148‑1 (MQ=255)
                                                                                                                            tGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAg                         <  2:2993424/148‑1 (MQ=255)
                                                                                                                              ccGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCa                       <  1:2344277/148‑1 (MQ=255)
                                                                                                                                  ccTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACAt                   >  2:2714677/1‑148 (MQ=255)
                                                                                                                                  ccTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACAt                   >  1:1438144/1‑148 (MQ=255)
                                                                                                                                  ccTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACAt                   >  2:1077070/1‑148 (MQ=255)
                                                                                                                                   ctctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATc                  >  2:918107/1‑148 (MQ=255)
                                                                                                                                   ctctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATc                  >  2:3150096/1‑148 (MQ=255)
                                                                                                                                    tctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCa                 >  1:733116/1‑148 (MQ=255)
                                                                                                                                    tctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCa                 >  2:2542868/1‑148 (MQ=255)
                                                                                                                                    tctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCa                 <  2:165357/148‑1 (MQ=255)
                                                                                                                                     ctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCaa                <  1:1688737/148‑1 (MQ=255)
                                                                                                                                        tCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAAccc             <  2:3148253/148‑1 (MQ=255)
                                                                                                                                        tCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAAccc             <  2:1509355/148‑1 (MQ=255)
                                                                                                                                         cTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCg            <  1:2376785/148‑1 (MQ=255)
                                                                                                                                           ggACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGaa          <  2:2864745/148‑1 (MQ=255)
                                                                                                                                            gACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGaaa         <  1:1548077/148‑1 (MQ=255)
                                                                                                                                               tGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAggg      >  2:3077701/1‑148 (MQ=255)
                                                                                                                                                gAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAgggg     >  1:1959451/1‑148 (MQ=255)
                                                                                                                                                gAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAgggg     >  2:277648/1‑148 (MQ=255)
                                                                                                                                                gAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAgggg     <  2:2671384/148‑1 (MQ=255)
                                                                                                                                                   gATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGGTTa  >  1:3145464/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGGTTA  >  CP009273/4024680‑4024974

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: