Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,999,632 G→T 100% intergenic (+73/‑5) fliT → / → amyA putative flagellar synthesis and assembly chaperone/cytoplasmic alpha‑amylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,999,6320GT98.6% 697.0 / ‑0.8 215intergenic (+73/‑5)fliT/amyAputative flagellar synthesis and assembly chaperone/cytoplasmic alpha‑amylase
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (100/112);  total (103/112)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 218 total aligned reads displayed.

GCAAAAATCGGTGTTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCGAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGCCGCCCGCCTATAAAGG  >  CP009273/1999487‑1999776
                                                                                                                                                 |                                                                                                                                                
      aTCGGTGTTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGc                                                                                                                                          >  1:889292/1‑148 (MQ=255)
       tCGGTGTTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCg                                                                                                                                         <  1:2702206/148‑1 (MQ=255)
       tCGGTGTTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCg                                                                                                                                         >  1:2560917/1‑148 (MQ=255)
        cGGTGTTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGt                                                                                                                                        >  1:1054003/1‑148 (MQ=255)
          gtgtTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTaa                                                                                                                                      <  1:2061792/148‑1 (MQ=255)
          gtgtTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTaa                                                                                                                                      <  1:2420996/148‑1 (MQ=255)
            gtTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATc                                                                                                                                    >  2:1637097/1‑148 (MQ=255)
            gtTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATc                                                                                                                                    >  2:3375499/1‑148 (MQ=255)
             ttAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAAAAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATcc                                                                                                                                   >  1:313540/1‑148 (MQ=255)
                aGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCAc                                                                                                                                >  1:3472991/1‑148 (MQ=255)
                 gTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACg                                                                                                                               <  1:1878254/148‑1 (MQ=255)
                  tGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGc                                                                                                                              <  1:2834621/148‑1 (MQ=255)
                  tGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGc                                                                                                                              >  1:1675633/1‑148 (MQ=255)
                       aTGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTa                                                                                                                         >  1:2043219/1‑148 (MQ=255)
                        tGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTATAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTAc                                                                                                                        >  2:2651651/1‑148 (MQ=255)
                         ggCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCTCGCTGTTAct                                                                                                                       <  2:36847/148‑2 (MQ=255)
                         ggCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACa                                                                                                                       >  2:803297/1‑148 (MQ=255)
                         ggCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACa                                                                                                                       <  2:3350144/148‑1 (MQ=255)
                         ggCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATACCACCCTGTTACa                                                                                                                       >  1:311803/1‑148 (MQ=255)
                          gCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACaa                                                                                                                      <  1:3301066/148‑1 (MQ=255)
                            gATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATg                                                                                                                    <  2:3340765/148‑1 (MQ=255)
                            gATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATg                                                                                                                    >  1:1391865/1‑148 (MQ=255)
                              tCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGtt                                                                                                                  <  2:2466903/148‑1 (MQ=255)
                               cAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGttt                                                                                                                 <  2:2552414/148‑1 (MQ=255)
                                aGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGtttt                                                                                                                <  1:2349538/148‑1 (MQ=255)
                                aGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGtttt                                                                                                                <  1:585528/148‑1 (MQ=255)
                                 gggCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCGAAGAATGCGTAATCCCACGCTGTTACAATGtttta                                                                                                               >  1:986600/1‑147 (MQ=255)
                                  ggcggcTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCa                                                                                                              >  2:1000315/1‑148 (MQ=255)
                                  ggcggcTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCa                                                                                                              >  1:3154444/1‑148 (MQ=255)
                                  ggcggcTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCa                                                                                                              >  2:32556/1‑148 (MQ=255)
                                     ggcTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTg                                                                                                           <  1:2162341/148‑1 (MQ=255)
                                       cTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGt                                                                                                         <  1:1455425/148‑1 (MQ=255)
                                        tttGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTa                                                                                                        <  1:233018/148‑1 (MQ=255)
                                         ttGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTAt                                                                                                       <  2:1939930/148‑1 (MQ=255)
                                          tgtgCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTAtt                                                                                                      <  2:1202994/148‑1 (MQ=255)
                                                 gCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGa                                                                                               <  1:748368/148‑1 (MQ=255)
                                                 gCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGa                                                                                               <  1:2521762/148‑1 (MQ=255)
                                                 gCTCCGCAGGAGAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGa                                                                                               <  2:2616260/148‑1 (MQ=255)
                                                  cTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGaa                                                                                              <  1:1600486/148‑1 (MQ=255)
                                                   tCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAg                                                                                             <  2:2136449/148‑1 (MQ=255)
                                                         ggATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTa                                                                                       <  1:519122/148‑1 (MQ=255)
                                                          gATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTaa                                                                                      >  2:386381/1‑148 (MQ=255)
                                                           aTAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAg                                                                                     <  1:3453977/148‑1 (MQ=255)
                                                             aaCCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGct                                                                                   >  2:880903/1‑148 (MQ=255)
                                                             aaCCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGct                                                                                   >  2:3208222/1‑148 (MQ=255)
                                                              aCCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGctc                                                                                  >  1:172276/1‑148 (MQ=255)
                                                              aCCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGctc                                                                                  <  1:454855/148‑1 (MQ=255)
                                                              aCCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGctc                                                                                  >  1:2143782/1‑148 (MQ=255)
                                                                 tctTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTgg                                                                               >  1:1690677/1‑148 (MQ=255)
                                                                 tctTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTgg                                                                               >  1:1567585/1‑148 (MQ=255)
                                                                  ctTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATTCGTAATCACACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGc                                                                              >  1:2547930/1‑148 (MQ=255)
                                                                   ttttGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGcc                                                                             >  1:1325929/1‑148 (MQ=255)
                                                                   ttttGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGcc                                                                             <  2:515247/148‑1 (MQ=255)
                                                                        aaTCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAAc                                                                        >  1:1226692/1‑148 (MQ=255)
                                                                           cTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTgg                                                                     <  1:394293/148‑1 (MQ=255)
                                                                           cTGAATGAGACGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGCCCTGAACTgg                                                                     >  2:722263/1‑148 (MQ=255)
                                                                             gAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGcc                                                                   >  1:543608/1‑148 (MQ=255)
                                                                                  aGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAgcg                                                              >  2:1705010/1‑148 (MQ=255)
                                                                                  aGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAgcg                                                              >  2:1469349/1‑148 (MQ=255)
                                                                                   gTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAgcgc                                                             >  2:1020789/1‑148 (MQ=255)
                                                                                   gTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAgcgc                                                             >  2:1166266/1‑148 (MQ=255)
                                                                                    tCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAAGGTTCCCCTTGACTGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAgcgcg                                                            <  2:1999172/148‑1 (MQ=255)
                                                                                     cGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAgcgcgc                                                           <  2:634006/148‑1 (MQ=255)
                                                                                       aTGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCg                                                         >  2:2877674/1‑148 (MQ=255)
                                                                                          gCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACg                                                      >  2:1264304/1‑148 (MQ=255)
                                                                                          gCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACg                                                      >  2:2480371/1‑148 (MQ=255)
                                                                                             cgcgAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGtt                                                   <  2:1220618/148‑1 (MQ=255)
                                                                                              gcgAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGttt                                                  <  2:1338082/148‑1 (MQ=255)
                                                                                              gcgAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGttt                                                  <  1:1020780/148‑1 (MQ=255)
                                                                                                gAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTa                                                <  2:377004/148‑1 (MQ=255)
                                                                                                 aataatCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTaa                                               <  2:2366693/148‑1 (MQ=255)
                                                                                                   taatCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATg                                             >  1:1405015/1‑148 (MQ=255)
                                                                                                   taatCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATg                                             <  2:1054003/148‑1 (MQ=255)
                                                                                                     atCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGat                                           >  2:1602961/1‑148 (MQ=255)
                                                                                                       ccGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGatat                                         >  1:2096668/1‑148 (MQ=255)
                                                                                                          aTTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTgg                                      <  2:2984932/148‑1 (MQ=255)
                                                                                                           ttACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGt                                     <  2:313540/148‑1 (MQ=255)
                                                                                                           ttACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGt                                     <  2:2790887/148‑1 (MQ=255)
                                                                                                             aCGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTAt                                   <  1:1093830/148‑1 (MQ=255)
                                                                                                                gCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCaa                                >  1:2918863/1‑148 (MQ=255)
                                                                                                                gCTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCaa                                >  1:3039896/1‑148 (MQ=255)
                                                                                                                 cTACTCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAat                               <  1:355403/148‑1 (MQ=255)
                                                                                                                 cTACGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAat                               <  1:2125658/148‑1 (MQ=255)
                                                                                                                    cGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATg                            <  1:1653081/148‑1 (MQ=255)
                                                                                                                    cGCTTCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATg                            <  1:2239588/148‑1 (MQ=255)
                                                                                                                       ttCTAATGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTc                         >  2:2918346/1‑148 (MQ=255)
                                                                                                                           aaTGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGt                     >  1:2607548/1‑148 (MQ=255)
                                                                                                                            aTGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATTGTCTGGtt                    >  1:3317619/1‑148 (MQ=255)
                                                                                                                            aTGTTCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGtt                    >  2:840141/1‑148 (MQ=255)
                                                                                                                                tCCCCTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTgccg                >  1:2448588/1‑148 (MQ=255)
                                                                                                                                 ccccTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTgccgc               >  2:2310324/1‑148 (MQ=255)
                                                                                                                                  cccTTGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTgccgcc              <  2:1889659/148‑1 (MQ=255)
                                                                                                                                     ttGAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGccgcccgc           <  1:803297/148‑1 (MQ=255)
                                                                                                                                       gAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGCCGCCCGCCt         <  1:1145719/148‑1 (MQ=255)
                                                                                                                                       gAATGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGCCGCCCGCCt         <  1:2651508/148‑1 (MQ=255)
                                                                                                                                        aaTGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGCCGCCCGCCta        >  2:2293572/1‑148 (MQ=255)
                                                                                                                                        aaTGGAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGCCGCCCGCCta        >  2:69816/1‑148 (MQ=255)
                                                                                                                                           ggAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGCCGCCCGCCTATaa     >  1:196389/1‑148 (MQ=255)
                                                                                                                                            gAGTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGCCGCCCGCCTATaaa    >  2:2555809/1‑148 (MQ=255)
                                                                                                                                              gTCTAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGCCGCCCGCCTATAAAgg  >  2:784339/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAAAAATCGGTGTTAAGTGCCTATGGCGATCAGGGCGGCTTTGTGCTGGCTCCGCAGGATAACCTCTTTTGAATCTGAATGAGTCGATGGCTCGCGAATAATCCGATTACGGCTACGCTTCTAATGTTCCCCTTGAATGGAGTCGAAGAATGCGTAATCCCACGCTGTTACAATGTTTTCACTGGTATTACCCGGAAGGCGGTAAGCTCTGGCCTGAACTGGCCGAGCGCGCCGACGGTTTTAATGATATTGGTATCAATATGGTCTGGTTGCCGCCCGCCTATAAAGG  >  CP009273/1999487‑1999776

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: