Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,736,722 C→A 100% G210G (GGC→GGA yiaK → 2,3‑diketo‑L‑gulonate reductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,736,7220CA99.5% 657.3 / ‑4.3 214G210G (GGC→GGAyiaK2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (116/96);  total (116/98)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 214 total aligned reads displayed.

CCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGC  >  CP009273/3736577‑3736866
                                                                                                                                                 |                                                                                                                                                
ccATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGt                                                                                                                                                <  2:1188760/148‑1 (MQ=255)
 cATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGACCCTGGAGtt                                                                                                                                               >  2:3184644/1‑148 (MQ=255)
  aTTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTa                                                                                                                                              <  2:1422561/148‑1 (MQ=255)
  aTTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGAGGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTa                                                                                                                                              >  2:25733/1‑148 (MQ=255)
   ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGCTGGGGGCTTTGATGATGAGGGCACTTTGACCAAAGAACCTGGAGTTAt                                                                                                                                             >  2:2423713/1‑148 (MQ=255)
   ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTAt                                                                                                                                             <  2:3752914/148‑1 (MQ=255)
   ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTAt                                                                                                                                             <  1:1163361/148‑1 (MQ=255)
   ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTAt                                                                                                                                             <  1:2239399/148‑1 (MQ=255)
   ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTAt                                                                                                                                             >  2:501646/1‑148 (MQ=255)
   ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTAt                                                                                                                                             <  1:1690417/148‑1 (MQ=255)
   ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTAt                                                                                                                                             <  1:1690434/148‑1 (MQ=255)
   ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTAt                                                                                                                                             <  1:1690542/148‑1 (MQ=255)
   ttccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTAt                                                                                                                                             >  2:2155098/1‑148 (MQ=255)
      cttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCga                                                                                                                                          >  2:3128617/1‑148 (MQ=255)
       ttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGGTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATagag                                                                                                                                         >  1:58703/1‑148 (MQ=255)
       ttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCgag                                                                                                                                         >  1:3459826/1‑148 (MQ=255)
            cGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAgaagaa                                                                                                                                    <  1:1314374/148‑1 (MQ=255)
                 aTCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATcgcc                                                                                                                               >  1:3774949/1‑148 (MQ=255)
                  tCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATcgccg                                                                                                                              >  2:2320649/1‑148 (MQ=255)
                  tCACAATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATAGAGAAGAATcgccg                                                                                                                              >  2:460905/1‑148 (MQ=255)
                     ccATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAt                                                                                                                           >  2:1536830/1‑148 (MQ=255)
                     ccATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAt                                                                                                                           >  2:1742388/1‑148 (MQ=255)
                       aTGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAttt                                                                                                                         <  1:145039/148‑1 (MQ=255)
                       aTGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAttt                                                                                                                         <  2:3598109/148‑1 (MQ=255)
                       aTGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAttt                                                                                                                         <  1:1255098/148‑1 (MQ=255)
                              atatGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCCATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCCCCGCAATTTGCCGAt                                                                                                                  >  1:3241987/1‑148 (MQ=255)
                              atatGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGAt                                                                                                                  >  1:3561569/1‑148 (MQ=255)
                              atatGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGAt                                                                                                                  >  2:1634381/1‑148 (MQ=255)
                                atgtcgatgtcgATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATgg                                                                                                                >  1:2942507/1‑148 (MQ=255)
                                atgtcgatgtcgATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATgg                                                                                                                >  1:1031329/1‑148 (MQ=255)
                                      atgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACt                                                                                                          >  2:759069/1‑148 (MQ=255)
                                          cgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGaa                                                                                                      <  1:3112652/148‑1 (MQ=255)
                                             tgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAgg                                                                                                   >  2:2895740/1‑148 (MQ=255)
                                              gtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGt                                                                                                  <  1:3334178/148‑1 (MQ=255)
                                                tctcTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGACCCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTc                                                                                                <  1:2999511/148‑1 (MQ=255)
                                                tctcTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGCATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTc                                                                                                <  1:3156685/148‑1 (MQ=255)
                                                tctcTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTc                                                                                                >  1:858412/1‑148 (MQ=255)
                                                tctcTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTc                                                                                                >  2:1883002/1‑148 (MQ=255)
                                                  tctTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTg                                                                                              >  1:2805733/1‑148 (MQ=255)
                                                     tACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCa                                                                                           <  2:3684256/148‑1 (MQ=255)
                                                         gCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTc                                                                                       <  1:3423326/148‑1 (MQ=255)
                                                           aTGTTAGATGGTCACCGTCTGGCAGGTCGTCCGCTCGCGGTGGATGGTGGCTTTGAGGATGCGGGCTATTTGACCAACGAACGTGGCGTTATCGAGAAGAATGGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGa                                                                                     <  2:2503651/148‑1 (MQ=255)
                                                              ttAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATtg                                                                                  >  2:284884/1‑148 (MQ=255)
                                                              ttAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATtg                                                                                  >  2:3794997/1‑148 (MQ=255)
                                                              ttAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATtg                                                                                  >  2:402515/1‑148 (MQ=255)
                                                                    gTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTgg                                                                            <  2:1713992/148‑1 (MQ=255)
                                                                    gTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTgg                                                                            <  2:1713026/148‑1 (MQ=255)
                                                                    gTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTgg                                                                            <  1:1899970/148‑1 (MQ=255)
                                                                    gTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGGGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTgg                                                                            <  2:1489741/148‑1 (MQ=255)
                                                                     ttAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTTCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGa                                                                           >  1:2388146/1‑148 (MQ=255)
                                                                     ttAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGa                                                                           >  2:1307173/1‑148 (MQ=255)
                                                                      tAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGat                                                                          >  2:3661873/1‑148 (MQ=255)
                                                                      tAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGat                                                                          >  2:1622212/1‑148 (MQ=255)
                                                                         ccGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATg                                                                       <  2:3398585/148‑1 (MQ=255)
                                                                         ccGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATg                                                                       >  1:2076462/1‑148 (MQ=255)
                                                                         ccGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATg                                                                       <  2:1832063/148‑1 (MQ=255)
                                                                          cGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGa                                                                      <  2:621452/148‑1 (MQ=255)
                                                                            tCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCGACTGGAAAGGTTCTGGCCTGTCGATTGTGCTGGTTATGATc                                                                    >  1:59450/1‑148 (MQ=255)
                                                                                gCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTa                                                                <  1:1946889/148‑1 (MQ=255)
                                                                                gCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTa                                                                <  1:18279/148‑1 (MQ=255)
                                                                                gCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTa                                                                <  1:506879/148‑1 (MQ=255)
                                                                                   ggTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGTTGAGGGCAATTTGACCAAAAAACCTGGAGTTTTCGAGAAGAATCGCCGCATTTTTCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTActc                                                             >  1:719325/1‑148 (MQ=255)
                                                                                     tcgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACtctc                                                           <  1:919363/148‑1 (MQ=255)
                                                                                     tcgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACtctc                                                           <  1:2460707/148‑1 (MQ=255)
                                                                                     tcgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACtctc                                                           <  1:1164787/148‑1 (MQ=255)
                                                                                     tcgtcAGCTCCCGGCCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACtctc                                                           <  2:1113246/148‑1 (MQ=255)
                                                                                      cgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTcc                                                          <  2:1981565/148‑1 (MQ=255)
                                                                                      cgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTcc                                                          <  2:1981570/148‑1 (MQ=255)
                                                                                        tcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCtt                                                        <  1:713954/148‑1 (MQ=255)
                                                                                        tcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCtt                                                        <  1:2629926/148‑1 (MQ=255)
                                                                                                   tCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGgcgc                                             <  2:899754/148‑1 (MQ=255)
                                                                                                      aTGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATc                                          >  1:980082/1‑148 (MQ=255)
                                                                                                        ggtggCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCg                                        <  1:3499801/148‑1 (MQ=255)
                                                                                                          tggCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGtt                                      >  1:3602880/1‑148 (MQ=255)
                                                                                                           ggCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTg                                     <  2:836550/148‑1 (MQ=255)
                                                                                                           ggCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTg                                     <  2:298350/148‑1 (MQ=255)
                                                                                                               ttGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGa                                 >  1:1112669/1‑148 (MQ=255)
                                                                                                               ttGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGa                                 >  1:2733778/1‑148 (MQ=255)
                                                                                                                   tgatgaGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTc                             <  1:17955/148‑1 (MQ=255)
                                                                                                                   tgatgaGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTc                             <  1:240640/148‑1 (MQ=255)
                                                                                                                          ggCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAgg                      >  2:2038107/1‑148 (MQ=255)
                                                                                                                          ggCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAgg                      >  1:1315142/1‑148 (MQ=255)
                                                                                                                          ggCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAgg                      >  2:108249/1‑148 (MQ=255)
                                                                                                                          ggCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAgg                      >  2:2300573/1‑148 (MQ=255)
                                                                                                                             aaTTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCAGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGaca                   <  1:2423713/148‑1 (MQ=255)
                                                                                                                              aTTTGACCAAAGAAGCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGTCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTGCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaa                  <  2:1555088/148‑1 (MQ=255)
                                                                                                                              aTTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaa                  <  2:3459826/148‑1 (MQ=255)
                                                                                                                              aTTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaa                  <  1:3095825/148‑1 (MQ=255)
                                                                                                                              aTTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaa                  <  2:58703/148‑1 (MQ=255)
                                                                                                                                ttGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaaca                >  2:2093505/1‑148 (MQ=255)
                                                                                                                                   aCCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcg             <  1:2204187/148‑1 (MQ=255)
                                                                                                                                    ccAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcga            <  1:3128617/148‑1 (MQ=255)
                                                                                                                                      aaaGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcgacg          >  1:2766018/1‑148 (MQ=255)
                                                                                                                                      aaaGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcgacg          >  1:2710046/1‑148 (MQ=255)
                                                                                                                                       aaGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCACATCCGTTGCCGAAGTCACCCAGGACAACAGcgacga         >  1:3719041/1‑148 (MQ=255)
                                                                                                                                        aGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGaa        >  2:12781/1‑148 (MQ=255)
                                                                                                                                         taaCGTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGTATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCTAAGTCACCCAGGACAACAGCGACGAAt       <  1:407485/147‑1 (MQ=255)
                                                                                                                                         gAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAAt       <  1:473792/148‑1 (MQ=255)
                                                                                                                                          aaCCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATa      <  2:1363995/148‑1 (MQ=255)
                                                                                                                                              tGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGc  >  1:2520373/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGC  >  CP009273/3736577‑3736866

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