Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,920,050 C→A 100% I128I (ATC→ATA ygdH → UPF0717 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,920,0500CA96.4% 630.8 / 12.3 210I128I (ATC→ATAygdHUPF0717 family protein
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (95/107);  total (99/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 212 total aligned reads displayed.

GAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCT  >  CP009273/2919907‑2920193
                                                                                                                                               |                                                                                                                                               
gAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTa                                                                                                                                             <  1:1374831/148‑1 (MQ=255)
  agagGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACAtaactaac                                                                                                                                           >  1:2892206/1‑148 (MQ=255)
    agGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCt                                                                                                                                         >  2:1501714/1‑148 (MQ=255)
    agGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCt                                                                                                                                         <  1:2517747/148‑1 (MQ=255)
    agGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCt                                                                                                                                         >  1:1722194/1‑148 (MQ=255)
     ggCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTg                                                                                                                                        <  1:2675338/148‑1 (MQ=255)
      gCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTgg                                                                                                                                       <  2:3030631/148‑1 (MQ=255)
        ttttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTc                                                                                                                                     <  1:2985132/148‑1 (MQ=255)
        ttttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTc                                                                                                                                     <  1:1311977/148‑1 (MQ=255)
        ttttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTc                                                                                                                                     <  1:936583/148‑1 (MQ=255)
        ttttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTc                                                                                                                                     >  1:1758906/1‑148 (MQ=255)
         tttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCt                                                                                                                                    <  2:2794060/148‑1 (MQ=255)
         tttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCt                                                                                                                                    <  2:1013376/148‑1 (MQ=255)
         tttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCt                                                                                                                                    <  2:1446633/148‑1 (MQ=255)
         tttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACAAACCTGGTCt                                                                                                                                    >  2:2578767/1‑148 (MQ=255)
          ttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCtt                                                                                                                                   >  1:2925417/1‑148 (MQ=255)
            gTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCtttt                                                                                                                                 <  1:292020/148‑1 (MQ=255)
            gTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCtttt                                                                                                                                 <  1:2133622/148‑1 (MQ=255)
                       aaTTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGt                                                                                                                      >  1:2980590/1‑148 (MQ=255)
                       aaTTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGt                                                                                                                      >  1:3509631/1‑148 (MQ=255)
                        attattCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTa                                                                                                                     >  1:2060089/1‑148 (MQ=255)
                         ttattCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTaa                                                                                                                    <  2:3191197/148‑1 (MQ=255)
                           attCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACg                                                                                                                  >  1:756679/1‑148 (MQ=255)
                             tCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCt                                                                                                                >  2:3484266/1‑148 (MQ=255)
                                  cTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTcgcgc                                                                                                           >  1:1130997/1‑148 (MQ=255)
                                    tttcAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCgcgcgc                                                                                                         >  1:1459578/4‑148 (MQ=255)
                                    ttGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCgcgcgc                                                                                                         >  2:1387910/1‑148 (MQ=255)
                                      gCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTg                                                                                                       <  1:2068630/148‑1 (MQ=255)
                                        aGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCa                                                                                                     <  1:2778834/148‑1 (MQ=255)
                                         ggCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCAt                                                                                                    <  1:1596287/148‑1 (MQ=255)
                                          gCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCAtg                                                                                                   <  2:271551/148‑1 (MQ=255)
                                           ccAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCAtgt                                                                                                  <  2:1341911/148‑1 (MQ=255)
                                             aaCTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTgg                                                                                                <  1:193658/148‑1 (MQ=255)
                                               cTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGt                                                                                              <  1:3450957/148‑1 (MQ=255)
                                                   ttCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAg                                                                                          >  1:3091019/1‑148 (MQ=255)
                                                         gTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCaa                                                                                    >  2:794845/1‑148 (MQ=255)
                                                         gTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCaa                                                                                    <  2:443934/148‑1 (MQ=255)
                                                            cTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAata                                                                                 >  2:3146324/1‑148 (MQ=255)
                                                            cTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAata                                                                                 >  2:2115990/1‑148 (MQ=255)
                                                              gCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtg                                                                               <  1:2619388/148‑1 (MQ=255)
                                                              gCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtg                                                                               <  1:1866700/148‑1 (MQ=255)
                                                                 tGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtggtg                                                                            <  2:881337/148‑1 (MQ=255)
                                                                  gACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtggtgg                                                                           <  1:314317/148‑1 (MQ=255)
                                                                  gACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtggtgg                                                                           <  1:1365040/148‑1 (MQ=255)
                                                                  gACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtggtgg                                                                           >  1:1578847/1‑148 (MQ=255)
                                                                   actttCTCTTAGTTTACGGTCAAATCCATAACACCTTTCGTTTTCCAAACCTGAATCTCGACAACTCCGTCCACATAACTACCCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCTCTGCATGTTGGTGAAGCGCCAAATTTtgtggt                                                                          <  2:3003208/145‑1 (MQ=255)
                                                                   aCATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATggtggt                                                                          >  2:623935/1‑148 (MQ=255)
                                                                    cATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTc                                                                         >  1:842951/1‑148 (MQ=255)
                                                                    cATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTc                                                                         >  2:1836980/1‑148 (MQ=255)
                                                                         tctTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTctgctg                                                                    <  2:2590489/148‑1 (MQ=255)
                                                                          ctTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTgg                                                                   >  1:3072666/1‑148 (MQ=255)
                                                                            tCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTgggg                                                                 >  2:2170750/1‑148 (MQ=255)
                                                                              gTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCg                                                               <  2:401120/148‑1 (MQ=255)
                                                                                     ggCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTc                                                        >  1:2153988/1‑148 (MQ=255)
                                                                                       cAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCaa                                                      <  1:2184660/148‑1 (MQ=255)
                                                                                          aTCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTa                                                   <  2:70059/148‑1 (MQ=255)
                                                                                          aTCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTa                                                   >  1:1501758/1‑148 (MQ=255)
                                                                                          aTCCAGAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTa                                                   <  2:68520/148‑1 (MQ=255)
                                                                                           tCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTaa                                                  <  2:2060089/148‑1 (MQ=255)
                                                                                           tCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTaa                                                  >  2:1901475/1‑148 (MQ=255)
                                                                                             cATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACg                                                <  1:687319/148‑1 (MQ=255)
                                                                                                aaCACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaa                                             >  1:1004210/1‑148 (MQ=255)
                                                                                                aaCACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaa                                             <  2:1883449/148‑1 (MQ=255)
                                                                                                aaCACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaa                                             <  2:556698/148‑1 (MQ=255)
                                                                                                aaCACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaa                                             >  2:2352292/1‑148 (MQ=255)
                                                                                                 acacCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaaa                                            >  1:2099608/1‑148 (MQ=255)
                                                                                                   acCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCTCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACg                                          <  2:3345075/148‑1 (MQ=255)
                                                                                                    ccGTTCGTTTTCCCAACCTGAATTTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGa                                         <  2:3050376/148‑1 (MQ=255)
                                                                                                    ccGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGa                                         >  1:1814643/1‑148 (MQ=255)
                                                                                                    ccGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGa                                         >  2:1063156/1‑148 (MQ=255)
                                                                                                       ttcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTa                                      >  2:1668296/1‑148 (MQ=255)
                                                                                                        tcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAt                                     <  1:2996707/148‑1 (MQ=255)
                                                                                                        tcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAt                                     <  2:3232514/148‑1 (MQ=255)
                                                                                                        tcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAt                                     >  2:3399660/1‑148 (MQ=255)
                                                                                                         cgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAtt                                    <  2:2892206/148‑1 (MQ=255)
                                                                                                           ttttCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTg                                  >  2:1080477/1‑148 (MQ=255)
                                                                                                                 cAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGcc                            <  2:1373541/148‑1 (MQ=255)
                                                                                                                  aaCCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTAAGccc                           >  1:3264701/1‑148 (MQ=255)
                                                                                                                   aCCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcg                          >  2:727495/1‑148 (MQ=255)
                                                                                                                   aCCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcg                          >  1:3222195/1‑148 (MQ=255)
                                                                                                                    ccTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgt                         <  2:855726/148‑1 (MQ=255)
                                                                                                                      tGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcg                       <  1:2513565/148‑1 (MQ=255)
                                                                                                                       gAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcgc                      >  2:832417/1‑148 (MQ=255)
                                                                                                                        aaTCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcgcg                     >  2:3321501/1‑148 (MQ=255)
                                                                                                                          tctcGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcgcgtc                   <  1:1887020/148‑1 (MQ=255)
                                                                                                                           ctcGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcgcgtcg                  >  1:3034388/1‑148 (MQ=255)
                                                                                                                               aCAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGaaa              >  2:517176/1‑148 (MQ=255)
                                                                                                                               aCAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGaaa              <  1:336360/148‑1 (MQ=255)
                                                                                                                                cAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAAc             <  1:48167/148‑1 (MQ=255)
                                                                                                                                 aaCTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAAcc            >  2:2277103/1‑148 (MQ=255)
                                                                                                                                   cTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAg          <  1:1066602/148‑1 (MQ=255)
                                                                                                                                   cTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAg          <  2:1112378/148‑1 (MQ=255)
                                                                                                                                    tCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGc         <  2:3335228/148‑1 (MQ=255)
                                                                                                                                    tCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGc         <  1:117099/148‑1 (MQ=255)
                                                                                                                                      cGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTg       <  2:273266/148‑1 (MQ=255)
                                                                                                                                        tCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTggg     >  1:3019922/1‑148 (MQ=255)
                                                                                                                                         ccACATCACAACCCTCCTCATCTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATCTGGTGGTCAGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGc    <  2:2088110/148‑1 (MQ=255)
                                                                                                                                         ccACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGc    >  2:3560316/1‑148 (MQ=255)
                                                                                                                                         ccACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGc    <  2:1501923/148‑1 (MQ=255)
                                                                                                                                           acaTAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCt  >  2:975220/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCT  >  CP009273/2919907‑2920193

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: