Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,387,121 C→T 5.8% A833T (GCT→ACT)  yhdP ← DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,387,1210CT 5.8% 1051.0 / 34.3 298A833T (GCT→ACT) yhdPDUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base C (109/169);  new base T (13/4);  total (124/174)
Fisher's exact test for biased strand distribution p-value = 3.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 298 total aligned reads displayed.

CACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACC  >  CP009273/3386974‑3387268
                                                                                                                                                   |                                                                                                                                                   
cACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGTCt                                                                                                                                                   >  2:5011518/1‑150 (MQ=255)
cACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCt                                                                                                                                                   <  1:2582609/150‑1 (MQ=255)
 aCCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTg                                                                                                                                                  <  1:1288526/150‑1 (MQ=255)
    aTTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTcc                                                                                                                                               <  2:2038916/150‑1 (MQ=255)
     ttCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCg                                                                                                                                              >  2:4758326/1‑150 (MQ=255)
      tCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCgg                                                                                                                                             <  2:4546723/150‑1 (MQ=255)
        aTCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACTGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGATTTTGACCGATCAACCAGCGGCTATTGAAATGATTATCCGCACCAGACTGACCgcta                                                                                                                                           >  2:2759933/1‑147 (MQ=255)
        aTCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCATCCTGTCCGGtt                                                                                                                                           >  1:3124187/1‑150 (MQ=255)
           ggcggcATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAAt                                                                                                                                        >  1:4525904/1‑150 (MQ=255)
            gcggcATATTGATTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAAtt                                                                                                                                       <  2:3245981/150‑1 (MQ=255)
            gcggcATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAAtt                                                                                                                                       <  2:2014279/150‑1 (MQ=255)
              ggcATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCCGCGGCTCTTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCa                                                                                                                                     >  2:3645679/1‑150 (MQ=255)
               gcATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCaa                                                                                                                                    <  2:792692/150‑1 (MQ=255)
                 atatTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAg                                                                                                                                  <  1:4334622/150‑1 (MQ=255)
                 atatTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAg                                                                                                                                  >  2:1012426/1‑150 (MQ=255)
                  tatTGAGTTCAACGCCACTTTGTTCCGGCTACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGc                                                                                                                                 <  1:4040369/150‑1 (MQ=255)
                    ttGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTg                                                                                                                               >  2:4616768/1‑150 (MQ=255)
                     tGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGt                                                                                                                              >  2:1821380/1‑150 (MQ=255)
                       aGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTg                                                                                                                            >  2:1854313/1‑150 (MQ=255)
                        gTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTga                                                                                                                           <  1:3606573/150‑1 (MQ=255)
                        gTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGACCGGTTAATTCAAAGCTGTTga                                                                                                                           >  2:2406703/1‑150 (MQ=255)
                          tCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCAGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTgaga                                                                                                                         >  1:1145979/1‑150 (MQ=255)
                            aaCGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGAtt                                                                                                                       <  1:1341388/150‑1 (MQ=255)
                            aaCGCCACTTTGTTCCGGCAACGGCGGGAGCGATTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAAACAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGAtt                                                                                                                       <  2:1100965/150‑1 (MQ=255)
                                ccACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCa                                                                                                                   <  2:2718828/150‑1 (MQ=255)
                                   cTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAACTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCa                                                                                                                <  2:1765840/150‑1 (MQ=255)
                                   cTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCa                                                                                                                <  2:2980847/150‑1 (MQ=255)
                                    tttGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCaa                                                                                                               >  2:3364797/1‑150 (MQ=255)
                                     ttGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAAc                                                                                                              >  1:580126/1‑150 (MQ=255)
                                       gTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCt                                                                                                            <  1:4225393/150‑1 (MQ=255)
                                        ttCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGTCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCtt                                                                                                           <  1:4154668/150‑1 (MQ=255)
                                        ttCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCtt                                                                                                           <  1:3296947/150‑1 (MQ=255)
                                          ccGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAAGAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTaa                                                                                                         <  2:2202828/150‑1 (MQ=255)
                                            ggCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACg                                                                                                       >  1:2097203/1‑150 (MQ=255)
                                            ggCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACg                                                                                                       >  2:4459564/1‑150 (MQ=255)
                                             gCAAAGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCTAGCGTCAGCTTTTTAACGAGCAACCAGCGTCTATTGAAATGATTATCCGCACCATCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGt                                                                                                      <  2:1451302/150‑1 (MQ=255)
                                              cAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGtt                                                                                                     <  2:161110/150‑1 (MQ=255)
                                               aaCGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGttt                                                                                                    <  1:5058006/150‑1 (MQ=255)
                                                   gcggGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCg                                                                                                <  1:3836129/150‑1 (MQ=255)
                                                     gggAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCt                                                                                              >  1:1385921/1‑150 (MQ=255)
                                                     gggAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCt                                                                                              >  1:2263143/1‑150 (MQ=255)
                                                        aGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGTCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGt                                                                                           >  1:1659687/1‑150 (MQ=255)
                                                        aGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGt                                                                                           <  1:2196212/150‑1 (MQ=255)
                                                              ttACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGGTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCa                                                                                     >  2:1377850/1‑150 (MQ=255)
                                                                    tCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTGGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCa                                                                               >  1:3373196/1‑150 (MQ=255)
                                                                      tGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAgg                                                                             <  2:3970481/150‑1 (MQ=255)
                                                                          gccCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCAAGTGGTTCACCCGCAGGTTTg                                                                         >  2:5015592/1‑150 (MQ=255)
                                                                            ccAAAGAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCCGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGc                                                                       <  2:3749439/150‑1 (MQ=255)
                                                                             cAAAGAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCt                                                                      <  1:3597433/150‑1 (MQ=255)
                                                                              aaaTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTa                                                                     >  2:4506050/1‑150 (MQ=255)
                                                                              aaaTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTa                                                                     >  1:4837394/1‑150 (MQ=255)
                                                                              aaaTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTa                                                                     <  1:1193321/150‑1 (MQ=255)
                                                                              aaaTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTa                                                                     <  2:3057258/150‑1 (MQ=255)
                                                                                aTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGTCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAAc                                                                   >  1:4972108/1‑150 (MQ=255)
                                                                                     aCGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGa                                                              >  2:2847296/1‑150 (MQ=255)
                                                                                     aCGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGa                                                              >  2:4937447/1‑150 (MQ=255)
                                                                                          cGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTa                                                         <  1:2189103/150‑1 (MQ=255)
                                                                                             gCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGt                                                      <  1:292060/150‑1 (MQ=255)
                                                                                             gCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCAGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGt                                                      <  1:2031878/150‑1 (MQ=255)
                                                                                              cGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTg                                                     <  1:189281/150‑1 (MQ=255)
                                                                                              cGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTg                                                     <  1:4456330/150‑1 (MQ=255)
                                                                                               gTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGa                                                    >  2:3545695/1‑150 (MQ=255)
                                                                                                 cAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACt                                                  <  2:4461136/150‑1 (MQ=255)
                                                                                                 cAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACt                                                  <  2:4356830/150‑1 (MQ=255)
                                                                                                      tttGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTca                                             >  1:5087676/1‑150 (MQ=255)
                                                                                                         gACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACAt                                          <  2:4890086/150‑1 (MQ=255)
                                                                                                         gACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACAt                                          >  2:3481073/1‑150 (MQ=255)
                                                                                                              aGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCttt                                     <  2:3902758/150‑1 (MQ=255)
                                                                                                                cAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAg                                   >  2:4288592/1‑150 (MQ=255)
                                                                                                                cAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAg                                   >  1:1061604/1‑150 (MQ=255)
                                                                                                                cAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAg                                   >  1:4172541/1‑150 (MQ=255)
                                                                                                                cAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAg                                   <  1:2423976/150‑1 (MQ=255)
                                                                                                                   ccAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATc                                >  2:3079782/1‑150 (MQ=255)
                                                                                                                   ccAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATc                                <  2:4392334/150‑1 (MQ=255)
                                                                                                                    cAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCg                               >  1:81372/1‑150 (MQ=255)
                                                                                                                      gCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGcc                             <  2:2138186/150‑1 (MQ=255)
                                                                                                                      gCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGcc                             <  2:4180308/150‑1 (MQ=255)
                                                                                                                       cGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCg                            <  1:3916049/150‑1 (MQ=255)
                                                                                                                            aatgtcaTGAGTATCCGCCCCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAg                       <  2:3130448/144‑1 (MQ=255)
                                                                                                                               gAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctc                    >  1:1321155/1‑150 (MQ=255)
                                                                                                                                aaaTGATTATCCGCACCAGTCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctct                   <  1:396492/150‑1 (MQ=255)
                                                                                                                                aaaTGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctct                   >  2:4789753/1‑150 (MQ=255)
                                                                                                                                 aaTGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTa                  >  2:4020476/1‑150 (MQ=255)
                                                                                                                                 aaTGATGATCCGCACCAGCCGGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTGGACCGCTAGTGGTTCACCCGCAGGTTTGGCTACCGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTa                  <  1:317500/150‑1 (MQ=255)
                                                                                                                                   tGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATg                <  1:503793/150‑1 (MQ=255)
                                                                                                                                    gATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGt               <  1:1102504/150‑1 (MQ=255)
                                                                                                                                    gATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGt               >  1:5192207/1‑150 (MQ=255)
                                                                                                                                    gATTATCCGCACCAGCCGGTCCGGTTAATTCAAAGCTGTGGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGt               >  1:1587925/1‑150 (MQ=255)
                                                                                                                                       tATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtat            >  1:4580287/1‑150 (MQ=255)
                                                                                                                                       tATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtat            <  1:2381420/150‑1 (MQ=255)
                                                                                                                                        aTCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtata           <  1:2823469/150‑1 (MQ=255)
                                                                                                                                         tCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAg          <  2:2509776/150‑1 (MQ=255)
                                                                                                                                         tCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAg          >  1:4752557/1‑150 (MQ=255)
                                                                                                                                          ccGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAgg         <  1:787497/150‑1 (MQ=255)
                                                                                                                                          ccGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAgg         <  2:5028903/150‑1 (MQ=255)
                                                                                                                                            gCACCAGTCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTc       >  1:2638252/1‑150 (MQ=255)
                                                                                                                                            gCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTc       >  1:830818/1‑150 (MQ=255)
                                                                                                                                                cAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCAc   >  2:1005641/1‑150 (MQ=255)
                                                                                                                                                 aGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCAcc  >  2:2070508/1‑150 (MQ=255)
                                                                                                                                                 aGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCAcc  >  2:700085/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACC  >  CP009273/3386974‑3387268

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: