Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,964,899 C→A 7.6% intergenic (+51/‑57) ygdR → / → tas DUF903 family verified lipoprotein/putative NADP(H)‑dependent aldo‑keto reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,964,8990CA 7.6% 891.8 / 34.6 238intergenic (+51/‑57)ygdR/tasDUF903 family verified lipoprotein/putative NADP(H)‑dependent aldo‑keto reductase
Reads supporting (aligned to +/- strand):  ref base C (118/101);  new base A (2/16);  total (121/117)
Fisher's exact test for biased strand distribution p-value = 4.06e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 292 total aligned reads displayed.

GAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGAC  >  CP009273/2964750‑2965048
                                                                                                                                                     |                                                                                                                                                     
gAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCtttcccttttc                                                                                                                                                       >  1:3539696/1‑150 (MQ=255)
gAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCtttcacttttc                                                                                                                                                       >  1:1120090/1‑150 (MQ=255)
   aTTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCtttcccttttccct                                                                                                                                                    >  1:900274/1‑150 (MQ=255)
    ttGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTttcccttttccctt                                                                                                                                                   <  1:1158497/150‑1 (MQ=255)
      gatgatgatACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTttcccttcc                                                                                                                                                 >  1:1824551/1‑150 (MQ=255)
       atgatgaaAACGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATCCAAATTAACCGTGAAGATGTTTAGAAAATTATAGAACGAAAACAAATAAGGTCAGCATCCGGCTGGCCTAAAGATAGTGCTCTTTCCCTTgtcccttccc                                                                                                                                                >  2:907341/1‑150 (MQ=255)
        tgatgatACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAGATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTtaccttccct                                                                                                                                               >  2:4725579/1‑150 (MQ=255)
        tgatgatACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTtcccttccct                                                                                                                                               <  1:66160/150‑1 (MQ=255)
         gatgatACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCTTCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCctc                                                                                                                                              <  2:3502118/150‑1 (MQ=255)
           tgatACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTACCTTCCCTCTg                                                                                                                                            <  2:737574/150‑1 (MQ=255)
             atACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGcc                                                                                                                                          <  2:2444459/150‑1 (MQ=255)
              tACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCa                                                                                                                                         <  2:359418/150‑1 (MQ=255)
               aCCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCAt                                                                                                                                        <  2:139217/150‑1 (MQ=255)
               aCCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCAt                                                                                                                                        <  2:4040158/150‑1 (MQ=255)
               aCCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCAt                                                                                                                                        >  1:2569950/1‑150 (MQ=255)
               aCCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCAt                                                                                                                                        <  2:2737205/150‑1 (MQ=255)
                ccGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCAtt                                                                                                                                       >  2:3156580/1‑150 (MQ=255)
                  gggCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCAtttt                                                                                                                                     >  1:4024109/1‑150 (MQ=255)
                  gggCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTCTCTCTTTCCCTTTTCCCTTCCCTCTGCCAtttt                                                                                                                                     <  2:3613271/150‑1 (MQ=255)
                   ggCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTACCTTCCCTCTGCCAttttt                                                                                                                                    <  2:701497/150‑1 (MQ=255)
                    gCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTta                                                                                                                                   <  1:5116879/150‑1 (MQ=255)
                     cTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTACCTTCCCTCTGCCATTTTtat                                                                                                                                  <  1:1439325/150‑1 (MQ=255)
                        gTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATCtt                                                                                                                               >  1:4576766/1‑150 (MQ=255)
                        gTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTACCTTCCCTCTGCCATTTTTATAtt                                                                                                                               <  2:1624599/150‑1 (MQ=255)
                            gTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCtt                                                                                                                           >  2:2257619/1‑150 (MQ=255)
                               aTCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATg                                                                                                                        >  1:3702466/1‑150 (MQ=255)
                                 cACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTc                                                                                                                      <  1:2721615/150‑1 (MQ=255)
                                   cGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGt                                                                                                                    <  2:1704128/150‑1 (MQ=255)
                                    gATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTg                                                                                                                   >  1:4411134/1‑150 (MQ=255)
                                       cAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGAtt                                                                                                                >  1:489986/1‑150 (MQ=255)
                                       cAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGAtt                                                                                                                <  2:2907232/150‑1 (MQ=255)
                                       cAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGAtt                                                                                                                >  2:4779960/1‑150 (MQ=255)
                                       cAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGGGAtt                                                                                                                >  1:3077291/1‑150 (MQ=255)
                                        aGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATta                                                                                                               <  2:3399107/150‑1 (MQ=255)
                                         gCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATtat                                                                                                              >  1:4019012/1‑150 (MQ=255)
                                           aaGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATaa                                                                                                            >  2:2644164/1‑150 (MQ=255)
                                           aaGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATaa                                                                                                            >  2:1211085/1‑150 (MQ=255)
                                            aGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATaaa                                                                                                           <  1:3420909/150‑1 (MQ=255)
                                               tAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAg                                                                                                        >  1:3767535/1‑150 (MQ=255)
                                                  cgcgATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGaa                                                                                                     <  1:727495/150‑1 (MQ=255)
                                                  cgcgATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTAATGTCGTGATTATAAAAAGGaa                                                                                                     >  2:4711307/1‑150 (MQ=255)
                                                   gcgATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGaaa                                                                                                    <  1:3505151/150‑1 (MQ=255)
                                                    tgATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAAc                                                                                                   >  1:1862774/2‑150 (MQ=255)
                                                         cAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATACGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTa                                                                                              <  2:2409724/150‑1 (MQ=255)
                                                         cAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTa                                                                                              <  2:4953778/150‑1 (MQ=255)
                                                         cAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTAAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTATTATATTCCTTATGTCGTGATTATAAAAAGGAAAAGGCTa                                                                                              >  2:332545/1‑150 (MQ=255)
                                                          aaaTTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTAt                                                                                             >  2:64270/1‑150 (MQ=255)
                                                            aTTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGc                                                                                           >  2:2576828/1‑150 (MQ=255)
                                                             ttAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCa                                                                                          <  1:1583964/150‑1 (MQ=255)
                                                             ttAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTACCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCa                                                                                          <  1:3290287/150‑1 (MQ=255)
                                                              tAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCaa                                                                                         >  1:4385389/1‑150 (MQ=255)
                                                                   gTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCa                                                                                    <  2:3185771/150‑1 (MQ=255)
                                                                   gTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCa                                                                                    >  1:1910884/1‑150 (MQ=255)
                                                                   gTGATGATGTTTCCAAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTACCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCa                                                                                    <  1:3449677/150‑1 (MQ=255)
                                                                    tgatgaTGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAc                                                                                   >  1:3258382/1‑150 (MQ=255)
                                                                     gatgatGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAcc                                                                                  <  1:1590143/150‑1 (MQ=255)
                                                                      atgatgTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCg                                                                                 >  1:4349054/1‑150 (MQ=255)
                                                                      atgatgTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCg                                                                                 >  1:3068014/1‑150 (MQ=255)
                                                                        gatgTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGta                                                                               >  2:3085669/1‑150 (MQ=255)
                                                                              tCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCACTATCACCGTATAcccc                                                                         >  2:2014541/1‑150 (MQ=255)
                                                                              tCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCACTATCACCGTATAcccc                                                                         >  2:4878579/1‑150 (MQ=255)
                                                                                 cAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCcaca                                                                      <  1:3512782/150‑1 (MQ=255)
                                                                                  aaaTTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAg                                                                     <  1:286714/150‑1 (MQ=255)
                                                                                    attattGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGtt                                                                   <  2:1299911/150‑1 (MQ=255)
                                                                                    attattGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGtt                                                                   >  2:4449045/1‑150 (MQ=255)
                                                                                     ttattGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTc                                                                  <  2:3395331/150‑1 (MQ=255)
                                                                                     ttattGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTc                                                                  <  2:5042545/150‑1 (MQ=255)
                                                                                         tGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTACCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTg                                                              <  2:560572/150‑1 (MQ=255)
                                                                                           aaCGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGa                                                            <  2:3692986/150‑1 (MQ=255)
                                                                                                 aaCAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCGTTTTGCCTTCCCTCTGAGATTTTTATATTGGTTATGTCGTGATTATAAAATGTAAACGGCTATTCAATATCACCGTATACCCCTCAGTTCGCTGGAAGTCAg                                                      >  2:3937561/1‑150 (MQ=255)
                                                                                                    aaaTAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCAc                                                   <  1:4105522/150‑1 (MQ=255)
                                                                                                     aaTAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACg                                                  <  2:1977685/150‑1 (MQ=255)
                                                                                                            tCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGc                                           >  1:1250226/1‑150 (MQ=255)
                                                                                                              aGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCtt                                         >  2:474979/1‑150 (MQ=255)
                                                                                                                 aTCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGc                                      >  2:4596715/1‑150 (MQ=255)
                                                                                                                 aTCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGc                                      >  2:1740799/1‑150 (MQ=255)
                                                                                                                        tGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGa                               <  2:3756637/150‑1 (MQ=255)
                                                                                                                              tAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTg                         >  2:4267232/1‑150 (MQ=255)
                                                                                                                                aGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGt                       >  1:4048020/1‑150 (MQ=255)
                                                                                                                                  aTTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTga                     <  2:2044348/150‑1 (MQ=255)
                                                                                                                                  aTTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTga                     <  2:4529201/150‑1 (MQ=255)
                                                                                                                                  aTTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTga                     >  1:2163411/1‑150 (MQ=255)
                                                                                                                                       tctcTTTCCATTTTCAATTCCATCTGCCTTTTTTATATTCCTTATTTCGTTATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCATCTCGCTGGGGCTTGGCACGATGACGTTTGGTgaacaga                <  1:999846/150‑1 (MQ=255)
                                                                                                                                        ctctTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTgaacagaa               <  1:1021639/150‑1 (MQ=255)
                                                                                                                                         tctTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTgaacagaac              <  1:214528/150‑1 (MQ=255)
                                                                                                                                          ctTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTgaacagaaca             >  2:140986/1‑150 (MQ=255)
                                                                                                                                          ctTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTgaacagaaca             >  1:3575650/1‑150 (MQ=255)
                                                                                                                                          ctTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTgaacagaaca             >  1:1700130/1‑150 (MQ=255)
                                                                                                                                           tttcccttttccctTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGaacagaacag            >  2:4801481/1‑150 (MQ=255)
                                                                                                                                             tcccttttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCg          <  2:4166835/150‑1 (MQ=255)
                                                                                                                                             tcccgtttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCg          <  2:168951/150‑1 (MQ=255)
                                                                                                                                               ccttttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGaa        >  1:1456525/1‑150 (MQ=255)
                                                                                                                                                cttttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAg       >  1:4161524/1‑150 (MQ=255)
                                                                                                                                                 ttttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGc      >  2:2595979/1‑150 (MQ=255)
                                                                                                                                                 ttttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGc      >  1:2266853/1‑150 (MQ=255)
                                                                                                                                                 ttttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAAAGGCTATTCAATATCAACGTATAACCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGc      >  1:4238378/1‑150 (MQ=255)
                                                                                                                                                  tttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGcc     <  2:111085/150‑1 (MQ=255)
                                                                                                                                                  tttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGcc     <  1:480834/150‑1 (MQ=255)
                                                                                                                                                     cccttccctCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGAc  <  2:3847395/150‑1 (MQ=255)
                                                                                                                                                     cccttccctCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGAc  <  2:3008465/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAGGTAACGCGATGCAAATTAACCGTGATGATGTTTCGCAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGAC  >  CP009273/2964750‑2965048

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: