Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,387,121 C→T 7.0% 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 7.0% 1164.0 / 50.0 327A833T (GCT→ACT) yhdPDUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base C (139/165);  new base T (15/8);  total (154/173)
Fisher's exact test for biased strand distribution p-value = 8.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

CACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCA  >  CP009273/3386974‑3387269
                                                                                                                                                   |                                                                                                                                                    
cACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGcc                                                                                                                                                     <  2:653168/149‑1 (MQ=255)
 aCCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCt                                                                                                                                                    >  2:5121073/1‑149 (MQ=255)
 aCCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCt                                                                                                                                                    >  2:4052583/1‑149 (MQ=255)
  ccATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTg                                                                                                                                                   >  1:2839127/1‑149 (MQ=255)
    aTTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTc                                                                                                                                                 >  1:2222048/1‑149 (MQ=255)
        aTCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCTGGt                                                                                                                                             <  1:1962845/149‑1 (MQ=255)
         tCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGtt                                                                                                                                            <  1:1182695/149‑1 (MQ=255)
         tCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGtt                                                                                                                                            >  1:4708310/1‑149 (MQ=255)
               gcATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCCGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCa                                                                                                                                      >  1:4054563/1‑149 (MQ=255)
                cATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCaa                                                                                                                                     <  1:4771397/149‑1 (MQ=255)
                 atatTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCaaa                                                                                                                                    <  1:1364741/149‑1 (MQ=255)
                  tatTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAg                                                                                                                                   <  1:3680238/149‑1 (MQ=255)
                  tatTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAg                                                                                                                                   >  1:4119239/1‑149 (MQ=255)
                  tatTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAg                                                                                                                                   <  2:1081344/149‑1 (MQ=255)
                   atTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGTCTGTCCGGTTAATTCAAAGc                                                                                                                                  >  2:5143300/1‑149 (MQ=255)
                   atTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGc                                                                                                                                  <  1:1450849/149‑1 (MQ=255)
                    ttGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGTCTGTCCGGTTAATTCAAAGCt                                                                                                                                 <  1:2612513/149‑1 (MQ=255)
                     tGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTg                                                                                                                                <  2:3436223/149‑1 (MQ=255)
                     tGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTg                                                                                                                                <  2:838851/149‑1 (MQ=255)
                      gAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGTAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGt                                                                                                                               >  2:4783900/1‑149 (MQ=255)
                      gAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGt                                                                                                                               <  1:1689877/149‑1 (MQ=255)
                      gAGTTCAACGCCACTTTGTTCCGGCAACGGCGAGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGt                                                                                                                               <  1:3189576/149‑1 (MQ=255)
                         ttCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTga                                                                                                                            <  2:1967124/149‑1 (MQ=255)
                          tCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTgag                                                                                                                           >  1:4316657/1‑149 (MQ=255)
                            aaCGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGAt                                                                                                                         >  2:3821606/1‑149 (MQ=255)
                              cGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTg                                                                                                                       >  2:910385/1‑149 (MQ=255)
                               gCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGc                                                                                                                      <  2:623624/149‑1 (MQ=255)
                               gCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGc                                                                                                                      >  2:4004056/1‑149 (MQ=255)
                                ccACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGcc                                                                                                                     >  2:1278165/1‑149 (MQ=255)
                                  aCTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCAt                                                                                                                   <  2:2496385/149‑1 (MQ=255)
                                   cTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATc                                                                                                                  >  2:863851/1‑149 (MQ=255)
                                   cTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATc                                                                                                                  >  2:562982/1‑149 (MQ=255)
                                   cTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATc                                                                                                                  <  1:5081150/149‑1 (MQ=255)
                                      tGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAAc                                                                                                               <  2:2100881/149‑1 (MQ=255)
                                         tCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCtt                                                                                                            <  1:3334216/149‑1 (MQ=255)
                                               aaCGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGtt                                                                                                      <  2:1856912/149‑1 (MQ=255)
                                               aaCGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGtt                                                                                                      <  2:3345791/149‑1 (MQ=255)
                                                 cggcggGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTa                                                                                                    <  1:2822151/149‑1 (MQ=255)
                                                 cggcggGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTa                                                                                                    <  1:1398314/149‑1 (MQ=255)
                                                     gggAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGc                                                                                                >  2:1383426/1‑149 (MQ=255)
                                                      ggAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCt                                                                                               >  1:3296245/1‑149 (MQ=255)
                                                      ggAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCt                                                                                               >  1:923191/1‑149 (MQ=255)
                                                      ggAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCt                                                                                               <  1:348033/149‑1 (MQ=255)
                                                       gAGCTTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTTTTGAGATTGCCATCAACCTTAACGTTTACCGCTa                                                                                              <  2:594282/149‑1 (MQ=255)
                                                          cGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGTCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTg                                                                                           >  2:5187067/1‑149 (MQ=255)
                                                          cGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTg                                                                                           <  2:4959813/149‑1 (MQ=255)
                                                           gTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTgg                                                                                          <  1:69602/149‑1 (MQ=255)
                                                                 cTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCAcc                                                                                    >  1:5145082/1‑149 (MQ=255)
                                                                  tGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCAccc                                                                                   >  2:2396681/1‑149 (MQ=255)
                                                                  tGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCAccc                                                                                   <  1:637465/149‑1 (MQ=255)
                                                                              aaaTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCt                                                                       <  1:4243041/149‑1 (MQ=255)
                                                                              aaaTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCt                                                                       >  2:3769377/1‑149 (MQ=255)
                                                                               aaTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTa                                                                      <  2:2281322/149‑1 (MQ=255)
                                                                                aTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTaa                                                                     >  1:2756536/1‑149 (MQ=255)
                                                                                aTAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTaa                                                                     <  2:431997/149‑1 (MQ=255)
                                                                                            aGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTaa                                                         <  2:3555638/149‑1 (MQ=255)
                                                                                             gCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAg                                                        <  1:4639072/149‑1 (MQ=255)
                                                                                             gCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAg                                                        >  1:57010/1‑149 (MQ=255)
                                                                                                tCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGa                                                     >  1:4564422/1‑149 (MQ=255)
                                                                                                 cAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGAc                                                    <  2:365585/149‑1 (MQ=255)
                                                                                                   gCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTg                                                  >  2:1468904/1‑149 (MQ=255)
                                                                                                   gCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTg                                                  >  1:745273/1‑149 (MQ=255)
                                                                                                    cTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGc                                                 <  1:219834/149‑1 (MQ=255)
                                                                                                       ttGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGAGTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTca                                              <  2:3205992/149‑1 (MQ=255)
                                                                                                         gACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTcaca                                            <  1:1675750/149‑1 (MQ=255)
                                                                                                         gACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTcaca                                            >  2:4553353/1‑149 (MQ=255)
                                                                                                           ccGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACAtt                                          <  2:699311/149‑1 (MQ=255)
                                                                                                            cGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTc                                         <  1:3013226/149‑1 (MQ=255)
                                                                                                            cGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTACTCACATTc                                         <  1:999576/149‑1 (MQ=255)
                                                                                                              aGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCtt                                       <  1:4946331/149‑1 (MQ=255)
                                                                                                                cAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTa                                     <  2:2813094/149‑1 (MQ=255)
                                                                                                                  aCCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGa                                   >  1:3861636/1‑149 (MQ=255)
                                                                                                                   ccAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGAt                                  <  1:848331/149‑1 (MQ=255)
                                                                                                                      gCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGc                               <  2:2544837/149‑1 (MQ=255)
                                                                                                                      gCGGATATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGc                               <  2:1740539/149‑1 (MQ=255)
                                                                                                                         gCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGt                            <  2:4436782/149‑1 (MQ=255)
                                                                                                                            aTTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCa                         >  2:30677/1‑149 (MQ=255)
                                                                                                                            aTTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCa                         >  2:1695048/1‑149 (MQ=255)
                                                                                                                              tGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGc                       >  1:2307589/1‑149 (MQ=255)
                                                                                                                                 aaTGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctct                    >  1:2481395/1‑149 (MQ=255)
                                                                                                                                  aTGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTa                   >  2:4150682/1‑149 (MQ=255)
                                                                                                                                    gATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATg                 >  1:2020120/1‑149 (MQ=255)
                                                                                                                                        aTCCGCACCAGTCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtat             <  1:1859217/149‑1 (MQ=255)
                                                                                                                                        aTCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtat             >  2:4210744/1‑149 (MQ=255)
                                                                                                                                        aTCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtat             <  2:1656135/149‑1 (MQ=255)
                                                                                                                                         tCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTTAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtata            >  2:2416178/1‑149 (MQ=255)
                                                                                                                                         tCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtata            >  2:1882466/1‑149 (MQ=255)
                                                                                                                                         tCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtata            <  1:2740324/149‑1 (MQ=255)
                                                                                                                                         tCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtata            <  1:1770390/149‑1 (MQ=255)
                                                                                                                                         tCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTtata            >  2:259969/1‑149 (MQ=255)
                                                                                                                                           cGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAgg          <  1:4105847/149‑1 (MQ=255)
                                                                                                                                            gCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGt         >  1:3565394/1‑149 (MQ=255)
                                                                                                                                            gCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGt         <  2:1607035/149‑1 (MQ=255)
                                                                                                                                             cACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTc        >  2:1919985/1‑149 (MQ=255)
                                                                                                                                               ccAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGc      >  1:1030436/1‑149 (MQ=255)
                                                                                                                                               ccAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGATCTATGTTATAGGTCGc      <  1:5043884/149‑1 (MQ=255)
                                                                                                                                                 aGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCAc    >  2:2913365/1‑149 (MQ=255)
                                                                                                                                                  gCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCAcc   >  2:4095905/1‑149 (MQ=255)
                                                                                                                                                  gCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCAcc   <  2:3547203/149‑1 (MQ=255)
                                                                                                                                                   ccTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCa  >  1:3421823/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCA  >  CP009273/3386974‑3387269

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: