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,3920GT99.0% 653.1 / ‑2.0 202A104E (GCG→GAG) surApeptidyl‑prolyl cis‑trans isomerase (PPIase)
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (102/98);  total (102/100)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 203 total aligned reads displayed.

GATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTG  >  CP009273/54247‑54536
                                                                                                                                                 |                                                                                                                                                
 aTGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGa                                                                                                                                               >  2:2932754/1‑148 (MQ=255)
  tGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGGTCTGTTTCTCGAt                                                                                                                                              >  2:1696839/1‑148 (MQ=255)
    gTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGt                                                                                                                                            <  1:472475/148‑1 (MQ=255)
    gTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGt                                                                                                                                            <  1:746779/148‑1 (MQ=255)
         gCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCa                                                                                                                                       <  2:3180105/148‑1 (MQ=255)
         gCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCa                                                                                                                                       <  2:1498162/148‑1 (MQ=255)
            acgacgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATc                                                                                                                                    >  2:2032448/1‑148 (MQ=255)
              gacgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGc                                                                                                                                  <  2:1193585/148‑1 (MQ=255)
               acgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGcc                                                                                                                                 <  1:605121/148‑1 (MQ=255)
                   aCCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGAt                                                                                                                             <  2:2431127/148‑1 (MQ=255)
                    ccTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATc                                                                                                                            <  2:2196733/148‑1 (MQ=255)
                       cGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAg                                                                                                                         <  2:2041596/148‑1 (MQ=255)
                       cGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAg                                                                                                                         >  1:2159046/1‑148 (MQ=255)
                       cGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAg                                                                                                                         <  2:47306/148‑1 (MQ=255)
                        gttgttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgc                                                                                                                        <  2:1962404/148‑1 (MQ=255)
                         ttgttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgca                                                                                                                       >  2:3359853/1‑147 (MQ=255)
                          tgttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgctg                                                                                                                      >  1:1143251/1‑148 (MQ=255)
                          tgttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgctg                                                                                                                      >  2:2051347/1‑148 (MQ=255)
                          tgttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgctg                                                                                                                      <  2:3326277/148‑1 (MQ=255)
                             tACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTc                                                                                                                   <  2:617050/148‑1 (MQ=255)
                               cGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCAt                                                                                                                 >  1:2786270/1‑148 (MQ=255)
                                gCACGTCAGAGATAATCATCTCTTGTCGGATCTGGTTACGATAGGTGTGGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCTTGTTGTTCTGTGTCTCTATGTTAGCAATCGCCTGATACAGCTGCTCATc                                                                                                                <  1:1841646/148‑1 (MQ=255)
                                 cACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCg                                                                                                               >  2:2111501/1‑148 (MQ=255)
                                   cTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGga                                                                                                             >  1:1305571/1‑148 (MQ=255)
                                     tCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGgaga                                                                                                           >  2:2727686/1‑148 (MQ=255)
                                      cAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGAt                                                                                                          <  1:303613/148‑1 (MQ=255)
                                        gagaTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGAttt                                                                                                        >  2:2500653/1‑148 (MQ=255)
                                        gagaTAATCATCACTTTGCGGATCAGCTTACGATAGGTGTTCTAGTTCAGTCCATCGTAAGCCACACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATCTAAGCAAACGCCAGATCCAGCTGCTCATCGGAGAttt                                                                                                        >  2:543426/1‑148 (MQ=255)
                                          gaTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTc                                                                                                      <  2:2340179/148‑1 (MQ=255)
                                           aTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCa                                                                                                     <  1:2568208/148‑1 (MQ=255)
                                             aaTCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACt                                                                                                   >  1:197572/1‑148 (MQ=255)
                                             aaTCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACt                                                                                                   >  1:1981623/1‑148 (MQ=255)
                                               tcatcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTcc                                                                                                 <  2:1151929/148‑1 (MQ=255)
                                                catcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTccc                                                                                                <  2:846194/148‑1 (MQ=255)
                                                  tctcTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAt                                                                                              <  1:3194412/148‑1 (MQ=255)
                                                  tctcTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAt                                                                                              >  1:82922/1‑148 (MQ=255)
                                                     ctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAtttt                                                                                           <  1:3123570/148‑1 (MQ=255)
                                                               cTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATc                                                                                 <  2:900387/148‑1 (MQ=255)
                                                               cTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATc                                                                                 >  2:127634/1‑148 (MQ=255)
                                                               cTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATc                                                                                 <  2:328137/148‑1 (MQ=255)
                                                                   ttACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCa                                                                             >  2:1939853/1‑148 (MQ=255)
                                                                   ttACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCa                                                                             >  2:1940577/1‑148 (MQ=255)
                                                                     aCGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAgg                                                                           <  2:1596532/148‑1 (MQ=255)
                                                                        aTAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatg                                                                        <  2:532674/148‑1 (MQ=255)
                                                                        aTAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGAACCAGCGTCATGTTGTTCTGTATCTCTATGTTAGCAATCTCCTGATCCAGCTGCTCATCGGAGATTTTAACTCCCATATTCTGCCCCATCTGCAGgatg                                                                        >  1:1470639/1‑148 (MQ=255)
                                                                          aGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatgat                                                                      >  1:306060/1‑148 (MQ=255)
                                                                          aGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatgat                                                                      >  2:2421555/1‑148 (MQ=255)
                                                                          aGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCAATCTGCAGgatgat                                                                      >  1:1863073/1‑148 (MQ=255)
                                                                           ggTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGAtt                                                                     <  1:2218217/148‑1 (MQ=255)
                                                                           ggTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGAtt                                                                     >  1:623438/1‑148 (MQ=255)
                                                                                tgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATc                                                                <  1:127634/148‑1 (MQ=255)
                                                                                 gtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATcc                                                               <  2:2636046/148‑1 (MQ=255)
                                                                                 gtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCACATTTTCTGCCCCATCTGCAGGATGATTTGATcc                                                               <  1:1556090/148‑1 (MQ=255)
                                                                                  tAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCACTCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCa                                                              <  2:1851175/148‑1 (MQ=255)
                                                                                  tAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCa                                                              >  1:1721186/1‑148 (MQ=255)
                                                                                  tAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCa                                                              >  2:40845/1‑148 (MQ=255)
                                                                                     ttCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGa                                                           >  2:2974431/1‑148 (MQ=255)
                                                                                     ttCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGa                                                           <  2:1968086/148‑1 (MQ=255)
                                                                                     ttCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGa                                                           <  1:3309836/148‑1 (MQ=255)
                                                                                      tCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAt                                                          <  2:3432146/148‑1 (MQ=255)
                                                                                       cAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATc                                                         <  1:2943844/148‑1 (MQ=255)
                                                                                       cAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATc                                                         <  1:2946455/148‑1 (MQ=255)
                                                                                        aGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCa                                                        >  2:2370959/1‑148 (MQ=255)
                                                                                        aGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCa                                                        <  2:2585961/148‑1 (MQ=255)
                                                                                         gTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaa                                                       >  2:3439266/1‑148 (MQ=255)
                                                                                          tCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaaa                                                      >  1:980623/1‑148 (MQ=255)
                                                                                          tCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaaa                                                      >  1:3468315/1‑148 (MQ=255)
                                                                                          tCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaaa                                                      >  2:88228/1‑148 (MQ=255)
                                                                                          tCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaaa                                                      <  1:2500653/148‑1 (MQ=255)
                                                                                          tCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaaa                                                      >  1:206211/1‑148 (MQ=255)
                                                                                           ccATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAAc                                                     >  2:3285704/1‑148 (MQ=255)
                                                                                            cATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACg                                                    >  2:570321/1‑148 (MQ=255)
                                                                                             aTCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGt                                                   <  2:619955/148‑1 (MQ=255)
                                                                                              tCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGtt                                                  >  2:3032430/1‑148 (MQ=255)
                                                                                                 tAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCa                                               <  2:2858567/148‑1 (MQ=255)
                                                                                                  aaGCCAGACGGCTGCGCATCTGATCGAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCAt                                              >  2:837529/1‑148 (MQ=255)
                                                                                                  aaGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCAt                                              >  1:1471045/1‑148 (MQ=255)
                                                                                                     ccAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGAt                                           <  2:1839315/148‑1 (MQ=255)
                                                                                                       aGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGAttt                                         <  2:1467554/148‑1 (MQ=255)
                                                                                                                cgcATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCa                                >  1:2314369/1‑148 (MQ=255)
                                                                                                                  cATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgc                              >  2:502312/1‑148 (MQ=255)
                                                                                                                  cATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgc                              >  2:766566/1‑148 (MQ=255)
                                                                                                                    tCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgt                            >  2:1450953/1‑148 (MQ=255)
                                                                                                                        aTCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgcg                        >  2:790638/1‑148 (MQ=255)
                                                                                                                          ccAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgcgtc                      >  1:2993321/1‑148 (MQ=255)
                                                                                                                          ccAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgcgtc                      <  2:2148310/148‑1 (MQ=255)
                                                                                                                          ccAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgcgtc                      >  1:3160661/1‑148 (MQ=255)
                                                                                                                            aGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcat                    >  2:3207244/1‑148 (MQ=255)
                                                                                                                            aGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcat                    >  2:3206583/1‑148 (MQ=255)
                                                                                                                              cGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcatca                  <  2:1738695/148‑1 (MQ=255)
                                                                                                                              cGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTAGCGtcatca                  >  1:368626/1‑148 (MQ=255)
                                                                                                                                     ttgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGtt           >  1:135828/1‑148 (MQ=255)
                                                                                                                                     ttgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGtt           >  1:2392345/1‑148 (MQ=255)
                                                                                                                                      tgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtg          >  1:1070755/1‑148 (MQ=255)
                                                                                                                                      tgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtg          <  2:2159046/148‑1 (MQ=255)
                                                                                                                                       gttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgc         <  1:2932754/148‑1 (MQ=255)
                                                                                                                                            gTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgcct    >  1:2628340/1‑148 (MQ=255)
                                                                                                                                            gTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgcct    >  1:3012931/1‑148 (MQ=255)
                                                                                                                                             tttCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgcctt   >  1:504313/1‑148 (MQ=255)
                                                                                                                                              ttCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgccttg  <  1:2555699/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTG  >  CP009273/54247‑54536

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: