Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,206,749 A→T 5.6% intergenic (+114/‑120) mcrA → / → icdC putative 5‑methylcytosine/5‑hydroxymethylcytosine‑specific restriction nuclease; 5‑methylcytosine DNA binding protein; e14 prophage gene/pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,206,7490AT 5.6% 546.9 / 18.4 160intergenic (+114/‑120)mcrA/icdCputative 5‑methylcytosine/5‑hydroxymethylcytosine‑specific restriction nuclease; 5‑methylcytosine DNA binding protein; e14 prophage gene/pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion
Reads supporting (aligned to +/- strand):  ref base A (93/58);  new base T (5/4);  total (98/62)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 164 total aligned reads displayed.

GCTTTACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGC  >  CP009273/1206602‑1206897
                                                                                                                                                   |                                                                                                                                                    
gCTTTACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAAc                                                                                                                                                    <  1:4135140/150‑1 (MQ=255)
 cTTTACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACt                                                                                                                                                   <  1:2587711/150‑1 (MQ=255)
 cTTTACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACt                                                                                                                                                   <  1:3796635/150‑1 (MQ=255)
    tACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCaa                                                                                                                                                >  1:4444798/1‑150 (MQ=255)
      cGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCaaaa                                                                                                                                              >  1:3604777/1‑150 (MQ=255)
      cGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCaaaa                                                                                                                                              >  1:1976288/1‑150 (MQ=255)
      cGTTAATATAAACCGATTACAGAAATAAAATTATTAATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCaaaa                                                                                                                                              >  1:1976113/1‑150 (MQ=255)
       gTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCaaaaa                                                                                                                                             >  1:1643400/1‑150 (MQ=255)
       gTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCaaaaa                                                                                                                                             >  1:1643318/1‑150 (MQ=255)
        ttAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAac                                                                                                                                            <  2:528989/150‑1 (MQ=255)
           atatAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAAcaca                                                                                                                                         <  2:415771/150‑1 (MQ=255)
            tataAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACAt                                                                                                                                        <  1:3479923/150‑1 (MQ=255)
             ataAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATg                                                                                                                                       <  2:2269103/150‑1 (MQ=255)
               aaaCCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGtt                                                                                                                                     >  2:2781818/1‑150 (MQ=255)
                aaCCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGttt                                                                                                                                    >  2:1177136/1‑150 (MQ=255)
                aaCCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGttt                                                                                                                                    >  2:3340491/1‑150 (MQ=255)
                     aTTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAAt                                                                                                                               <  1:2969302/150‑1 (MQ=255)
                     aTTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAAt                                                                                                                               <  1:2968510/150‑1 (MQ=255)
                          aGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATg                                                                                                                          >  2:3566565/1‑150 (MQ=255)
                                aaaaTTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAATACACATGTTTATAATCAATGAGTtat                                                                                                                    >  2:4589920/1‑150 (MQ=255)
                                 aaaTTATTTATTAAAGTCACATTTAAGTCGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAAAGAGTtata                                                                                                                   >  1:3876657/1‑150 (MQ=255)
                                 aaaTTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTtata                                                                                                                   >  1:186860/1‑150 (MQ=255)
                                 aaaTTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTtata                                                                                                                   >  1:3876668/1‑150 (MQ=255)
                                     tatttattAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGaaa                                                                                                               <  1:1738843/150‑1 (MQ=255)
                                      atttattAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAAt                                                                                                              <  1:1716569/150‑1 (MQ=255)
                                       tttattAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATg                                                                                                             >  1:2084607/1‑150 (MQ=255)
                                        ttattaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGc                                                                                                            >  2:691869/1‑150 (MQ=255)
                                         tattaAAGTCACATTTAAGACGTAATACCCTACAGGTTAAAAATTTTCTCTGATCTTAACTTATGCAAATGTTAACTTCTATTTTTATGCTAAAAATGGTTATAAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCt                                                                                                           >  1:1630598/1‑150 (MQ=255)
                                         tattaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCt                                                                                                           >  1:2448881/1‑150 (MQ=255)
                                         tattaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTATGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCt                                                                                                           >  1:5132520/1‑150 (MQ=255)
                                          attaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTa                                                                                                          <  2:2244741/150‑1 (MQ=255)
                                            taAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAg                                                                                                        >  2:2458272/1‑150 (MQ=255)
                                             aaaGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAgg                                                                                                       <  2:1247787/150‑1 (MQ=255)
                                             aaaGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAgg                                                                                                       >  1:5046037/1‑150 (MQ=255)
                                               aGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGc                                                                                                     <  1:4969709/150‑1 (MQ=255)
                                                gTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATCGAAATGCTAAGGGCt                                                                                                    >  2:3142192/1‑150 (MQ=255)
                                                   acaTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAAt                                                                                                 <  2:2592597/150‑1 (MQ=255)
                                                      tttAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAg                                                                                              <  1:4778440/150‑1 (MQ=255)
                                                      tttAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAg                                                                                              <  2:195552/150‑1 (MQ=255)
                                                      tttAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAg                                                                                              <  1:2318001/150‑1 (MQ=255)
                                                       ttAAGACGTAATACACTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTCTAGACATGCTCAGGGCTCATGAGt                                                                                             >  2:985053/1‑150 (MQ=255)
                                                        tAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTTCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGtt                                                                                            <  1:4894342/150‑1 (MQ=255)
                                                        tAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGtt                                                                                            <  2:4252825/150‑1 (MQ=255)
                                                            aCGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTatat                                                                                        <  2:2942148/150‑1 (MQ=255)
                                                            aCGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTatat                                                                                        <  1:2067657/150‑1 (MQ=255)
                                                            aCGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTatat                                                                                        <  1:4552329/150‑1 (MQ=255)
                                                            aCGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTatat                                                                                        <  1:2067661/150‑1 (MQ=255)
                                                            aCGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTatat                                                                                        <  2:2942125/150‑1 (MQ=255)
                                                             cGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGATATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATg                                                                                       <  1:4172657/150‑1 (MQ=255)
                                                               tAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCa                                                                                     >  2:3199559/1‑150 (MQ=255)
                                                               tAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCa                                                                                     >  2:941305/1‑150 (MQ=255)
                                                                aaTACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCaa                                                                                    <  2:4918980/150‑1 (MQ=255)
                                                                            ggTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAAtt                                                                        >  1:4681062/1‑150 (MQ=255)
                                                                              tAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTAt                                                                      <  1:5021247/150‑1 (MQ=255)
                                                                               aaaaaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATg                                                                     <  1:1026969/150‑1 (MQ=255)
                                                                                aaaaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGt                                                                    >  2:5084077/1‑150 (MQ=255)
                                                                                 aaaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGtt                                                                   <  2:375634/150‑1 (MQ=255)
                                                                                  aaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                  >  2:987806/1‑150 (MQ=255)
                                                                                  aaTTTTCACTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                  >  2:3444463/1‑150 (MQ=255)
                                                                                    ttttCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCt                                                                <  1:241803/150‑1 (MQ=255)
                                                                                     tttCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTa                                                               <  2:1816107/150‑1 (MQ=255)
                                                                                       tctcTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATg                                                             >  1:1033775/1‑150 (MQ=255)
                                                                                          ctGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCa                                                          >  2:3902894/1‑150 (MQ=255)
                                                                                            gATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGa                                                        >  2:3560834/1‑150 (MQ=255)
                                                                                            gATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGa                                                        >  1:2002047/1‑150 (MQ=255)
                                                                                             aTCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGAt                                                       <  2:1933552/150‑1 (MQ=255)
                                                                                                ttAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGt                                                    <  2:5061355/150‑1 (MQ=255)
                                                                                                ttAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGt                                                    <  1:2666623/150‑1 (MQ=255)
                                                                                                 tAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGtt                                                   <  2:4414571/150‑1 (MQ=255)
                                                                                                   aCTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTATATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTAc                                                 <  1:5096657/150‑1 (MQ=255)
                                                                                                       cTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGAtt                                             >  1:1193391/1‑150 (MQ=255)
                                                                                                           aaaTGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGt                                         >  1:1076743/1‑150 (MQ=255)
                                                                                                                  aaCTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaa                                  <  2:4852307/150‑1 (MQ=255)
                                                                                                                  aaCTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaa                                  <  2:4852273/150‑1 (MQ=255)
                                                                                                                    cTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaact                                >  1:1211900/1‑150 (MQ=255)
                                                                                                                     tGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaacta                               <  2:3656266/150‑1 (MQ=255)
                                                                                                                       cTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAt                             <  2:3557648/150‑1 (MQ=255)
                                                                                                                         aTTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgc                           >  2:4593249/1‑150 (MQ=255)
                                                                                                                          tttttATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgct                          <  1:4425605/150‑1 (MQ=255)
                                                                                                                          tttttATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgct                          <  2:743986/150‑1 (MQ=255)
                                                                                                                            tttATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGATAAGGTCTAATGAGTTATATGCAAATAAGTAAAAAAATGTTTATATGTCAGATAGTTACGAATAAGTCATCTAACTAAtgctgc                        >  1:3665744/1‑150 (MQ=255)
                                                                                                                            tttATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgctgc                        >  2:3835936/1‑150 (MQ=255)
                                                                                                                             ttATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTgcg                       >  1:3902989/1‑150 (MQ=255)
                                                                                                                                 gCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCat                   <  2:2656061/150‑1 (MQ=255)
                                                                                                                                   tAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                 >  1:4254602/1‑150 (MQ=255)
                                                                                                                                   tAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                 >  1:4254587/1‑150 (MQ=255)
                                                                                                                                   tAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                 >  2:1556817/1‑150 (MQ=255)
                                                                                                                                   tAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                 >  2:1556880/1‑150 (MQ=255)
                                                                                                                                    aaaaaTGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGATGCGCCAGATg                >  2:929954/1‑150 (MQ=255)
                                                                                                                                      aaaTGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATggg              >  1:1548808/1‑150 (MQ=255)
                                                                                                                                      aaaTGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATggg              >  1:1548831/1‑150 (MQ=255)
                                                                                                                                      aaaTGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATggg              >  1:1375163/1‑150 (MQ=255)
                                                                                                                                        aTGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGtt            <  2:3508569/150‑1 (MQ=255)
                                                                                                                                          ggTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTgg          >  1:2289932/1‑150 (MQ=255)
                                                                                                                                          ggTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTCATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTgg          >  1:2289891/1‑150 (MQ=255)
                                                                                                                                             tATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACt       >  1:4028280/1‑150 (MQ=255)
                                                                                                                                              aTCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTg      <  2:1339655/150‑1 (MQ=255)
                                                                                                                                               tCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGa     <  1:1177136/150‑1 (MQ=255)
                                                                                                                                                  aaaCTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGc  <  1:2781818/150‑1 (MQ=255)
                                                                                                                                                  aaaCTCAAAAACACATGTTTATAATCAATGAGTTATAAAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGc  >  1:1059419/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTTTACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGC  >  CP009273/1206602‑1206897

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: