Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,189,569 C→A 38.0% D182Y (GAT→TAT)  hldE ← fused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,189,5690CA38.0% 270.8 / 168.9 193D182Y (GAT→TAT) hldEfused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base C (68/51);  new base A (49/25);  total (117/76)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

TCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAG  >  CP009273/3189422‑3189716
                                                                                                                                                   |                                                                                                                                                   
tCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATc                                                                                                                                                     <  2:3475683/148‑1 (MQ=255)
   cTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGt                                                                                                                                                  >  2:2970102/1‑148 (MQ=255)
    tAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTa                                                                                                                                                 >  1:486066/1‑148 (MQ=255)
    tAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTa                                                                                                                                                 >  1:3298760/1‑148 (MQ=255)
       cAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCt                                                                                                                                              <  2:2078550/148‑1 (MQ=255)
       cAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCt                                                                                                                                              >  1:2075885/1‑148 (MQ=255)
       cAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGAGCTCAAAATCGGTACCt                                                                                                                                              <  1:2067545/148‑1 (MQ=255)
        agagCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCtt                                                                                                                                             <  2:203628/148‑1 (MQ=255)
         gagCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCttt                                                                                                                                            <  2:2756610/148‑1 (MQ=255)
          agCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCtttt                                                                                                                                           >  2:2791899/1‑148 (MQ=255)
          agCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCtttt                                                                                                                                           >  2:3385547/1‑148 (MQ=255)
          agCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCtttt                                                                                                                                           >  1:1982006/1‑148 (MQ=255)
          agCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCtttt                                                                                                                                           >  2:3607549/1‑148 (MQ=255)
          agCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCtttt                                                                                                                                           >  1:1847095/1‑148 (MQ=255)
          agCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCtttt                                                                                                                                           <  1:1503869/148‑1 (MQ=255)
          agCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCtttt                                                                                                                                           >  1:2820082/1‑148 (MQ=255)
                gagTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatca                                                                                                                                     >  2:1592069/1‑148 (MQ=255)
                   ttCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatcaatc                                                                                                                                  <  1:3624208/148‑1 (MQ=255)
                   ttCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatcaatc                                                                                                                                  <  1:2296709/148‑1 (MQ=255)
                    tCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatcaatca                                                                                                                                 >  1:3166519/1‑148 (MQ=255)
                     cGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAg                                                                                                                                >  1:509129/1‑148 (MQ=255)
                        aaTCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCAc                                                                                                                             <  2:239495/148‑1 (MQ=255)
                                 tCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACAcc                                                                                                                    >  2:2198738/1‑148 (MQ=255)
                                 tCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACAcc                                                                                                                    <  1:3418144/148‑1 (MQ=255)
                                 tCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACAcc                                                                                                                    <  1:3417570/148‑1 (MQ=255)
                                   aGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCg                                                                                                                  >  2:792040/1‑148 (MQ=255)
                                   aGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCg                                                                                                                  <  1:2828303/148‑1 (MQ=255)
                                     tttCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCt                                                                                                                <  2:1769052/148‑1 (MQ=255)
                                     tttCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCt                                                                                                                >  1:1705284/1‑148 (MQ=255)
                                     tttCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCt                                                                                                                >  1:1941388/1‑148 (MQ=255)
                                          tGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTAcg                                                                                                           >  2:3653379/1‑148 (MQ=255)
                                             cgcgCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACgcgc                                                                                                        <  2:2246486/148‑1 (MQ=255)
                                               cgcTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCa                                                                                                      <  2:956243/148‑1 (MQ=255)
                                                gcTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAg                                                                                                     <  1:1491222/148‑1 (MQ=255)
                                                  tCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGtt                                                                                                   >  1:500356/1‑148 (MQ=255)
                                                   caacaaTCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTg                                                                                                  <  1:1958072/148‑1 (MQ=255)
                                                    aacaaTCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTgg                                                                                                 >  1:3607110/1‑148 (MQ=255)
                                                     acaaTCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGa                                                                                                >  2:681038/1‑148 (MQ=255)
                                                       aaTCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGatc                                                                                              >  1:1345254/1‑148 (MQ=255)
                                                        aTCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCCAAATCGGTACCTTTTGGATCAATCAGCACCGGACCACCCGCTATACGCGCCAGATGGatca                                                                                             >  1:3205536/1‑148 (MQ=255)
                                                        aTCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGatca                                                                                             >  1:1381639/1‑148 (MQ=255)
                                                         tctcTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGatcat                                                                                            >  2:100079/1‑148 (MQ=255)
                                                          ctctTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGatcatc                                                                                           <  1:3491707/148‑1 (MQ=255)
                                                          ctctTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGatcatc                                                                                           <  2:264475/148‑1 (MQ=255)
                                                             ttcttcGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATctgc                                                                                        <  2:2655388/148‑1 (MQ=255)
                                                              tcttcGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATctgct                                                                                       >  1:825036/1‑148 (MQ=255)
                                                              tcttcGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATctgct                                                                                       >  1:99202/1‑148 (MQ=255)
                                                               cttcGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATctgctg                                                                                      <  2:1513094/148‑1 (MQ=255)
                                                                ttcGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGt                                                                                     <  1:3092885/148‑1 (MQ=255)
                                                                   ggTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACg                                                                                  <  2:1826489/148‑1 (MQ=255)
                                                                   ggTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACg                                                                                  <  2:137155/148‑1 (MQ=255)
                                                                    gTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGc                                                                                 >  2:1583117/1‑148 (MQ=255)
                                                                    gTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGc                                                                                 <  1:1671095/148‑1 (MQ=255)
                                                                      cTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCtt                                                                               >  2:3482427/1‑148 (MQ=255)
                                                                      cTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCtt                                                                               >  2:418385/1‑148 (MQ=255)
                                                                            tttACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAgc                                                                         >  2:804863/1‑148 (MQ=255)
                                                                               aCCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCa                                                                      >  1:3375125/1‑148 (MQ=255)
                                                                                 cGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCAcc                                                                    <  1:3303658/148‑1 (MQ=255)
                                                                                 cGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCAcc                                                                    <  1:2282714/148‑1 (MQ=255)
                                                                                   acaacaGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCtt                                                                  >  1:874279/1‑148 (MQ=255)
                                                                                     aacaGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTg                                                                >  2:2685054/1‑148 (MQ=255)
                                                                                       caGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGc                                                              >  2:3621070/1‑148 (MQ=255)
                                                                                       caGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGc                                                              >  2:2101418/1‑148 (MQ=255)
                                                                                        aGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCTGGGGGGCCCCCGCGGGCCGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCg                                                             >  2:1743893/1‑148 (MQ=255)
                                                                                        aGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCGGGCCGGGGGGGGTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCg                                                             >  2:1744401/1‑148 (MQ=255)
                                                                                          cTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTa                                                           <  1:1355802/148‑1 (MQ=255)
                                                                                           ttCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAg                                                          >  2:3095050/1‑148 (MQ=255)
                                                                                            tCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGt                                                         >  1:3337788/1‑148 (MQ=255)
                                                                                                aTTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGa                                                     >  2:1379861/1‑148 (MQ=255)
                                                                                                 ttCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGaa                                                    >  2:54232/1‑148 (MQ=255)
                                                                                                 ttCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGaa                                                    <  2:3166519/148‑1 (MQ=255)
                                                                                                 ttCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGaa                                                    >  2:1051725/1‑148 (MQ=255)
                                                                                                 ttCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGaa                                                    >  1:2950611/1‑148 (MQ=255)
                                                                                                 ttCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGaa                                                    <  1:2410629/148‑1 (MQ=255)
                                                                                                  tCCGAGAGATTCGGCGTCAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGaaa                                                   >  1:3792878/1‑148 (MQ=255)
                                                                                                    cGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGc                                                 <  2:2114275/148‑1 (MQ=255)
                                                                                                    cGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGc                                                 >  1:2020330/1‑148 (MQ=255)
                                                                                                      agagaTTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCAc                                               >  2:3130794/1‑148 (MQ=255)
                                                                                                          aTTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAgc                                           >  1:611707/1‑148 (MQ=255)
                                                                                                          aTTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAgc                                           >  1:103414/1‑148 (MQ=255)
                                                                                                          aTTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAgc                                           >  1:103717/1‑148 (MQ=255)
                                                                                                          aTTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAgc                                           >  1:1980553/1‑148 (MQ=255)
                                                                                                                 gTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCaa                                    >  2:2091935/1‑148 (MQ=255)
                                                                                                                  ttAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAAt                                   <  2:2855926/148‑1 (MQ=255)
                                                                                                                   tAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATc                                  >  1:2892125/1‑148 (MQ=255)
                                                                                                                   tAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATc                                  >  2:2138706/1‑148 (MQ=255)
                                                                                                                       aGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAAc                              <  2:2874116/148‑1 (MQ=255)
                                                                                                                         cGTAGCGCCGCGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTc                            <  1:2791143/148‑1 (MQ=255)
                                                                                                                              cgccgcGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAgcgc                       <  1:2375764/148‑1 (MQ=255)
                                                                                                                              cgccgcGGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAgcgc                       >  2:613160/1‑148 (MQ=255)
                                                                                                                                   cgGTAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGAt                  <  1:3663431/148‑1 (MQ=255)
                                                                                                                                     gTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTa                <  1:1291603/148‑1 (MQ=255)
                                                                                                                                      tAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTaa               <  2:3298760/148‑1 (MQ=255)
                                                                                                                                      tAGCGCTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTaa               >  1:742861/1‑148 (MQ=255)
                                                                                                                                        gcgcTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATc             <  2:1190117/148‑1 (MQ=255)
                                                                                                                                          gcTCAAAATAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCg           <  1:3318659/148‑1 (MQ=255)
                                                                                                                                               aaaTAGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGt      >  2:731475/1‑148 (MQ=255)
                                                                                                                                                 aTCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGc    >  2:2494736/1‑148 (MQ=255)
                                                                                                                                                   cGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAg  <  2:1503417/148‑1 (MQ=255)
                                                                                                                                                   cGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAg  <  2:873933/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAG  >  CP009273/3189422‑3189716

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: