Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,538,254 C→A 13.2% P5Q (CCG→CAG)  gntX → DNA catabolic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,538,2540CA13.2% 1053.7 / 99.9 319P5Q (CCG→CAG) gntXDNA catabolic protein
Reads supporting (aligned to +/- strand):  ref base C (131/146);  new base A (25/17);  total (156/163)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 328 total aligned reads displayed.

GTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTCCCATCTTCCCTGCGGTCGCTGCC  >  CP009273/3538105‑3538403
                                                                                                                                                     |                                                                                                                                                     
gTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTAcc                                                                                                                                                       <  2:273651/150‑1 (MQ=255)
gTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTAcc                                                                                                                                                       <  1:1679959/150‑1 (MQ=255)
 tCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCg                                                                                                                                                      >  1:1373838/1‑150 (MQ=255)
 tCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCg                                                                                                                                                      >  1:3966955/1‑150 (MQ=255)
 tCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCg                                                                                                                                                      >  2:4062714/1‑150 (MQ=255)
  cAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCgg                                                                                                                                                     <  2:3192742/150‑1 (MQ=255)
  cAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCgg                                                                                                                                                     <  2:3865484/150‑1 (MQ=255)
   aaTGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCggg                                                                                                                                                    <  2:995745/150‑1 (MQ=255)
   aaTGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCggg                                                                                                                                                    <  1:4060791/150‑1 (MQ=255)
   aaTGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCggg                                                                                                                                                    >  1:1344792/1‑150 (MQ=255)
    aTGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGa                                                                                                                                                   >  2:3717874/1‑150 (MQ=255)
     tGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGAt                                                                                                                                                  <  1:2865054/150‑1 (MQ=255)
      gCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGAtt                                                                                                                                                 >  2:1604456/1‑150 (MQ=255)
      gCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGAtt                                                                                                                                                 <  2:1524241/150‑1 (MQ=255)
        aaCGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTAt                                                                                                                                               <  2:2015882/150‑1 (MQ=255)
           gCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCt                                                                                                                                            >  1:1947016/1‑150 (MQ=255)
              acacTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGc                                                                                                                                         >  2:2826024/1‑150 (MQ=255)
               cacTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCt                                                                                                                                        <  2:1342504/150‑1 (MQ=255)
                 cTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTAt                                                                                                                                      <  2:3534164/150‑1 (MQ=255)
                  ttCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATg                                                                                                                                     <  1:2618355/150‑1 (MQ=255)
                    cGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGcc                                                                                                                                   <  2:3326316/150‑1 (MQ=255)
                     ggCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCg                                                                                                                                  <  1:1568183/150‑1 (MQ=255)
                        aTTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAAt                                                                                                                               >  1:4396749/1‑150 (MQ=255)
                        aTTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAAt                                                                                                                               <  1:3162510/150‑1 (MQ=255)
                        aTTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACAGGGATTATGCTGGCTATGCCGAAt                                                                                                                               <  2:1186415/150‑1 (MQ=255)
                          tCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGc                                                                                                                             >  1:661631/1‑150 (MQ=255)
                             gTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCAc                                                                                                                          >  1:3649328/1‑150 (MQ=255)
                                     ccGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTTGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTa                                                                                                                  <  2:3411839/150‑1 (MQ=255)
                                      cGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAg                                                                                                                 <  2:3180698/150‑1 (MQ=255)
                                      cGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACAGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAg                                                                                                                 >  1:799074/1‑150 (MQ=255)
                                          cagcagcaCAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACAGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCa                                                                                                             >  2:2123633/1‑150 (MQ=255)
                                          cagcagcaCAAGATGAACAGTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCa                                                                                                             <  2:539029/150‑1 (MQ=255)
                                            gcagcaCAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCAtt                                                                                                           <  2:2911682/150‑1 (MQ=255)
                                              agcaCAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTgg                                                                                                         >  1:624516/1‑150 (MQ=255)
                                              agcaCAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTgg                                                                                                         <  2:3435827/150‑1 (MQ=255)
                                              agcaCAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACAGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTgg                                                                                                         <  2:4352092/150‑1 (MQ=255)
                                               gcaCAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTggg                                                                                                        >  2:3140380/1‑150 (MQ=255)
                                               gcaCAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTggg                                                                                                        >  2:634588/1‑150 (MQ=255)
                                                 acaAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTggggg                                                                                                      <  2:867912/150‑1 (MQ=255)
                                                   aaGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGAt                                                                                                    >  2:4147198/1‑150 (MQ=255)
                                                    aGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGAtt                                                                                                   >  2:597974/1‑150 (MQ=255)
                                                    aGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGAtt                                                                                                   <  2:1873843/150‑1 (MQ=255)
                                                       tGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACAGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGt                                                                                                <  1:2548911/150‑1 (MQ=255)
                                                          aCATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCg                                                                                             <  2:3343597/150‑1 (MQ=255)
                                                          aCATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCg                                                                                             <  2:4233196/150‑1 (MQ=255)
                                                          aCATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCg                                                                                             >  1:24217/1‑150 (MQ=255)
                                                           cATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCgg                                                                                            >  2:555243/1‑150 (MQ=255)
                                                            aTTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCAGTGGGGGATTTGTTCGGt                                                                                           >  2:2145524/1‑150 (MQ=255)
                                                             ttCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTc                                                                                          >  1:2665865/1‑150 (MQ=255)
                                                               ccccTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTg                                                                                        <  1:1455969/150‑1 (MQ=255)
                                                               ccccTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTg                                                                                        <  1:141199/150‑1 (MQ=255)
                                                                    gACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCAc                                                                                   <  1:428293/150‑1 (MQ=255)
                                                                     aCCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCAcg                                                                                  <  1:2269133/150‑1 (MQ=255)
                                                                     aCCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCAcg                                                                                  <  1:497379/150‑1 (MQ=255)
                                                                      ccTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCAcgc                                                                                 >  1:3365032/1‑150 (MQ=255)
                                                                        tttGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACAGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACgcgc                                                                               >  2:2161775/1‑150 (MQ=255)
                                                                         ttGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGcc                                                                              <  2:399693/150‑1 (MQ=255)
                                                                              cTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCg                                                                         <  1:4257104/150‑1 (MQ=255)
                                                                               tGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGc                                                                        <  2:1362312/150‑1 (MQ=255)
                                                                                  caccaGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGcaca                                                                     <  2:2935814/150‑1 (MQ=255)
                                                                                  caccaGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGcaca                                                                     <  2:673132/150‑1 (MQ=255)
                                                                                     caGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGAt                                                                  <  1:3581730/150‑1 (MQ=255)
                                                                                     caGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGAt                                                                  <  1:1023750/150‑1 (MQ=255)
                                                                                      aGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATa                                                                 <  2:1437580/150‑1 (MQ=255)
                                                                                       gATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATaa                                                                <  2:2612280/150‑1 (MQ=255)
                                                                                        aTGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATaaa                                                               <  1:701134/150‑1 (MQ=255)
                                                                                        aTGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATaaa                                                               <  1:3359058/150‑1 (MQ=255)
                                                                                          gTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAAc                                                             <  1:2003521/150‑1 (MQ=255)
                                                                                          gTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAAc                                                             >  1:4209612/1‑150 (MQ=255)
                                                                                          gTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACAGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCAATGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAAc                                                             >  2:1282938/1‑150 (MQ=255)
                                                                                           ttATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACg                                                            <  1:111264/150‑1 (MQ=255)
                                                                                           ttATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACg                                                            <  1:111629/150‑1 (MQ=255)
                                                                                             aTTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGtt                                                          >  2:1731356/1‑150 (MQ=255)
                                                                                               tCATCCGCTATTGTACTATTTTGACTAACAAGGATGAACATAAGCTAACAGTAACGGGATTAAGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGGTCGGTCAGTTCACGCGCCAACCGCACAGAAAAAActgaat                                                        >  2:2027713/1‑145 (MQ=255)
                                                                                                  tCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGcc                                                     <  2:2428244/150‑1 (MQ=255)
                                                                                                       tATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACaa                                                >  1:2993195/1‑150 (MQ=255)
                                                                                                        aTTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAAt                                               <  1:1266352/150‑1 (MQ=255)
                                                                                                        aTTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAAt                                               >  2:1027492/1‑150 (MQ=255)
                                                                                                        aTTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAAt                                               >  2:1685353/1‑150 (MQ=255)
                                                                                                          tGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAAtgt                                             <  2:3365633/150‑1 (MQ=255)
                                                                                                           gTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAAtgtg                                            <  1:2624298/150‑1 (MQ=255)
                                                                                                            ttCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTgg                                           >  1:4096019/1‑150 (MQ=255)
                                                                                                                    tGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCg                                   <  1:2804499/150‑1 (MQ=255)
                                                                                                                     gACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGc                                  >  2:3134488/1‑150 (MQ=255)
                                                                                                                     gACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGc                                  >  1:1842095/1‑150 (MQ=255)
                                                                                                                      aCTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGcc                                 >  2:1950996/1‑150 (MQ=255)
                                                                                                                       cTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCca                                >  1:3330311/1‑150 (MQ=255)
                                                                                                                          aCAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCcacac                             <  1:3674810/150‑1 (MQ=255)
                                                                                                                          aCAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCcacac                             >  1:2013012/1‑150 (MQ=255)
                                                                                                                             aGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACt                          >  2:763649/1‑150 (MQ=255)
                                                                                                                             aGGATGAACATATGCTAACAGGACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACt                          <  1:2640383/150‑1 (MQ=255)
                                                                                                                                aTGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTccc                       >  1:3737126/1‑150 (MQ=255)
                                                                                                                                aTGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTccc                       <  1:1145313/150‑1 (MQ=255)
                                                                                                                                     cATATGCTAACAGTACAGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTCCCATCtt                  >  2:3803525/1‑150 (MQ=255)
                                                                                                                                        atGCTAACAGTACAGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCATACACTCCCATCTTccc               <  2:92716/150‑1 (MQ=255)
                                                                                                                                          gCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTCCCATCTTCCCTg             >  1:931165/1‑150 (MQ=255)
                                                                                                                                                aGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTCCCATCTTCCCTGCGGTCg       >  1:57863/1‑150 (MQ=255)
                                                                                                                                                aGTACAGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTCCCATCTTCCCTGCGGTCg       >  2:453156/1‑150 (MQ=255)
                                                                                                                                                     cGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTCCCATCTTCCCTGCGGTCGctgcc  >  2:1020747/1‑150 (MQ=255)
                                                                                                                                                     cGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTCCCATCTTCCCTGCGGTCGctgcc  >  2:2356438/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTCATCCGCTATTGTTCTCTTTTGACTTACAAGGATGAACATATGCTAACAGTACCGGGATTATGCTGGCTATGCCGAATGCCACTGGCGTTAGGTCATTGGGGGATTTGTTCGGTCTGTTCACGCGCCACCCGCACAGATAAAACGTTATGCCCACAATGTGGATTACCCGCCACACACTCCCATCTTCCCTGCGGTCGCTGCC  >  CP009273/3538105‑3538403

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: