Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,206,749 A→T 29.4% 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,7490AT29.4% 100.5 / 58.1 68intergenic (+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 (24/24);  new base T (5/15);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 6.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

TTTACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAG  >  CP009273/1206604‑1206896
                                                                                                                                                 |                                                                                                                                                   
tttACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTc                                                                                                                                                 <  1:1192440/150‑1 (MQ=255)
 ttACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCa                                                                                                                                                >  1:109784/1‑150 (MQ=255)
        aaTATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAacac                                                                                                                                         >  1:71567/1‑150 (MQ=255)
          tataAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACAt                                                                                                                                       <  2:1893785/150‑1 (MQ=255)
          tataAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACAt                                                                                                                                       <  2:4144581/150‑1 (MQ=255)
           ataAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATg                                                                                                                                      >  2:3127199/1‑150 (MQ=255)
           ataAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATg                                                                                                                                      >  2:3127246/1‑150 (MQ=255)
                 cGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTtata                                                                                                                                <  1:2312186/150‑1 (MQ=255)
                 cGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTtata                                                                                                                                <  1:2312187/150‑1 (MQ=255)
                 cGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTtata                                                                                                                                <  1:2312193/150‑1 (MQ=255)
                  gATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATaa                                                                                                                               >  2:3986316/1‑150 (MQ=255)
                  gATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCGAAAAATGGGTATCAAAACTCAAAAAAACGTGTTTCTaa                                                                                                                               >  2:3985684/1‑150 (MQ=255)
                       cAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAAt                                                                                                                          >  1:2478040/1‑150 (MQ=255)
                       cAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAAt                                                                                                                          >  1:2829627/1‑150 (MQ=255)
                        aGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATg                                                                                                                         >  2:2194604/1‑150 (MQ=255)
                            aTAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGtt                                                                                                                     <  1:1326530/150‑1 (MQ=255)
                                            aaGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAggg                                                                                                     <  1:3213202/150‑1 (MQ=255)
                                            aaGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAggg                                                                                                     <  1:2673634/150‑1 (MQ=255)
                                            aaGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAggg                                                                                                     <  1:2672922/150‑1 (MQ=255)
                                                cacaTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTaa                                                                                                 >  1:3675540/1‑150 (MQ=255)
                                                cacaTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTaa                                                                                                 >  1:3675804/1‑150 (MQ=255)
                                                 acaTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAAt                                                                                                <  1:3344557/150‑1 (MQ=255)
                                                 acaTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAAt                                                                                                <  1:200893/150‑1 (MQ=255)
                                                  caTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATg                                                                                               >  1:219982/1‑150 (MQ=255)
                                                              aaTACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCaa                                                                                   <  1:2596711/150‑1 (MQ=255)
                                                                tACCCTCCAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAAt                                                                                 <  1:2612096/150‑1 (MQ=255)
                                                                       cAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTaaa                                                                          <  1:3485407/150‑1 (MQ=255)
                                                                       cAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTaaa                                                                          <  1:3485952/150‑1 (MQ=255)
                                                                        aGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTaaaa                                                                         >  2:2083729/1‑150 (MQ=255)
                                                                             aaaaaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATg                                                                    <  2:71567/150‑1 (MQ=255)
                                                                              aaaaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGt                                                                   >  1:2204666/1‑150 (MQ=255)
                                                                                aaTTTTCTCTGATCTTAACTTCTGCAAATTTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                 <  2:1563052/150‑1 (MQ=255)
                                                                                aaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                 <  2:1566638/150‑1 (MQ=255)
                                                                                aaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                 <  2:1566618/150‑1 (MQ=255)
                                                                                aaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                 <  2:1566496/150‑1 (MQ=255)
                                                                                aaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                 <  2:1566493/150‑1 (MQ=255)
                                                                                aaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                 <  2:1563191/150‑1 (MQ=255)
                                                                                aaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                 >  2:2234145/1‑150 (MQ=255)
                                                                                 aTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGc                                                                <  1:4311590/150‑1 (MQ=255)
                                                                                 aTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGc                                                                <  1:4311955/150‑1 (MQ=255)
                                                                                   tttCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTa                                                              <  1:1015723/150‑1 (MQ=255)
                                                                                    ttCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTAt                                                             <  2:1249680/150‑1 (MQ=255)
                                                                                    ttCTCTGATCTTAACTTCTGCAAATGTTAACTGCTAATTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTAt                                                             <  2:1249684/150‑1 (MQ=255)
                                                                                      ctctGATCTTAACTTCTGCAAATTTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGt                                                           <  1:55634/150‑1 (MQ=255)
                                                                                                aaCTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTa                                                 <  1:1152301/150‑1 (MQ=255)
                                                                                                aaCTACTGCAAATGTTAACTGATATTTTTATGCTAAAAATGGTTATCAAAACTCAATAACACATTTTTATTATCAATGAGTTATAGTAATTCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTTTGTCAGATAGTTa                                                 <  2:2605949/150‑1 (MQ=255)
                                                                                                    tCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAAATAGTAAAATTATGTTGCTATGTCAGATAGTTACGAt                                             <  1:4427477/150‑1 (MQ=255)
                                                                                                       gCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTa                                          <  1:4531476/150‑1 (MQ=255)
                                                                                                        cAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAAAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAg                                         >  2:4042855/1‑150 (MQ=255)
                                                                                                        cAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAAAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAg                                         >  2:1138494/1‑150 (MQ=255)
                                                                                                        cAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAAAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAg                                         >  2:1138499/1‑150 (MQ=255)
                                                                                                            tGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAATACGCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCAt                                     >  1:2371134/1‑150 (MQ=255)
                                                                                                                   tGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaacta                              >  1:3405969/1‑150 (MQ=255)
                                                                                                                       aTTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgc                          <  1:4277755/150‑1 (MQ=255)
                                                                                                                            tATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTgcgc                     >  1:3067876/1‑150 (MQ=255)
                                                                                                                                 tAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                >  2:1022000/1‑150 (MQ=255)
                                                                                                                                 tAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                <  1:3127199/150‑1 (MQ=255)
                                                                                                                                 tAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                <  1:3127246/150‑1 (MQ=255)
                                                                                                                                 tAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                >  2:1021956/1‑150 (MQ=255)
                                                                                                                                 tAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                >  2:1021967/1‑150 (MQ=255)
                                                                                                                                  aaaaaTGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATg               <  2:2345649/150‑1 (MQ=255)
                                                                                                                                  aaaaaTGGTTATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATg               <  2:1097154/150‑1 (MQ=255)
                                                                                                                                      aTGGTTATCAAAACTCAAAAACACATGTTTATATTCAATGAGTTATAGAATTGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCTGAGAGTTACGATTTAGTCATGTCACTCATGCTGCGCCGTATggggt           >  2:2590349/1‑148 (MQ=255)
                                                                                                                                      aTGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACTATTTAGTCATCTAACTAATGCTGCGCCATATGGGtt           >  2:2590432/1‑150 (MQ=255)
                                                                                                                                          ttATCAATACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCACATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGAc       >  1:3869098/1‑150 (MQ=255)
                                                                                                                                          ttATCAATACTCAAAAACACATGTTTATAATCAAGGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGAc       >  1:3869083/1‑150 (MQ=255)
                                                                                                                                          ttATCAAAACTCAAAAACACATGATTATAATCAATGAGTTATAGAAATTCTAAGTGCTAATCAGGTATATGCAAATTAGTAAAATTATGATGCTATGTCATATAGTTACTTTTTAGTCATCTAACTAAAGCTTCGCCATATGGGTTTGaa       >  1:3868988/1‑149 (MQ=255)
                                                                                                                                             tCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGa    <  1:975697/150‑1 (MQ=255)
                                                                                                                                               aaaaCTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAg  >  1:2981419/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTTACGTTAATATAAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAG  >  CP009273/1206604‑1206896

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: