Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 857,186 G→T 12.0% L293I (CTC→ATC)  ybiW ← putative pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273857,1860GT12.0% 585.9 / 39.0 175L293I (CTC→ATC) ybiWputative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (69/85);  new base T (10/11);  total (79/96)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 189 total aligned reads displayed.

GACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAG  >  CP009273/857039‑857333
                                                                                                                                                   |                                                                                                                                                   
gACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTagag                                                                                                                                                     >  1:83357/1‑148 (MQ=255)
gACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTagag                                                                                                                                                     <  1:1219105/148‑1 (MQ=255)
gACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTagag                                                                                                                                                     >  1:1159368/1‑148 (MQ=255)
 aCTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAgaga                                                                                                                                                    >  2:236944/1‑148 (MQ=255)
 aCTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGGCGCGGCGATAGTACGGGTAgaga                                                                                                                                                    >  2:236989/1‑148 (MQ=255)
  cTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTATGGGTAGAGAt                                                                                                                                                   >  1:671759/1‑148 (MQ=255)
  cTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGAt                                                                                                                                                   <  1:3304724/148‑1 (MQ=255)
      ccGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTg                                                                                                                                               <  2:333564/148‑1 (MQ=255)
       cGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTgg                                                                                                                                              <  2:2302352/148‑1 (MQ=255)
         cAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTc                                                                                                                                            <  1:1909346/148‑1 (MQ=255)
          agagGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTcc                                                                                                                                           >  2:178054/1‑148 (MQ=255)
             ggCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATa                                                                                                                                        <  1:1272180/148‑1 (MQ=255)
             ggCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATa                                                                                                                                        >  1:236881/1‑148 (MQ=255)
              gCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATAc                                                                                                                                       >  2:2981477/1‑148 (MQ=255)
               cTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACg                                                                                                                                      >  2:1394116/1‑148 (MQ=255)
                 tttGAGTGTGAGCCGGAGCTGATCTTTTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGAAGCATCTCGATGGCGTGTTCGCGATCCCGCGTCTGGTTGAGTTCAACGTCGCGGCGATATTACGGGTAGAGATACTGGTCCATACGAc                                                                                                                                    <  1:1124690/148‑1 (MQ=255)
                 tttGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGAc                                                                                                                                    <  1:1124763/148‑1 (MQ=255)
                 tttGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACATCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGAc                                                                                                                                    >  1:367727/1‑148 (MQ=255)
                 tttGAGTGTGAGCCGGAGCGGATCTTGTTCACGTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGAc                                                                                                                                    <  1:2031568/148‑1 (MQ=255)
                  ttGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGAcc                                                                                                                                   <  1:545450/148‑1 (MQ=255)
                    gAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCaa                                                                                                                                 <  2:1279237/148‑1 (MQ=255)
                    gAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCaa                                                                                                                                 <  2:1279271/148‑1 (MQ=255)
                    gAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCaa                                                                                                                                 <  1:2028708/148‑1 (MQ=255)
                    gAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCaa                                                                                                                                 <  2:1398732/148‑1 (MQ=255)
                       tgtgAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACg                                                                                                                              >  2:948245/1‑148 (MQ=255)
                          gAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATa                                                                                                                           <  1:1603118/148‑1 (MQ=255)
                               ggAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAg                                                                                                                      >  1:2998706/1‑148 (MQ=255)
                               ggAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAg                                                                                                                      <  2:410975/148‑1 (MQ=255)
                                 aGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTg                                                                                                                    <  1:2537017/148‑1 (MQ=255)
                                      aTCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGt                                                                                                               >  2:492045/1‑148 (MQ=255)
                                       tCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGtt                                                                                                              <  1:2671755/148‑1 (MQ=255)
                                        cTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACAGTTa                                                                                                             >  1:3209773/1‑148 (MQ=255)
                                         ttGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTNAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAg                                                                                                            <  1:1682127/148‑1 (MQ=255)
                                          tGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGa                                                                                                           <  2:792788/148‑1 (MQ=255)
                                            ttCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGAtt                                                                                                         >  1:1903991/1‑148 (MQ=255)
                                                 ttCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATc                                                                                                    <  1:2669839/148‑1 (MQ=255)
                                                 ttCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATc                                                                                                    <  1:1016342/148‑1 (MQ=255)
                                                      gcagTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTg                                                                                                                    <  2:3240216/127‑1 (MQ=255)
                                                      gcagTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTg                                                                                                                    >  1:3240188/1‑127 (MQ=255)
                                                      gcagTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTg                                                                                                                    >  1:3240216/1‑127 (MQ=255)
                                                       cagTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCaaa                                                                                              <  2:475379/148‑1 (MQ=255)
                                                       cagTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCaaa                                                                                              <  2:2033418/148‑1 (MQ=255)
                                                       cagTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCaaa                                                                                              <  1:2375170/148‑1 (MQ=255)
                                                       cagTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCaaa                                                                                              >  2:1196408/1‑148 (MQ=255)
                                                        agTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCaaaa                                                                                             <  1:1170525/148‑1 (MQ=255)
                                                         gTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAAt                                                                                            <  1:2670610/148‑1 (MQ=255)
                                                          tttCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATc                                                                                           >  2:2366974/1‑148 (MQ=255)
                                                          tttCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATc                                                                                           <  2:1328196/148‑1 (MQ=255)
                                                              agccagcAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACt                                                                                       >  1:276947/1‑148 (MQ=255)
                                                                ccagcAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTgg                                                                                     <  1:2849142/148‑1 (MQ=255)
                                                                ccagcAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTgg                                                                                     <  1:2816220/148‑1 (MQ=255)
                                                                    cagcTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGa                                                                                 >  2:2107441/1‑148 (MQ=255)
                                                                      gcTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAg                                                                               >  1:1035742/1‑148 (MQ=255)
                                                                       cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGt                                                                              >  1:2747769/1‑148 (MQ=255)
                                                                       cTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGt                                                                              <  1:2865114/148‑1 (MQ=255)
                                                                         aTGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTaa                                                                            <  2:1436319/148‑1 (MQ=255)
                                                                          tGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAac                                                                           >  2:2495425/1‑148 (MQ=255)
                                                                          tGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAac                                                                           <  2:3248934/148‑1 (MQ=255)
                                                                            cagcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAacac                                                                         <  1:295509/148‑1 (MQ=255)
                                                                              gcaTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAg                                                                       >  1:278553/1‑148 (MQ=255)
                                                                               caTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGTTTCGATCTGCACAATCAACTGGATGAAGTAACACAGt                                                                      >  2:1923919/1‑148 (MQ=255)
                                                                               caTCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGt                                                                      <  2:1139356/148‑1 (MQ=255)
                                                                                    cGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAg                                                                 <  2:2942510/148‑1 (MQ=255)
                                                                                    cGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAg                                                                 >  2:817120/1‑148 (MQ=255)
                                                                                      aTGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAgcg                                                               >  2:2510094/1‑148 (MQ=255)
                                                                                      aTGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAgcg                                                               >  2:1707210/1‑148 (MQ=255)
                                                                                       tGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAgcgc                                                              >  1:1797279/1‑148 (MQ=255)
                                                                                          cGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTg                                                           >  1:1955923/1‑148 (MQ=255)
                                                                                              ttCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAg                                                       >  1:1489752/1‑148 (MQ=255)
                                                                                              ttCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAg                                                       >  1:1400639/1‑148 (MQ=255)
                                                                                                cgcgATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGaa                                                     <  2:1022213/148‑1 (MQ=255)
                                                                                                     tCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCt                                                >  2:633769/1‑148 (MQ=255)
                                                                                                      ccAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTg                                               >  1:490594/1‑148 (MQ=255)
                                                                                                        aGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTgcg                                             >  1:728877/1‑148 (MQ=255)
                                                                                                          cGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTgcggc                                           >  1:2137187/1‑148 (MQ=255)
                                                                                                            tCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGcggcgg                                         <  2:1141087/148‑1 (MQ=255)
                                                                                                             cTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCggcggc                                        <  2:2651329/148‑1 (MQ=255)
                                                                                                                gTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCtgg                                     <  2:2459615/148‑1 (MQ=255)
                                                                                                                gTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCtgg                                     <  1:10663/148‑1 (MQ=255)
                                                                                                                gTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCtgg                                     <  1:2505909/148‑1 (MQ=255)
                                                                                                                gTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCtgg                                     <  1:355843/148‑1 (MQ=255)
                                                                                                                 ttGAGTTCCACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCtggt                                    <  1:3093830/148‑1 (MQ=255)
                                                                                                                 ttGAGTTCAAGGTCGCGGAGATAGGACGGGTAGAGATACTGGTCCATACGACAAAACGATACAGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTCACACAGTTGCAGAGCCTGCCAGAAACTCTGCGGCTGCttgt                                    >  2:1626292/1‑148 (MQ=255)
                                                                                                                  tGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCtggtg                                   >  1:143308/1‑148 (MQ=255)
                                                                                                                   gAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCtggtgg                                  >  2:176875/1‑148 (MQ=255)
                                                                                                                         aaCGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATa                            <  1:1394116/148‑1 (MQ=255)
                                                                                                                          aCGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATaa                           >  1:1542262/1‑148 (MQ=255)
                                                                                                                           cGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAg                          >  1:1336840/1‑148 (MQ=255)
                                                                                                                           cGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAg                          >  2:1863609/1‑148 (MQ=255)
                                                                                                                                   cgATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGt                  <  1:2412482/148‑1 (MQ=255)
                                                                                                                                   cgATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGt                  >  2:3179880/1‑148 (MQ=255)
                                                                                                                                   cgATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGt                  >  2:1538687/1‑148 (MQ=255)
                                                                                                                                    gATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGtt                 <  2:1144191/148‑1 (MQ=255)
                                                                                                                                       aGTACGGGTAGATATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTc              >  1:85430/1‑148 (MQ=255)
                                                                                                                                       aGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTc              <  1:178054/148‑1 (MQ=255)
                                                                                                                                         tACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTg            >  1:39688/1‑148 (MQ=255)
                                                                                                                                         tACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTg            >  1:464436/1‑148 (MQ=255)
                                                                                                                                           cGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGcc          <  1:1734750/148‑1 (MQ=255)
                                                                                                                                             ggTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCAt        <  1:2409618/148‑1 (MQ=255)
                                                                                                                                                  agaTACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCca   >  2:1852828/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATACGACCAAACGATACTGAGTGACCGTTAGATTCGATCTGCAAAATCAACTGGATGAAGTAACACAGTTGCAGCGCCTGCCAGAAAGTCTGCGGCGGCTGGTGGGCGATAAGATCGCAGTTTTCTGCCATCGCCAG  >  CP009273/857039‑857333

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: