Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,024,827 G→T 100% D103Y (GAT→TAT)  pepQ → proline dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,024,8270GT99.2% 866.2 / ‑2.1 272D103Y (GAT→TAT) pepQproline dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (120/149);  total (122/150)
Fisher's exact test for biased strand distribution p-value = 5.89e-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 272 total aligned reads displayed.

CAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGG  >  CP009273/4024683‑4024970
                                                                                                                                                |                                                                                                                                               
cAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtg                                                                                                                                              <  2:246516/148‑1 (MQ=255)
 aaTTCAAACCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtgt                                                                                                                                             >  1:3074315/1‑148 (MQ=255)
  aTTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtgtg                                                                                                                                            >  2:1458751/1‑148 (MQ=255)
   ttCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTgg                                                                                                                                           <  2:3488396/148‑1 (MQ=255)
   ttCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTgg                                                                                                                                           <  2:2811976/148‑1 (MQ=255)
    tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                          >  1:2369145/1‑148 (MQ=255)
    tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                          <  1:2123518/148‑1 (MQ=255)
       aaGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGt                                                                                                                                       >  1:2906294/1‑148 (MQ=255)
       aaGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGt                                                                                                                                       >  2:2499190/1‑148 (MQ=255)
       aaGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGt                                                                                                                                       <  2:1594944/148‑1 (MQ=255)
        aGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTg                                                                                                                                      >  1:3271697/1‑148 (MQ=255)
        aGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTg                                                                                                                                      >  1:3439879/1‑148 (MQ=255)
        aGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTg                                                                                                                                      <  2:622016/148‑1 (MQ=255)
        aGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTg                                                                                                                                      >  2:1964462/1‑148 (MQ=255)
         gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGa                                                                                                                                     >  2:1693576/1‑148 (MQ=255)
         gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGa                                                                                                                                     <  1:3562328/148‑1 (MQ=255)
          cGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGAt                                                                                                                                    >  1:3706386/1‑148 (MQ=255)
          cGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGAt                                                                                                                                    <  2:3628897/148‑1 (MQ=255)
                 gCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTg                                                                                                                             <  2:3091967/148‑1 (MQ=255)
                      tAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGaa                                                                                                                        >  2:2725232/1‑148 (MQ=255)
                       aaCTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGaaa                                                                                                                       <  2:2375392/148‑1 (MQ=255)
                        aCTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAg                                                                                                                      <  1:3473469/148‑1 (MQ=255)
                         cTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGc                                                                                                                     <  1:807829/148‑1 (MQ=255)
                            aGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGa                                                                                                                  <  1:2754563/148‑1 (MQ=255)
                             ggTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGAt                                                                                                                 <  2:2547842/148‑1 (MQ=255)
                             ggTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGAt                                                                                                                 <  2:3634124/148‑1 (MQ=255)
                              gTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATg                                                                                                                <  2:494021/148‑1 (MQ=255)
                                gCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTAATGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGc                                                                                                              <  1:2411817/148‑1 (MQ=255)
                                     aCTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTgg                                                                                                         >  2:2070749/1‑148 (MQ=255)
                                     aCTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTgg                                                                                                         <  2:3448687/148‑1 (MQ=255)
                                      ctgctgGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGt                                                                                                        <  1:116099/148‑1 (MQ=255)
                                       tgctgggtgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTa                                                                                                       <  2:847321/148‑1 (MQ=255)
                                            gttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctg                                                                                                  <  1:3767355/148‑1 (MQ=255)
                                             ttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgc                                                                                                 >  2:3090917/1‑148 (MQ=255)
                                             ttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgc                                                                                                 <  1:2899620/148‑1 (MQ=255)
                                                ctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgc                                                                                              >  1:3080941/1‑148 (MQ=255)
                                                  ggtggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcct                                                                                            >  1:2985431/1‑148 (MQ=255)
                                                   gtggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctg                                                                                           <  1:1209964/148‑1 (MQ=255)
                                                           cGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcgcgc                                                                                   >  1:1588006/1‑148 (MQ=255)
                                                             tGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCgg                                                                                 >  1:180488/1‑148 (MQ=255)
                                                                 cAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAat                                                                             <  1:1034133/148‑1 (MQ=255)
                                                                   aGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAatat                                                                           <  1:395247/148‑1 (MQ=255)
                                                                    gCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGACTGCTGCGCGCGGCAATATc                                                                          >  1:2345336/1‑148 (MQ=255)
                                                                     ccGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCg                                                                         >  1:667430/1‑148 (MQ=255)
                                                                        aaaCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGtt                                                                      >  2:1395254/1‑148 (MQ=255)
                                                                        aaaCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGtt                                                                      >  2:3300612/1‑148 (MQ=255)
                                                                         aaCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTta                                                                     <  2:2843796/148‑1 (MQ=255)
                                                                          aCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTtat                                                                    <  1:1796761/148‑1 (MQ=255)
                                                                           cTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTtata                                                                   >  2:2220868/1‑148 (MQ=255)
                                                                            tgtgGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTatat                                                                  <  2:3194280/148‑1 (MQ=255)
                                                                              tgGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCg                                                                >  1:488997/1‑148 (MQ=255)
                                                                                gTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGt                                                              <  1:1415508/148‑1 (MQ=255)
                                                                                  tCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCACTATCGGTTATATCGGTcc                                                            >  2:678787/1‑148 (MQ=255)
                                                                                  tCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTcc                                                            >  2:2130516/1‑148 (MQ=255)
                                                                                  tCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTcc                                                            >  2:389679/1‑148 (MQ=255)
                                                                                    tATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCgg                                                          <  1:3264599/148‑1 (MQ=255)
                                                                                    tATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCgg                                                          <  1:1361382/148‑1 (MQ=255)
                                                                                      tCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTg                                                        <  1:1483136/148‑1 (MQ=255)
                                                                                          ccGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCgg                                                    >  1:2691412/1‑148 (MQ=255)
                                                                                          ccGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCgg                                                    >  1:1060001/1‑148 (MQ=255)
                                                                                           cGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGa                                                   <  2:1839236/148‑1 (MQ=255)
                                                                                           cGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGa                                                   <  2:121597/148‑1 (MQ=255)
                                                                                           cGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGa                                                   <  1:931286/148‑1 (MQ=255)
                                                                                              ttGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACg                                                <  1:1557540/148‑1 (MQ=255)
                                                                                              ttGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACg                                                >  2:1099083/1‑148 (MQ=255)
                                                                                               tGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGt                                               <  2:3360219/148‑1 (MQ=255)
                                                                                               tGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGt                                               >  1:2762754/1‑148 (MQ=255)
                                                                                                gATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTg                                              <  2:1500103/148‑1 (MQ=255)
                                                                                                     cTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTg                                         >  2:2265910/1‑148 (MQ=255)
                                                                                                       ggCACCACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCa                                       <  2:2084937/148‑1 (MQ=255)
                                                                                                       ggCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCa                                       >  1:1991944/1‑148 (MQ=255)
                                                                                                        gCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCaa                                      <  2:1755368/148‑1 (MQ=255)
                                                                                                           caACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTg                                   <  1:666584/148‑1 (MQ=255)
                                                                                                            aaCGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTgg                                  <  2:1910805/148‑1 (MQ=255)
                                                                                                              cGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGt                                >  1:1623527/1‑148 (MQ=255)
                                                                                                              cGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGt                                <  2:3012469/148‑1 (MQ=255)
                                                                                                                tCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAt                              <  1:3112471/148‑1 (MQ=255)
                                                                                                                tCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAt                              <  2:1574474/148‑1 (MQ=255)
                                                                                                                   aaCCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGa                           >  2:3438252/1‑148 (MQ=255)
                                                                                                                    aCCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAg                          >  1:150184/1‑148 (MQ=255)
                                                                                                                    aCCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAg                          >  1:1299030/1‑148 (MQ=255)
                                                                                                                     ccGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCACCTGGGTATTGAgg                         >  2:397639/1‑148 (MQ=255)
                                                                                                                      cGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGc                        >  1:3478589/1‑148 (MQ=255)
                                                                                                                           ccGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCa                   <  1:3146340/148‑1 (MQ=255)
                                                                                                                           ccGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCa                   <  1:1283026/148‑1 (MQ=255)
                                                                                                                            cGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCaa                  <  2:3713814/148‑1 (MQ=255)
                                                                                                                              aCCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACa                >  2:598103/1‑148 (MQ=255)
                                                                                                                               ccTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACAt               >  1:3544431/1‑148 (MQ=255)
                                                                                                                                 tctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCa             <  2:780709/148‑1 (MQ=255)
                                                                                                                                  ctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCaa            >  2:3524018/1‑148 (MQ=255)
                                                                                                                                  ctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCaa            >  2:3244003/1‑148 (MQ=255)
                                                                                                                                  ctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCaa            <  1:1964462/148‑1 (MQ=255)
                                                                                                                                  ctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAGTGAGGCCAGCAACATCaa            >  2:2079948/1‑148 (MQ=255)
                                                                                                                                   tgtCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAAc           <  2:3131753/146‑1 (MQ=255)
                                                                                                                                         gACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGaaa     <  1:991993/148‑1 (MQ=255)
                                                                                                                                         gACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGaaa     <  1:2707728/148‑1 (MQ=255)
                                                                                                                                           cTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAgg   <  2:2226608/148‑1 (MQ=255)
                                                                                                                                           cTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAgg   <  1:1703418/148‑1 (MQ=255)
                                                                                                                                           cTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAgg   <  1:1053333/148‑1 (MQ=255)
                                                                                                                                            tGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAggg  >  1:2749644/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGG  >  CP009273/4024683‑4024970

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: