Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,122,013 G→T 100% T119T (ACC→ACA wcaF ← putative acyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,122,0130GT98.8% 806.6 / 0.5 246T119T (ACC→ACAwcaFputative acyl transferase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (103/140);  total (105/141)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-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 248 total aligned reads displayed.

GGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAAT  >  CP009273/2121868‑2122157
                                                                                                                                                 |                                                                                                                                                
 gCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATtgtg                                                                                                                                               >  1:2533477/1‑148 (MQ=255)
    aaCCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGaaa                                                                                                                                            <  2:3122897/148‑1 (MQ=255)
    aaCCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGaaa                                                                                                                                            <  1:34487/148‑1 (MQ=255)
     aCCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAAt                                                                                                                                           >  2:1628603/1‑148 (MQ=255)
      ccACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATg                                                                                                                                          <  2:2500045/148‑1 (MQ=255)
       cacaTTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGt                                                                                                                                         >  2:1778201/1‑148 (MQ=255)
       cacaTTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGt                                                                                                                                         <  1:569575/148‑1 (MQ=255)
         caTTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTg                                                                                                                                       >  1:1293641/1‑148 (MQ=255)
         caTTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTg                                                                                                                                       <  1:3339009/148‑1 (MQ=255)
            ttGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACt                                                                                                                                    >  2:3449557/1‑148 (MQ=255)
             tGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCCCGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACACTAGGCGTGGCGTTAATTGTGAAATGTTGACtt                                                                                                                                   <  1:3488834/148‑1 (MQ=255)
                  gAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATg                                                                                                                              <  2:1415506/148‑1 (MQ=255)
                  gAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATg                                                                                                                              >  1:2117179/1‑148 (MQ=255)
                   aaGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATgg                                                                                                                             >  1:453129/1‑148 (MQ=255)
                   aaGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATgg                                                                                                                             >  2:2991969/1‑148 (MQ=255)
                   aaGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATgg                                                                                                                             <  2:3397517/148‑1 (MQ=255)
                     gCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTTTGAAATGTTGACTTGCATGGTc                                                                                                                           >  1:1453129/1‑148 (MQ=255)
                     gCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTc                                                                                                                           >  2:707553/1‑148 (MQ=255)
                       gATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGt                                                                                                                         <  1:471389/148‑1 (MQ=255)
                        aTTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTg                                                                                                                        <  2:623454/148‑1 (MQ=255)
                        aTTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTg                                                                                                                        >  1:1928011/1‑148 (MQ=255)
                         tttAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTgg                                                                                                                       <  1:2299648/148‑1 (MQ=255)
                          ttAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGc                                                                                                                      >  1:2403916/1‑148 (MQ=255)
                             aaaaCACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTAc                                                                                                                   >  1:3319799/1‑148 (MQ=255)
                              aaaCACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACACTAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTAcc                                                                                                                  >  2:1088673/1‑148 (MQ=255)
                              aaaCACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTAcc                                                                                                                  >  2:3625807/1‑148 (MQ=255)
                              aaaCACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTAcc                                                                                                                  >  2:603430/1‑148 (MQ=255)
                               aaCACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCg                                                                                                                 >  1:1780038/1‑148 (MQ=255)
                                acacTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCgg                                                                                                                <  1:3604672/148‑1 (MQ=255)
                                acacTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGGTGACTTGCATGGTCGTGGCTACCgg                                                                                                                >  2:747620/1‑148 (MQ=255)
                                   cTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGc                                                                                                             >  2:1322731/1‑148 (MQ=255)
                                   cTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGc                                                                                                             >  2:3237638/1‑148 (MQ=255)
                                     gCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCat                                                                                                           >  1:2691114/1‑148 (MQ=255)
                                       ttCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCataa                                                                                                         <  2:2512580/148‑1 (MQ=255)
                                        tCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCataaa                                                                                                        >  1:1964110/1‑148 (MQ=255)
                                        tCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCataaa                                                                                                        <  2:2233076/148‑1 (MQ=255)
                                        tCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCataaa                                                                                                        <  2:838961/148‑1 (MQ=255)
                                          gTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCataaata                                                                                                      <  2:2383877/148‑1 (MQ=255)
                                             cACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACt                                                                                                   <  1:489248/148‑1 (MQ=255)
                                             cACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACt                                                                                                   <  2:3760237/148‑1 (MQ=255)
                                               cccACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACttt                                                                                                 <  1:3015577/148‑1 (MQ=255)
                                                ccACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACtttt                                                                                                <  2:3652564/148‑1 (MQ=255)
                                                ccACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACtttt                                                                                                <  2:135749/148‑1 (MQ=255)
                                                 cACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACttttt                                                                                               >  2:2858313/1‑148 (MQ=255)
                                                  acgacgGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACtttttt                                                                                              >  1:2620279/1‑148 (MQ=255)
                                                   cgacgGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGTTCGTGGCTACCGGTGCATAAATAACTTTTTTg                                                                                             <  2:12392/148‑1 (MQ=255)
                                                     acgGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCg                                                                                           >  2:1476279/1‑148 (MQ=255)
                                                     acgGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCg                                                                                           <  1:2382054/148‑1 (MQ=255)
                                                     acgGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGACAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCg                                                                                           <  1:1433118/148‑1 (MQ=255)
                                                        gTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCACTCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGata                                                                                        <  1:1026477/148‑1 (MQ=255)
                                                           ccGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCa                                                                                     <  1:3518368/148‑1 (MQ=255)
                                                               cGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGa                                                                                 <  2:3251862/148‑1 (MQ=255)
                                                                      tgtgACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTgcg                                                                          >  2:2498814/1‑148 (MQ=255)
                                                                         gACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCa                                                                       >  1:2011188/1‑148 (MQ=255)
                                                                         gACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCa                                                                       <  1:681137/148‑1 (MQ=255)
                                                                          aCGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCaa                                                                      <  1:1532337/148‑1 (MQ=255)
                                                                          aCGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCaa                                                                      >  2:2184956/1‑148 (MQ=255)
                                                                           cGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAAt                                                                     <  2:3746097/148‑1 (MQ=255)
                                                                             ccTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATgg                                                                   <  1:1559351/148‑1 (MQ=255)
                                                                                 gggCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTAt                                                               <  2:1216883/148‑1 (MQ=255)
                                                                                 gggCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTAt                                                               <  1:126965/148‑1 (MQ=255)
                                                                                    caacaaAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTc                                                            >  2:1726480/1‑148 (MQ=255)
                                                                                       caaAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCAcc                                                         >  2:3101544/1‑148 (MQ=255)
                                                                                        aaaGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCg                                                        <  2:1844380/148‑1 (MQ=255)
                                                                                        aaaGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCg                                                        <  2:481933/148‑1 (MQ=255)
                                                                                        aaaGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCg                                                        >  2:2263945/1‑148 (MQ=255)
                                                                                        aaaGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCg                                                        <  2:1844371/148‑1 (MQ=255)
                                                                                        aaaGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCg                                                        <  2:2620279/148‑1 (MQ=255)
                                                                                          aGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAg                                                      >  2:2383006/1‑148 (MQ=255)
                                                                                          aGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAg                                                      >  2:1356123/1‑148 (MQ=255)
                                                                                            aCATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAggg                                                    <  1:86512/148‑1 (MQ=255)
                                                                                            aCATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAggg                                                    >  1:2463508/1‑148 (MQ=255)
                                                                                             cATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGt                                                   <  1:3769680/148‑1 (MQ=255)
                                                                                             cATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGt                                                   <  1:2392442/148‑1 (MQ=255)
                                                                                               tCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGtat                                                 >  2:76538/1‑148 (MQ=255)
                                                                                                   ttGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATaa                                             <  1:784929/148‑1 (MQ=255)
                                                                                                   ttGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATaa                                             >  2:3409735/1‑148 (MQ=255)
                                                                                                     gccagccaGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAAt                                           <  2:890475/148‑1 (MQ=255)
                                                                                                       cagccagcATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTg                                         <  1:1084956/148‑1 (MQ=255)
                                                                                                        agccagcATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGa                                        <  1:3202230/148‑1 (MQ=255)
                                                                                                        agccagcATGTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGa                                        <  2:3574942/148‑1 (MQ=255)
                                                                                                            agcATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTc                                    >  1:2946812/1‑148 (MQ=255)
                                                                                                             gcATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCa                                   <  2:2161092/148‑1 (MQ=255)
                                                                                                              cATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCAt                                  <  2:1315627/148‑1 (MQ=255)
                                                                                                              cATTTCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCAt                                  <  2:1964110/148‑1 (MQ=255)
                                                                                                                tttCTCGCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCg                                <  1:3638026/148‑1 (MQ=255)
                                                                                                                      gCCAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGAcc                          <  1:354457/148‑1 (MQ=255)
                                                                                                                       ccAATCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGAccc                         >  2:659472/1‑148 (MQ=255)
                                                                                                                          aTCACAATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCAcg                      >  2:1612491/1‑148 (MQ=255)
                                                                                                                             acaATAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGt                   >  1:2341479/1‑148 (MQ=255)
                                                                                                                                 tAGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATc               >  1:963910/1‑148 (MQ=255)
                                                                                                                                  aGGCGTGGCGTTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCa              <  1:2590543/148‑1 (MQ=255)
                                                                                                                                    gcgtggcgtTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCAcc            <  2:20937/148‑1 (MQ=255)
                                                                                                                                      gtggcgtTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTa          <  1:1088673/148‑1 (MQ=255)
                                                                                                                                        ggcgtTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAg        <  1:603430/148‑1 (MQ=255)
                                                                                                                                         gcgtTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAgg       <  1:3427921/148‑1 (MQ=255)
                                                                                                                                         gcgtTAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAgg       <  1:2991969/148‑1 (MQ=255)
                                                                                                                                             tAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTaa   >  2:352208/1‑148 (MQ=255)
                                                                                                                                             tAATTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTaa   <  1:1747750/148‑1 (MQ=255)
                                                                                                                                              aaTTGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAAt  <  2:2691114/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAATAGGCGTGGCGTTAATGGTGAAATGTTGACTTGCATGGTCGTGGCTACCGGTGCATAAATAACTTTTTTGCGATATCACCGAATGTGCGCCAATGGTTATTTCACCGAGGGTATATAAATTGACGTCATCGCCGACCCACGCGTAATCACCTAAGGTTAAT  >  CP009273/2121868‑2122157

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: