Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 411,446 C→T 5.9% Q321Q (CAG→CAA sbcD ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273411,4460CT 5.9% 1135.2 / 39.4 305Q321Q (CAG→CAAsbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (145/142);  new base T (9/9);  total (154/151)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 316 total aligned reads displayed.

GCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTT  >  CP009273/411297‑411595
                                                                                                                                                     |                                                                                                                                                     
gCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTgcgc                                                                                                                                                       >  2:3854286/1‑150 (MQ=255)
     aCGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAAAGCCTGGATTTTGCGTTGAat                                                                                                                                                  >  2:4939563/1‑150 (MQ=255)
       ggCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGTTGAatat                                                                                                                                                >  2:5046132/1‑150 (MQ=255)
          gATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCAt                                                                                                                                             <  2:2465231/150‑1 (MQ=255)
          gATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCAt                                                                                                                                             >  2:22347/1‑150 (MQ=255)
          gATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCAt                                                                                                                                             >  1:2530945/1‑150 (MQ=255)
                acacCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGAt                                                                                                                                       <  2:1565159/150‑1 (MQ=255)
                 cacCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATa                                                                                                                                      >  2:4702600/1‑150 (MQ=255)
                  acCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGTTGAATATCATGCAGATAc                                                                                                                                     <  2:1727875/150‑1 (MQ=255)
                   ccTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACt                                                                                                                                    <  1:3015354/150‑1 (MQ=255)
                     tctTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtca                                                                                                                                  >  2:4022222/1‑150 (MQ=255)
                      ctTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtcat                                                                                                                                 >  2:3213873/1‑150 (MQ=255)
                       ttCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtcatc                                                                                                                                <  1:1355630/150‑1 (MQ=255)
                          gACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagt                                                                                                                             >  1:4797220/1‑150 (MQ=255)
                          gACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagt                                                                                                                             >  2:3491710/1‑150 (MQ=255)
                           aCGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagta                                                                                                                            <  2:4279381/150‑1 (MQ=255)
                             gCTGAGTTCGCTGAGGTTTTCACGCTGTTGGCTGGCTAACACACGCTCTCGCAGTTCACTACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTATGCGCTGAATATAATGCAGATACTCATCagtagt                                                                                                                          >  2:149849/1‑150 (MQ=255)
                              cTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTg                                                                                                                         <  2:1926161/150‑1 (MQ=255)
                              cTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTg                                                                                                                         >  1:862661/1‑150 (MQ=255)
                                   ttCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTc                                                                                                                    <  1:2895575/150‑1 (MQ=255)
                                   ttCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTc                                                                                                                    <  1:365403/150‑1 (MQ=255)
                                      gCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGat                                                                                                                 >  2:2043366/1‑150 (MQ=255)
                                        tGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGatat                                                                                                               <  1:2849636/150‑1 (MQ=255)
                                        tGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGatat                                                                                                               >  1:3807553/1‑150 (MQ=255)
                                           gggTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATcca                                                                                                            >  2:2968545/1‑150 (MQ=255)
                                               ttCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATccagcca                                                                                                        >  2:1313442/1‑150 (MQ=255)
                                               ttCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATccagcca                                                                                                        <  1:1018337/150‑1 (MQ=255)
                                                      tGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGt                                                                                                 <  1:1969386/150‑1 (MQ=255)
                                                       gTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTg                                                                                                <  2:2878684/150‑1 (MQ=255)
                                                          ggctggctAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCt                                                                                             >  1:422100/1‑150 (MQ=255)
                                                            ctggctAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTcc                                                                                           <  1:5138695/150‑1 (MQ=255)
                                                             tggctAACACTCTCTCGCGCTGTTCACGACTCCTACGTACCAGCAATACTTCGACAGTCAATGATTCGGTTAATGCCTGTATTTTGCTCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCt                                                                                          <  2:3295067/150‑1 (MQ=255)
                                                             tggctAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCt                                                                                          <  2:659776/150‑1 (MQ=255)
                                                             tggctAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCt                                                                                          <  2:2257136/150‑1 (MQ=255)
                                                              ggctAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTg                                                                                         <  1:1108051/150‑1 (MQ=255)
                                                               gctAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGc                                                                                        <  1:3213873/150‑1 (MQ=255)
                                                                  aaCACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGAt                                                                                     >  1:4372822/1‑150 (MQ=255)
                                                                  aaCACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGAt                                                                                     >  1:1820306/1‑150 (MQ=255)
                                                                    cacacGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATAc                                                                                   >  2:4831841/1‑150 (MQ=255)
                                                                     acacGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACa                                                                                  >  1:3750601/1‑150 (MQ=255)
                                                                       acGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATc                                                                                <  1:4938586/150‑1 (MQ=255)
                                                                       acGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATc                                                                                >  2:476833/1‑150 (MQ=255)
                                                                          cTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATcgcg                                                                             >  2:4599024/1‑150 (MQ=255)
                                                                            cgcgCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGcc                                                                           >  2:81863/1‑150 (MQ=255)
                                                                            cgcgCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGcc                                                                           <  2:3368793/150‑1 (MQ=255)
                                                                            cgcgCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGcc                                                                           <  1:1999364/150‑1 (MQ=255)
                                                                             gcgcTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCa                                                                          <  1:4991230/150‑1 (MQ=255)
                                                                             gcgcTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCa                                                                          <  1:3962093/150‑1 (MQ=255)
                                                                               gcTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACt                                                                        <  2:740506/150‑1 (MQ=255)
                                                                               gcTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACt                                                                        >  2:1000910/1‑150 (MQ=255)
                                                                                cTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTg                                                                       >  2:4537092/1‑150 (MQ=255)
                                                                                 tGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGt                                                                      >  1:1450268/1‑150 (MQ=255)
                                                                                  gTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGtt                                                                     >  2:2533433/1‑150 (MQ=255)
                                                                                  gTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGtt                                                                     >  2:4368068/1‑150 (MQ=255)
                                                                                      aCGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAg                                                                 >  1:4041268/1‑150 (MQ=255)
                                                                                        gACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCt                                                               >  2:3134002/1‑150 (MQ=255)
                                                                                        gACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCt                                                               <  1:3413847/150‑1 (MQ=255)
                                                                                         aCTCCGACGTACCAGCAATACTTCGTCAGGCTATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCTGATACTCATCTGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGAGATACATCGCGCCACTGTTACAGCtg                                                              <  1:4492881/150‑1 (MQ=255)
                                                                                           tCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCtgtg                                                            >  1:3358426/1‑150 (MQ=255)
                                                                                            ccGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGc                                                           <  2:4250466/150‑1 (MQ=255)
                                                                                               aCGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGt                                                        >  1:948883/1‑150 (MQ=255)
                                                                                                   aCCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATc                                                    >  1:2408165/1‑150 (MQ=255)
                                                                                                   aCCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATc                                                    <  2:3707278/150‑1 (MQ=255)
                                                                                                   aCCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATc                                                    <  2:2357944/150‑1 (MQ=255)
                                                                                                   aCCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATc                                                    >  2:3620699/1‑150 (MQ=255)
                                                                                                   aCCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATc                                                    >  2:1377318/1‑150 (MQ=255)
                                                                                                   aCCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATc                                                    >  1:3501749/1‑150 (MQ=255)
                                                                                                    ccAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCg                                                   >  1:3101058/1‑150 (MQ=255)
                                                                                                      aGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGAc                                                 <  2:2254518/150‑1 (MQ=255)
                                                                                                       gCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACg                                                >  2:51655/1‑150 (MQ=255)
                                                                                                          aTACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCCGACAGGTGGCTCCTGCGCTCCATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCa                                             >  2:391840/1‑150 (MQ=255)
                                                                                                          aTACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCa                                             >  1:614102/1‑150 (MQ=255)
                                                                                                                 gACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGcc                                      >  2:4175938/1‑150 (MQ=255)
                                                                                                                  aCAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCt                                     <  2:1217799/150‑1 (MQ=255)
                                                                                                                  aCAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCt                                     <  1:3652435/150‑1 (MQ=255)
                                                                                                                   cAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCtt                                    >  2:882493/1‑150 (MQ=255)
                                                                                                                   cAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCtt                                    <  2:311517/150‑1 (MQ=255)
                                                                                                                        aaTGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGc                               >  1:4559725/1‑150 (MQ=255)
                                                                                                                        aaTGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGc                               >  1:1513242/1‑150 (MQ=255)
                                                                                                                          tGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCAc                             <  2:2434251/150‑1 (MQ=255)
                                                                                                                           gATTCGGTTAATGCCTGGATTTTGCTCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGAAATCCAGCCAGACAGGTGGCTCCTGCGATACATAGCGCCACTGTTCCAGCTGAGCGGTAATCGACGCCAGATCGCCTTTCAGCACt                            >  1:848148/1‑150 (MQ=255)
                                                                                                                              tCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGcc                         <  2:2300425/150‑1 (MQ=255)
                                                                                                                              tCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGcc                         <  1:4581411/150‑1 (MQ=255)
                                                                                                                                ggTTAATGCCTGGATTTTGCGTTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCAt                       >  2:2759597/1‑150 (MQ=255)
                                                                                                                                ggTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCAt                       <  1:761963/150‑1 (MQ=255)
                                                                                                                                     aTGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGGTATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGtt                  <  1:3622002/150‑1 (MQ=255)
                                                                                                                                     aTGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGtt                  <  1:627716/150‑1 (MQ=255)
                                                                                                                                     aTGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGtt                  <  1:3707047/150‑1 (MQ=255)
                                                                                                                                        ccTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCg               <  1:696610/150‑1 (MQ=255)
                                                                                                                                          tGGATTTTGCGTTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGtt             >  2:1829545/1‑150 (MQ=255)
                                                                                                                                          tGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGtt             <  2:2651066/150‑1 (MQ=255)
                                                                                                                                          tGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGtt             >  1:1251949/1‑150 (MQ=255)
                                                                                                                                           ggATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTa            >  1:2509930/1‑150 (MQ=255)
                                                                                                                                            gATTTTGCGTTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTAc           <  1:2285298/150‑1 (MQ=255)
                                                                                                                                             aTTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTAcc          <  1:2152633/150‑1 (MQ=255)
                                                                                                                                                ttGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGt       <  2:3977324/150‑1 (MQ=255)
                                                                                                                                                 tGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTa      <  2:2585103/150‑1 (MQ=255)
                                                                                                                                                  gcgcTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTAc     >  1:2058981/1‑150 (MQ=255)
                                                                                                                                                  gcgcTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTAc     >  1:265396/1‑150 (MQ=255)
                                                                                                                                                    gcTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGt   <  1:1667785/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTGGCTAACACACGCTCGCGCTGTTCACGACTCCGACGTACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTT  >  CP009273/411297‑411595

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: