Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,790 A→T 37.8% L94Q (CTG→CAG)  pepE ← (alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,219,7900AT37.8% 196.8 / 244.9 206L94Q (CTG→CAG) pepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (66/62);  new base T (35/43);  total (101/105)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

CCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAG  >  CP009273/4219643‑4219936
                                                                                                                                                   |                                                                                                                                                  
 cTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAg                                                                                                                                                   >  2:290026/1‑148 (MQ=255)
 cTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAg                                                                                                                                                   >  2:44808/1‑148 (MQ=255)
        aTCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTgg                                                                                                                                            >  2:2396283/1‑148 (MQ=255)
          cGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGaa                                                                                                                                          >  2:2667002/1‑148 (MQ=255)
          cGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGaa                                                                                                                                          >  2:2304/1‑148 (MQ=255)
            aCAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAg                                                                                                                                        <  2:1987928/148‑1 (MQ=255)
            aCAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAg                                                                                                                                        <  1:2049986/148‑1 (MQ=255)
             cAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGt                                                                                                                                       >  2:1045047/1‑148 (MQ=255)
             cAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCACCTGGAAAGt                                                                                                                                       >  2:64720/1‑148 (MQ=255)
             cAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGt                                                                                                                                       <  1:1621120/148‑1 (MQ=255)
              aaTCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTa                                                                                                                                      <  2:3131137/148‑1 (MQ=255)
                 cGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTc                                                                                                                                   >  1:1510151/1‑148 (MQ=255)
                  ggCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTcc                                                                                                                                  <  1:2549719/148‑1 (MQ=255)
                    cATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCccg                                                                                                                                >  2:2536872/1‑148 (MQ=255)
                       atCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCccgccg                                                                                                                             >  1:2722706/1‑148 (MQ=255)
                       atCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCccgccg                                                                                                                             <  2:875874/148‑1 (MQ=255)
                       atCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCccgccg                                                                                                                             >  2:1133621/1‑148 (MQ=255)
                        tCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCccgccgc                                                                                                                            >  2:3533206/1‑148 (MQ=255)
                         cGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCccgccgcc                                                                                                                           >  2:2703712/1‑148 (MQ=255)
                           ttGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGa                                                                                                                         >  2:1209376/1‑148 (MQ=255)
                           ttGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGa                                                                                                                         >  2:2980776/1‑148 (MQ=255)
                              gtggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACaa                                                                                                                      >  2:3343151/1‑148 (MQ=255)
                               tggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAAt                                                                                                                     <  2:3073301/148‑1 (MQ=255)
                               tggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAAt                                                                                                                     <  1:79663/148‑1 (MQ=255)
                                ggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATc                                                                                                                    <  2:701586/148‑1 (MQ=255)
                                ggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATc                                                                                                                    <  2:2761107/148‑1 (MQ=255)
                                 gtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCa                                                                                                                   <  1:2800346/148‑1 (MQ=255)
                                   aCGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCCCg                                                                                                                 >  2:3380162/1‑148 (MQ=255)
                                       aTAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGAtct                                                                                                             >  1:1178839/1‑148 (MQ=255)
                                        tAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGAtctc                                                                                                            >  1:625652/1‑148 (MQ=255)
                                              ggCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAtt                                                                                                      <  2:81972/148‑1 (MQ=255)
                                              ggCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAtt                                                                                                      <  2:1363727/148‑1 (MQ=255)
                                              ggCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAtt                                                                                                      <  1:2304/148‑1 (MQ=255)
                                                caagcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAtttt                                                                                                    <  1:2964967/148‑1 (MQ=255)
                                                caagcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAtttt                                                                                                    <  1:2520445/148‑1 (MQ=255)
                                                caagcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAtttt                                                                                                    <  1:1192036/148‑1 (MQ=255)
                                                    caagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAAt                                                                                                <  1:1859779/148‑1 (MQ=255)
                                                    caagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCCACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAAt                                                                                                >  1:2550568/1‑148 (MQ=255)
                                                    caagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAAt                                                                                                <  2:2032323/148‑1 (MQ=255)
                                                      agGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcg                                                                                              <  2:1709880/148‑1 (MQ=255)
                                                       ggTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgc                                                                                             >  2:2172326/1‑148 (MQ=255)
                                                        gTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgcg                                                                                            >  1:843268/1‑148 (MQ=255)
                                                        gTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgcg                                                                                            >  2:3242870/1‑148 (MQ=255)
                                                         ttAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCgg                                                                                           <  1:205986/148‑1 (MQ=255)
                                                          tAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGc                                                                                          >  1:1064981/1‑148 (MQ=255)
                                                           aGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCa                                                                                         >  1:733087/1‑148 (MQ=255)
                                                            gcgcCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCaa                                                                                        >  1:3218306/1‑148 (MQ=255)
                                                              gcCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACg                                                                                      <  1:1142552/148‑1 (MQ=255)
                                                                cTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACggg                                                                                    <  2:587536/148‑1 (MQ=255)
                                                                 tGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGa                                                                                   >  1:3538631/1‑148 (MQ=255)
                                                                 tGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGa                                                                                   <  2:377885/148‑1 (MQ=255)
                                                                     cTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGa                                                                               <  1:460848/148‑1 (MQ=255)
                                                                        cAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACaa                                                                            >  2:906159/1‑148 (MQ=255)
                                                                            cAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCt                                                                        <  1:367610/148‑1 (MQ=255)
                                                                             aaTATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTa                                                                       >  1:823030/1‑148 (MQ=255)
                                                                               tataCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATg                                                                     <  2:133941/148‑1 (MQ=255)
                                                                                 taCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGaa                                                                   <  1:2472653/148‑1 (MQ=255)
                                                                                  aCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAAt                                                                  <  1:2068661/148‑1 (MQ=255)
                                                                                  aCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAAt                                                                  >  2:2884944/1‑148 (MQ=255)
                                                                                  aCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAAt                                                                  >  2:3515394/1‑148 (MQ=255)
                                                                                       gcgcCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCgg                                                             <  2:2550568/148‑1 (MQ=255)
                                                                                           cACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGAc                                                         >  2:469469/1‑148 (MQ=255)
                                                                                           cACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGAc                                                         >  2:1503576/1‑148 (MQ=255)
                                                                                           cACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGAc                                                         <  1:2211554/148‑1 (MQ=255)
                                                                                            aCGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACa                                                        >  1:3486659/1‑148 (MQ=255)
                                                                                              gTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGa                                                      <  1:2020410/148‑1 (MQ=255)
                                                                                              gTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGa                                                      <  1:2020493/148‑1 (MQ=255)
                                                                                              gTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGa                                                      >  2:3000216/1‑148 (MQ=255)
                                                                                               tttAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGaa                                                     >  1:1412753/1‑148 (MQ=255)
                                                                                                ttAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGaaa                                                    >  1:3121146/1‑148 (MQ=255)
                                                                                                   aCCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTAGTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAacac                                                 <  2:139922/148‑1 (MQ=255)
                                                                                                    ccACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAcc                                                <  1:2218950/148‑1 (MQ=255)
                                                                                                    ccACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAcc                                                >  2:472068/1‑148 (MQ=255)
                                                                                                     cACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAccc                                               >  1:3042981/1‑148 (MQ=255)
                                                                                                      aCGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCAATGAATACCGGTGACAGAAACACCCa                                              >  2:1170530/1‑148 (MQ=255)
                                                                                                          cAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCgg                                          >  2:2691980/1‑148 (MQ=255)
                                                                                                           aGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGa                                         >  2:2128353/1‑148 (MQ=255)
                                                                                                           aGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCACCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGa                                         <  2:1454487/148‑1 (MQ=255)
                                                                                                            gTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAg                                        >  1:2350042/1‑148 (MQ=255)
                                                                                                              aaTTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGGAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCg                                      <  2:675644/148‑1 (MQ=255)
                                                                                                                ttGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCgag                                    >  2:2394504/1‑148 (MQ=255)
                                                                                                                 tGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCgaga                                   <  2:2190151/148‑1 (MQ=255)
                                                                                                                 tGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCgaga                                   <  2:2364520/148‑1 (MQ=255)
                                                                                                                        agcagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAg                            >  2:3424426/1‑148 (MQ=255)
                                                                                                                         gcagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGc                           >  2:1555654/1‑148 (MQ=255)
                                                                                                                          cagcCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                          >  2:2222885/1‑148 (MQ=255)
                                                                                                                             cccGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                       >  1:695173/1‑148 (MQ=255)
                                                                                                                                 agcTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTct                   >  1:3309267/2‑147 (MQ=255)
                                                                                                                                  gcTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGc                  <  2:1230122/148‑1 (MQ=255)
                                                                                                                                   cTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCt                 <  2:3100969/148‑1 (MQ=255)
                                                                                                                                      gcgGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCtgtg              <  1:2186952/148‑1 (MQ=255)
                                                                                                                                      gcgGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCtgtg              >  2:2475964/1‑148 (MQ=255)
                                                                                                                                        ggCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTa            <  1:1310694/148‑1 (MQ=255)
                                                                                                                                         gCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTaa           >  2:1863636/1‑148 (MQ=255)
                                                                                                                                           aCTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatc         >  2:1353433/1‑148 (MQ=255)
                                                                                                                                              gTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatcatc      <  2:1403095/148‑1 (MQ=255)
                                                                                                                                               tttCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc     <  2:1178839/148‑1 (MQ=255)
                                                                                                                                               tttCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc     >  1:469361/1‑148 (MQ=255)
                                                                                                                                               tttCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc     <  1:3533019/148‑1 (MQ=255)
                                                                                                                                                ttCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATccc    <  1:2103966/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAG  >  CP009273/4219643‑4219936

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: ‑