Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 54,392 G→T 100% A104E (GCG→GAG)  surA ← peptidyl‑prolyl cis‑trans isomerase (PPIase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927354,3920GT97.5% 678.3 / 2.0 214A104E (GCG→GAG) surApeptidyl‑prolyl cis‑trans isomerase (PPIase)
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base T (110/98);  total (116/98)
Fisher's exact test for biased strand distribution p-value = 3.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 217 total aligned reads displayed.

AGGATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGC  >  CP009273/54245‑54537
                                                                                                                                                   |                                                                                                                                                 
aGGATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTcgc                                                                                                                                                  <  1:1624819/149‑1 (MQ=255)
 ggATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTctcg                                                                                                                                                 >  2:2300467/1‑149 (MQ=255)
 ggATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTctcg                                                                                                                                                 >  2:2602276/1‑149 (MQ=255)
   aTGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGAt                                                                                                                                               >  1:3982732/1‑149 (MQ=255)
    tGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATg                                                                                                                                              >  2:4103780/1‑149 (MQ=255)
    tGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCAGGTGGGTCTGTTTCTCGATg                                                                                                                                              >  2:3882712/1‑149 (MQ=255)
      gTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGtt                                                                                                                                            <  1:525237/149‑1 (MQ=255)
       tGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTa                                                                                                                                           <  1:3610914/149‑1 (MQ=255)
       tGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTa                                                                                                                                           <  2:2794673/149‑1 (MQ=255)
        gATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAg                                                                                                                                          >  2:1801734/1‑149 (MQ=255)
        gATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAg                                                                                                                                          >  2:3437677/1‑149 (MQ=255)
             gacgacGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATc                                                                                                                                     >  1:551253/1‑149 (MQ=255)
              acgacgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCg                                                                                                                                    >  2:1184494/1‑149 (MQ=255)
                gacgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGcc                                                                                                                                  <  1:2455737/149‑1 (MQ=255)
                gacgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGcc                                                                                                                                  >  2:3402273/1‑149 (MQ=255)
                 acgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCt                                                                                                                                 <  2:3193879/149‑1 (MQ=255)
                     aCCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATc                                                                                                                             >  2:3674974/1‑149 (MQ=255)
                     aCCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATc                                                                                                                             >  2:481919/1‑149 (MQ=255)
                         cGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgc                                                                                                                         <  2:310098/149‑1 (MQ=255)
                         cGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgc                                                                                                                         >  1:3702871/1‑149 (MQ=255)
                         cGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgc                                                                                                                         >  1:3702806/1‑149 (MQ=255)
                                aCGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCAt                                                                                                                  <  1:2051281/149‑1 (MQ=255)
                                    aCTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATTGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGga                                                                                                              <  1:64925/149‑1 (MQ=255)
                                     cTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGgag                                                                                                             <  2:3328519/149‑1 (MQ=255)
                                      ttCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGgaga                                                                                                            >  1:197036/1‑149 (MQ=255)
                                        cAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGAtt                                                                                                          <  1:2271070/149‑1 (MQ=255)
                                         agagATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGAttt                                                                                                         >  1:301000/1‑149 (MQ=255)
                                            gaTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCa                                                                                                      >  2:1628356/1‑149 (MQ=255)
                                            gaTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCa                                                                                                      <  1:2332425/149‑1 (MQ=255)
                                             aTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCAc                                                                                                     <  1:587069/149‑1 (MQ=255)
                                             aTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCAc                                                                                                     <  1:977526/149‑1 (MQ=255)
                                                atcatcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTcc                                                                                                  <  2:2124786/149‑1 (MQ=255)
                                                   atcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAt                                                                                               >  1:688180/1‑149 (MQ=255)
                                                   atcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAt                                                                                               >  2:1293162/1‑149 (MQ=255)
                                                     ctctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAttt                                                                                             <  1:1307409/149‑1 (MQ=255)
                                                      tctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAtttt                                                                                            <  1:1070631/149‑1 (MQ=255)
                                                       ctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGTAGATTTTCACTCCCATTTTc                                                                                           <  2:442793/149‑1 (MQ=255)
                                                         ttGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTg                                                                                         >  2:3709902/1‑149 (MQ=255)
                                                          tGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGc                                                                                        <  2:1164324/149‑1 (MQ=255)
                                                              gATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCa                                                                                    >  2:2024367/1‑149 (MQ=255)
                                                               aTCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCAt                                                                                   >  1:164398/1‑149 (MQ=255)
                                                               aTCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCAt                                                                                   <  1:3170980/149‑1 (MQ=255)
                                                                 cTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCt                                                                                 <  2:1827007/149‑1 (MQ=255)
                                                                 cTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCt                                                                                 >  1:2586190/1‑149 (MQ=255)
                                                                 cTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCt                                                                                 <  2:3604398/149‑1 (MQ=255)
                                                                   ggTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCTTGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGc                                                                               >  2:659188/1‑149 (MQ=255)
                                                                    gTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCa                                                                              >  1:447135/1‑149 (MQ=255)
                                                                          aTAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatga                                                                        >  1:2695578/1‑149 (MQ=255)
                                                                          aTAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTACCATTTTCTGCCCCATCTGCAGgatga                                                                        <  1:1548106/149‑1 (MQ=255)
                                                                           tAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatgat                                                                       >  2:1559376/1‑149 (MQ=255)
                                                                             ggTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGAttt                                                                     >  1:3315055/1‑149 (MQ=255)
                                                                             ggTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGAttt                                                                     >  1:1264211/1‑149 (MQ=255)
                                                                                 ttgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATc                                                                 <  1:2651315/149‑1 (MQ=255)
                                                                                 ttgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATc                                                                 >  1:125142/1‑149 (MQ=255)
                                                                                 ttgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATc                                                                 >  1:2169576/1‑149 (MQ=255)
                                                                                  tgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATcc                                                                <  1:344437/149‑1 (MQ=255)
                                                                                  tgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATcc                                                                <  2:436989/149‑1 (MQ=255)
                                                                                   gtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCa                                                               <  2:502115/149‑1 (MQ=255)
                                                                                    tAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCAt                                                              >  2:1730875/1‑149 (MQ=255)
                                                                                       ttCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTTCCCCATCTGCAGGATGATTTGATCCATGAt                                                           >  1:1243230/1‑149 (MQ=255)
                                                                                       ttCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAt                                                           <  1:1518011/149‑1 (MQ=255)
                                                                                           gTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaaa                                                       >  1:1461019/1‑149 (MQ=255)
                                                                                           gTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaaa                                                       >  1:3735887/1‑149 (MQ=255)
                                                                                           gTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCACATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaaa                                                       >  2:374377/1‑149 (MQ=255)
                                                                                            tCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAAc                                                      >  1:814644/1‑149 (MQ=255)
                                                                                              cATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGt                                                    <  1:841081/149‑1 (MQ=255)
                                                                                              cATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGt                                                    >  2:2045112/1‑149 (MQ=255)
                                                                                                tCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCACTCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTc                                                  >  2:2999248/1‑149 (MQ=255)
                                                                                                tCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTc                                                  >  1:1699076/1‑149 (MQ=255)
                                                                                                 cGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTcc                                                 <  2:3634373/149‑1 (MQ=255)
                                                                                                  gTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCa                                                >  1:1123537/1‑149 (MQ=255)
                                                                                                      gCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGAt                                            >  1:2856008/1‑149 (MQ=255)
                                                                                                         aGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTtg                                         <  2:4059670/149‑1 (MQ=255)
                                                                                                          gACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTtgg                                        >  1:488276/1‑149 (MQ=255)
                                                                                                                 gcgcATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCa                                 <  2:1320667/149‑1 (MQ=255)
                                                                                                                   gcATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgc                               <  1:1184494/149‑1 (MQ=255)
                                                                                                                    cATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcg                              >  2:1619280/1‑149 (MQ=255)
                                                                                                                    cATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcg                              <  1:725286/149‑1 (MQ=255)
                                                                                                                     aTCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgt                             >  1:2514622/1‑149 (MQ=255)
                                                                                                                      tCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtc                            <  2:1752463/149‑1 (MQ=255)
                                                                                                                      tCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtc                            <  2:3410024/149‑1 (MQ=255)
                                                                                                                       cTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcg                           >  2:769848/1‑149 (MQ=255)
                                                                                                                              aGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcatc                    >  2:1260070/1‑149 (MQ=255)
                                                                                                                              aGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcatc                    <  2:3982732/149‑1 (MQ=255)
                                                                                                                               gCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcatca                   >  2:380123/1‑149 (MQ=255)
                                                                                                                                cGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAg                  <  1:3493496/149‑1 (MQ=255)
                                                                                                                                cGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAg                  <  1:1439037/149‑1 (MQ=255)
                                                                                                                                cGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAg                  >  2:829211/1‑149 (MQ=255)
                                                                                                                                  tCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGa                >  1:391621/1‑149 (MQ=255)
                                                                                                                                  tCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCGGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGa                <  2:2501473/149‑1 (MQ=255)
                                                                                                                                    aTGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAg              <  1:1378279/149‑1 (MQ=255)
                                                                                                                                    aTGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAg              >  1:3221629/1‑149 (MQ=255)
                                                                                                                                    aTGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAg              >  2:800705/1‑149 (MQ=255)
                                                                                                                                     tgttgtTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGt             >  1:1341284/1‑149 (MQ=255)
                                                                                                                                     tgttgtTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGt             >  2:2061878/1‑149 (MQ=255)
                                                                                                                                     tgttgtTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGt             <  1:1144216/149‑1 (MQ=255)
                                                                                                                                      gttgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGtt            >  1:2871802/1‑149 (MQ=255)
                                                                                                                                              gTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgcctt    <  2:1401721/149‑1 (MQ=255)
                                                                                                                                               tttCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgccttg   >  1:2409298/1‑149 (MQ=255)
                                                                                                                                               tttCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGAGTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgccttg   >  2:2502601/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGGATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGC  >  CP009273/54245‑54537

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