Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,207,958 G→T 5.6% intergenic (‑499/+201) iraM ← / ← ycgX RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,207,9580GT 5.6% 1194.5 / 4.9 383intergenic (‑499/+201)iraM/ycgXRpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base G (116/228);  new base T (35/4);  total (151/232)
Fisher's exact test for biased strand distribution p-value = 8.34e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 386 total aligned reads displayed.

GAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATT  >  CP009273/1207809‑1208107
                                                                                                                                                     |                                                                                                                                                     
gagaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATg                                                                                                                                                       >  2:1209179/1‑150 (MQ=255)
 agaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGt                                                                                                                                                      >  2:3072024/1‑150 (MQ=255)
  gaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGtt                                                                                                                                                     >  2:2390446/1‑150 (MQ=255)
     tAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAc                                                                                                                                                  <  1:1877729/150‑1 (MQ=255)
     tAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAc                                                                                                                                                  <  1:1309930/150‑1 (MQ=255)
      aaGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAcc                                                                                                                                                 <  2:290068/150‑1 (MQ=255)
         gTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGtt                                                                                                                                              >  1:3077531/1‑150 (MQ=255)
         gTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGtt                                                                                                                                              <  2:2536589/150‑1 (MQ=255)
         gTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAAAATGTTAACCGtt                                                                                                                                              >  2:2330341/1‑150 (MQ=255)
            tCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTGTTTTTTTTAAATATTTTAAACGTtata                                                                                                                                           >  2:4935061/1‑150 (MQ=255)
            tCAATTTCTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTtata                                                                                                                                           <  1:3828375/150‑1 (MQ=255)
             cAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTCTaa                                                                                                                                          <  1:1023608/150‑1 (MQ=255)
              aaTTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaaa                                                                                                                                         <  2:3599951/150‑1 (MQ=255)
               aTTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaaaa                                                                                                                                        >  2:1954589/1‑150 (MQ=255)
                 ttttGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATTTTAAACGTTATAAAATa                                                                                                                                      >  1:511855/1‑150 (MQ=255)
                   ttGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAAc                                                                                                                                    <  2:4296615/150‑1 (MQ=255)
                    tGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACg                                                                                                                                   <  2:3216764/150‑1 (MQ=255)
                      aGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATACCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAACGTTATAAAATAAAGTa                                                                                                                                 >  1:1193564/1‑150 (MQ=255)
                        ttATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATc                                                                                                                               <  2:1395921/150‑1 (MQ=255)
                        ttATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATc                                                                                                                               <  2:3803098/150‑1 (MQ=255)
                          atataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCaa                                                                                                                             >  2:630426/1‑150 (MQ=255)
                           tataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAACGTTATAAAATAACGTATaaaa                                                                                                                            >  2:5061592/1‑150 (MQ=255)
                            ataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGTAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTTTCaaaa                                                                                                                           >  2:4618649/1‑150 (MQ=255)
                              aGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGt                                                                                                                         <  2:1070410/150‑1 (MQ=255)
                              aGGAATGATTTAAACTGTTCCTGTCTAAAATACATATAACCGTATGATGAATAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGt                                                                                                                         <  1:4809860/150‑1 (MQ=255)
                                gAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCa                                                                                                                       <  1:205385/150‑1 (MQ=255)
                                  aTGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAg                                                                                                                     <  1:1630502/150‑1 (MQ=255)
                                   tGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAACGTTATAAAATAACGTATCAAAATTCaatt                                                                                                                    >  2:3666012/1‑148 (MQ=255)
                                    gATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTg                                                                                                                   <  1:495534/150‑1 (MQ=255)
                                     aTTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGa                                                                                                                  <  1:2190310/150‑1 (MQ=255)
                                     aTTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAACGTTATAAAATAACGTATAAAAAGTCAAGTGa                                                                                                                  >  1:747088/1‑150 (MQ=255)
                                       ttAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATc                                                                                                                <  1:1224024/150‑1 (MQ=255)
                                         aaCCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATcac                                                                                                              >  2:867663/1‑150 (MQ=255)
                                          aCCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATcaca                                                                                                             <  2:1472065/150‑1 (MQ=255)
                                              gTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTGGTTTTTATAAATATGTTAAACGGTATAAAATAAAGGATCAAAAGTAAAGTGATCACAttta                                                                                                         >  1:3480622/1‑149 (MQ=255)
                                               ttCCTGGATAAAATACATATAACCTGATGATGAATAAACAAAAATACTTGTTCGTTAAGTATTGAAACGGACTTGTGGCTCGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCa                                                                                                        <  1:2958697/150‑1 (MQ=255)
                                                  cTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAat                                                                                                     <  1:4713241/150‑1 (MQ=255)
                                                    ggCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAatat                                                                                                   <  1:330406/150‑1 (MQ=255)
                                                    ggCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAatat                                                                                                   <  1:326165/150‑1 (MQ=255)
                                                       tAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCaa                                                                                                <  2:4849646/150‑1 (MQ=255)
                                                         aaaTACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGt                                                                                              <  1:1343753/150‑1 (MQ=255)
                                                              cATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATa                                                                                         <  1:4036737/150‑1 (MQ=255)
                                                               atatAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAg                                                                                        <  1:3564458/150‑1 (MQ=255)
                                                                  taACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTGGGTTTTATAAATATGTTAAACCTTATAAAAAAACGTATCAAAAGTCAAGTGATCACAATTAAAATATAAATTTTATAGTAt                                                                                     >  1:648801/1‑150 (MQ=255)
                                                                     ccGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAg                                                                                  <  2:1622499/150‑1 (MQ=255)
                                                                        gatgatGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTTTTAATATTTTAACCGTTATAAAATAAAGTATAAAAAGTAAAGTGATCACATTTAAAAAATAAATTTGATAttattattag                                                                               >  1:3417727/1‑146 (MQ=255)
                                                                        gatgatGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCt                                                                               <  2:5296144/150‑1 (MQ=255)
                                                                        gatgatGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTGGTTTTTTTAAAAATTTTAAACGATATAAAATAAAATATAAAAAGTCAATTGTTAACATTTAAAATATAAAGTTTGTAGTATTAGtgt                                                                               >  1:1490774/1‑148 (MQ=255)
                                                                          tgatgaCTAAACCAAAATACATGTGCGTTAAGTATTTAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTgg                                                                             <  2:2766828/150‑1 (MQ=255)
                                                                           gatgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAACCTTATAAAATAACGTATCAAAAGTCAATTGATAACATTTCAAATATAAAGTTGATATTATTAGtattgt                                                                            >  2:4096022/1‑145 (MQ=255)
                                                                                cTAAACCAAAATACCTGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGatta                                                                       <  1:269307/150‑1 (MQ=255)
                                                                                  aaaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGattatt                                                                     <  1:3970889/150‑1 (MQ=255)
                                                                                   aaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttattt                                                                    <  1:2524401/150‑1 (MQ=255)
                                                                                    aCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttattta                                                                   <  2:76438/150‑1 (MQ=255)
                                                                                     ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttatttat                                                                  <  2:3313723/150‑1 (MQ=255)
                                                                                       aaaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCCAGTGATCACATTTCAAATATCAAGTTGAGAGTATTAGTCTGGTGATTATTTATgg                                                                >  2:3478403/1‑150 (MQ=255)
                                                                                       aaaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATgg                                                                <  1:3117799/150‑1 (MQ=255)
                                                                                       aaaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAAAATGTTAAACGTTATAAAAAAAAGTATAAAAAGTAAAGTGATAAAATTTCAAATATAAAGTTGATATTATTATTATTGTGAttatttattg                                                                >  1:2232517/1‑148 (MQ=255)
                                                                                        aaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTGGTTTTTATAAATATTTTAACCGGTATAAAAAAAAGTATCAAAAGTCAAGTGTTCACAATTAAAATATCAAGTTGATAGTATTATTCTGGTGATTATTTATtgg                                                               >  1:838574/1‑150 (MQ=255)
                                                                                         aaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATAAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTATGGTGAttatttatttgt                                                              >  2:5233805/1‑146 (MQ=255)
                                                                                         aaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTGTTTTTTTTAAATAAGTTAAACGGTAAAAAAAAAAATAAAAAAACAAAAGTGATAAAATTAAAAATATAAAAATGAAAATATTAATAATGTGATTATTTATGGGt                                                              >  1:2601110/1‑150 (MQ=255)
                                                                                          gtACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTTACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTAAAATATAAAGTTGATAGTATTAGTATGGTGATTATTTATTGGTg                                                             >  2:2817417/2‑150 (MQ=255)
                                                                                           tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGa                                                            <  2:1037176/150‑1 (MQ=255)
                                                                                            aCATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAAAGTATCAAAATTCAATTGATAACATTTCAAATATCAAATTGATATTATTATTCTGGTGATTATTTATTGGttaa                                                           >  2:227232/1‑147 (MQ=255)
                                                                                              aTGTGCGTTTAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACaa                                                         <  2:4964136/150‑1 (MQ=255)
                                                                                               tgtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAAt                                                        <  1:1549550/150‑1 (MQ=255)
                                                                                               tgtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAAt                                                        <  2:698820/150‑1 (MQ=255)
                                                                                               tgtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAAt                                                        <  1:4561111/150‑1 (MQ=255)
                                                                                                gtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATa                                                       <  1:2361540/150‑1 (MQ=255)
                                                                                                gtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATa                                                       <  1:5009435/150‑1 (MQ=255)
                                                                                                gtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATa                                                       <  2:4342508/150‑1 (MQ=255)
                                                                                                 tgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaa                                                      <  1:1781056/150‑1 (MQ=255)
                                                                                                 tgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaa                                                      <  2:2607868/150‑1 (MQ=255)
                                                                                                  gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                                     <  1:3190023/150‑1 (MQ=255)
                                                                                                    gTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaaaa                                                   <  2:4433712/150‑1 (MQ=255)
                                                                                                    gTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAACGTTATAAAATAAAGTATCAAAAATCAAGTGATCACAATTAAAATATCAAGATGATATTATTAGTATGGTGATTATTTATGGGTGaaaaaaaaaa                                                   >  1:5258745/1‑150 (MQ=255)
                                                                                                       aaGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGAc                                                <  1:4452557/150‑1 (MQ=255)
                                                                                                         gTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAg                                              <  1:709825/150‑1 (MQ=255)
                                                                                                          tATTGTACCGGTCGTGTGGCACGGCCGTTGTTTTTATAAATCTGTTAACCGTTATAAAATATCGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTTTTAGTCTGGTGATTATTTATGGGTGACATTATAAAGTCAGt                                             <  1:2160009/150‑1 (MQ=255)
                                                                                                             tGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTAtt                                          <  1:2815346/150‑1 (MQ=255)
                                                                                                                aaCGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAat                                       <  1:1108376/150‑1 (MQ=255)
                                                                                                                   ggACGTGTGGCACGGCCGTTGTTTTTATTAATAATTTAAACGTTATAAAAAAAAGTATAAAAAATAAAGTGATAACAATTAAAAAAAAAAGTTGTTAGTATTATTATGGTTTTTATTTTTTTGTGTAAAAAAAAAAACAGTATTaataat                                    >  1:3249986/1‑150 (MQ=255)
                                                                                                                   ggACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatcat                                    <  1:3215083/150‑1 (MQ=255)
                                                                                                                    gACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAAGATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatcatc                                   <  2:393643/150‑1 (MQ=255)
                                                                                                                     aCGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATcc                                  <  1:4114005/150‑1 (MQ=255)
                                                                                                                        tgtgGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATa                               <  2:3854483/150‑1 (MQ=255)
                                                                                                                            gCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAATTAAAGTGATCACATTTCAAATATAAAGTTTATAGTATTAGTCTGGTGATTATTTATTGGTGACAATAAAAAGACAGTATTAATCATCCATAgaga                           >  2:1777173/1‑150 (MQ=255)
                                                                                                                             cACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGAt                          <  1:3752693/150‑1 (MQ=255)
                                                                                                                               cGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAg                        <  1:192830/150‑1 (MQ=255)
                                                                                                                                   cGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGtctc                    >  1:2855824/1‑150 (MQ=255)
                                                                                                                                    gttgttTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTctct                   >  2:296675/1‑150 (MQ=255)
                                                                                                                                             ataCATATGTTAACCGTTATACAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTa          <  1:425020/150‑1 (MQ=255)
                                                                                                                                               aaaTATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACTGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTAtt        >  1:4354253/1‑150 (MQ=255)
                                                                                                                                               aaaTATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTAtt        <  1:2107552/150‑1 (MQ=255)
                                                                                                                                               aaaTATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTAtt        >  2:2050961/1‑150 (MQ=255)
                                                                                                                                               aaaTATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTAtt        <  1:4938134/150‑1 (MQ=255)
                                                                                                                                                  tatGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTcc     >  2:3157362/1‑150 (MQ=255)
                                                                                                                                                  tatGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTcc     <  2:2563321/150‑1 (MQ=255)
                                                                                                                                                   atGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCa    <  2:512292/150‑1 (MQ=255)
                                                                                                                                                     gTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCAtt  <  2:4007762/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATT  >  CP009273/1207809‑1208107

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