Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,007,918 G→T 6.5% R167S (CGC→AGC)  yqeC ← putative selenium‑dependent hydroxylase accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,007,9180GT 6.5% 861.0 / 22.8 228R167S (CGC→AGC) yqeCputative selenium‑dependent hydroxylase accessory protein
Reads supporting (aligned to +/- strand):  ref base G (123/90);  new base T (3/12);  total (126/102)
Fisher's exact test for biased strand distribution p-value = 6.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 231 total aligned reads displayed.

AAAAACGGTTAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAAT  >  CP009273/3007771‑3008065
                                                                                                                                                   |                                                                                                                                                   
ataaaCGGTTTATGTACCAGACTCGCGGAAAACCTTGCTGTACGATTTTAAACGCCCCATGAGGGTTGCGATCCAGCGCAACGATATCGCTCTGTTTCAAGGTTGCATCAGGTGTTAACCGGGTAATATTTGAAATATGCGACCAGCg                                                                                                                                                     <  2:307947/146‑1 (MQ=255)
aaaaaCGGTTAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCg                                                                                                                                                     <  2:748919/148‑1 (MQ=255)
 aaaaCGGTTAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGa                                                                                                                                                    >  1:3815738/1‑148 (MQ=255)
 aaaaCGGTTAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGa                                                                                                                                                    <  2:2865093/148‑1 (MQ=255)
    aCGGTTAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCTATgg                                                                                                                                                 >  2:176096/1‑148 (MQ=255)
       gTTAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACa                                                                                                                                              >  1:1155952/1‑148 (MQ=255)
       gTTAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACa                                                                                                                                              <  2:3035028/148‑1 (MQ=255)
        ttAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACAt                                                                                                                                             >  1:3489954/1‑148 (MQ=255)
         tAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACAtt                                                                                                                                            >  1:3196318/1‑148 (MQ=255)
          aaTGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCTATGGACAttt                                                                                                                                           <  2:2931187/148‑1 (MQ=255)
           aTGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACAtttt                                                                                                                                          >  2:1929458/1‑148 (MQ=255)
             gAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGATGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGTCCAGCGATGGACATTTTCt                                                                                                                                        >  2:71611/1‑148 (MQ=255)
             gAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCt                                                                                                                                        <  2:3410980/148‑1 (MQ=255)
                ccAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCtgtg                                                                                                                                     <  1:1758634/148‑1 (MQ=255)
                ccAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCtgtg                                                                                                                                     <  1:2029469/148‑1 (MQ=255)
                 cAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCTATGGACATTTTCTGTGc                                                                                                                                    <  1:777224/148‑1 (MQ=255)
                    aCTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCa                                                                                                                                 >  2:3757780/1‑148 (MQ=255)
                     cTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCAc                                                                                                                                <  1:2015845/148‑1 (MQ=255)
                     cTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCAc                                                                                                                                <  1:1462373/148‑1 (MQ=255)
                     cTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCAc                                                                                                                                >  2:2751568/1‑148 (MQ=255)
                     cTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCAc                                                                                                                                >  2:3399084/1‑148 (MQ=255)
                          cGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTcg                                                                                                                           >  2:1901515/1‑148 (MQ=255)
                            aCAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTcgcg                                                                                                                         >  1:1570380/1‑148 (MQ=255)
                            aCAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTcgcg                                                                                                                         >  2:820223/1‑148 (MQ=255)
                             cAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCgcgc                                                                                                                        <  2:1102590/148‑1 (MQ=255)
                             cAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCgcgc                                                                                                                        >  2:2014357/1‑148 (MQ=255)
                                ccTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCa                                                                                                                     >  1:1611959/1‑148 (MQ=255)
                                ccTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCa                                                                                                                     >  2:3803961/1‑148 (MQ=255)
                                ccTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCa                                                                                                                     >  1:3670966/1‑148 (MQ=255)
                                  ttGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCaaa                                                                                                                   <  1:2970150/148‑1 (MQ=255)
                                  ttGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCaaa                                                                                                                   <  1:258664/148‑1 (MQ=255)
                                  ttGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCaaa                                                                                                                   <  1:1065303/148‑1 (MQ=255)
                                   tGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGTACATTTTCTGTGCTCACTTTCGCGCCCAAAg                                                                                                                  >  1:3053580/1‑148 (MQ=255)
                                   tGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAg                                                                                                                  >  1:2108456/1‑148 (MQ=255)
                                        tACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGt                                                                                                             <  2:3135863/148‑1 (MQ=255)
                                        tACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGt                                                                                                             >  1:429528/1‑148 (MQ=255)
                                         aCGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGtc                                                                                                            >  1:2151480/1‑148 (MQ=255)
                                         aCGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGtc                                                                                                            >  2:2160333/1‑148 (MQ=255)
                                          cGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGtcc                                                                                                           <  2:1408649/148‑1 (MQ=255)
                                            tttttAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGtcctc                                                                                                         <  1:227319/148‑1 (MQ=255)
                                             ttttAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGtcctcc                                                                                                        >  1:611287/1‑148 (MQ=255)
                                              tttAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTccc                                                                                                       <  2:57615/148‑1 (MQ=255)
                                                tAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCAt                                                                                                     <  2:633664/148‑1 (MQ=255)
                                                   aCGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATcac                                                                                                  <  1:2833910/148‑1 (MQ=255)
                                                   aCGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATcac                                                                                                  >  1:705042/1‑148 (MQ=255)
                                                   aCGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATcac                                                                                                  >  1:952726/1‑148 (MQ=255)
                                                    cGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATcaca                                                                                                 <  1:3151370/148‑1 (MQ=255)
                                                          cTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATc                                                                                           <  1:852352/148‑1 (MQ=255)
                                                              gggTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACgc                                                                                       <  1:1146264/148‑1 (MQ=255)
                                                               ggTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACgca                                                                                      <  1:3399101/148‑1 (MQ=255)
                                                               ggTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACgca                                                                                      <  1:3364935/148‑1 (MQ=255)
                                                                  ggCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACgcagca                                                                                   <  1:93487/148‑1 (MQ=255)
                                                                  ggCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACgcagca                                                                                   >  2:219591/1‑148 (MQ=255)
                                                                    cGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAAc                                                                                 <  2:601597/148‑1 (MQ=255)
                                                                      aaCCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTg                                                                               >  2:1694822/1‑148 (MQ=255)
                                                                      aaCCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTg                                                                               >  1:1869261/1‑148 (MQ=255)
                                                                              cAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCACCTGCTTTTggg                                                                       >  2:2446139/1‑148 (MQ=255)
                                                                              cAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTggg                                                                       >  1:2554693/1‑148 (MQ=255)
                                                                                 cGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTAt                                                                    >  1:1321474/1‑148 (MQ=255)
                                                                                    gaTCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCa                                                                 >  1:3054552/1‑148 (MQ=255)
                                                                                     aTCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCac                                                                >  2:3727888/1‑147 (MQ=255)
                                                                                        gCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGt                                                             >  1:3582986/1‑148 (MQ=255)
                                                                                            aGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGt                                                         >  1:3774276/1‑148 (MQ=255)
                                                                                             gTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTg                                                        <  1:1547982/148‑1 (MQ=255)
                                                                                             gTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTg                                                        <  1:1547648/148‑1 (MQ=255)
                                                                                                     gTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAg                                                >  2:1780844/1‑148 (MQ=255)
                                                                                                         cATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGgcgc                                            >  1:3209459/1‑148 (MQ=255)
                                                                                                         cATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGgcgc                                            <  1:3347836/148‑1 (MQ=255)
                                                                                                          aTCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCt                                           <  1:2469618/148‑1 (MQ=255)
                                                                                                               gtgtTAACCCGGTAATATCAGCAAACTGCGACCAGCTATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTaa                                      <  1:2285144/148‑1 (MQ=255)
                                                                                                                 gtTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACg                                    >  1:3082330/1‑148 (MQ=255)
                                                                                                                   tAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGc                                  >  2:544892/1‑148 (MQ=255)
                                                                                                                       ccGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATGAGGCGCTTTTAACGGCATTc                              <  1:2842243/148‑1 (MQ=255)
                                                                                                                       ccGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTc                              <  2:312002/148‑1 (MQ=255)
                                                                                                                         ggTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCa                            >  2:3811391/1‑148 (MQ=255)
                                                                                                                           tAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAcg                          >  2:3563729/1‑148 (MQ=255)
                                                                                                                           tAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAcg                          >  2:3503080/1‑148 (MQ=255)
                                                                                                                           tAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAcg                          <  1:1979644/148‑1 (MQ=255)
                                                                                                                             atatCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAcgcg                        <  2:1155745/148‑1 (MQ=255)
                                                                                                                             atatCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAcgcg                        <  2:2573701/148‑1 (MQ=255)
                                                                                                                              tatCAGCAAACTGCGACCAGCTATGGACATGTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGa                       >  2:1807938/1‑148 (MQ=255)
                                                                                                                              tatCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCCGGCGCTTTTAACGGCATTCCACGCGa                       >  2:183941/1‑148 (MQ=255)
                                                                                                                              tatCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGa                       >  2:3359151/1‑148 (MQ=255)
                                                                                                                               atCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAg                      <  2:2151477/148‑1 (MQ=255)
                                                                                                                                tCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGc                     <  2:3490295/148‑1 (MQ=255)
                                                                                                                                tCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGc                     >  2:3815786/1‑148 (MQ=255)
                                                                                                                                   gCAAACTGCGACCAGCTATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCAt                  <  1:12356/148‑1 (MQ=255)
                                                                                                                                   gCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGccct                  >  1:1301390/1‑146 (MQ=255)
                                                                                                                                       aCTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGc              >  2:765852/1‑148 (MQ=255)
                                                                                                                                         tGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCtt            >  2:53001/1‑148 (MQ=255)
                                                                                                                                          gCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTc           >  2:1863943/1‑148 (MQ=255)
                                                                                                                                          gCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTc           >  1:542119/1‑148 (MQ=255)
                                                                                                                                            gACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCaa         <  2:1547883/148‑1 (MQ=255)
                                                                                                                                             aCCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAAt        <  2:952726/148‑1 (MQ=255)
                                                                                                                                              ccAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCTAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTGTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATg       >  1:3670552/1‑148 (MQ=255)
                                                                                                                                              ccAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATg       >  2:1459904/1‑148 (MQ=255)
                                                                                                                                                 gCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATgaga    <  2:1300926/148‑1 (MQ=255)
                                                                                                                                                  cGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATgaga  >  2:2583633/1‑147 (MQ=255)
                                                                                                                                                  cGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGaa   >  2:3717843/1‑148 (MQ=255)
                                                                                                                                                  cGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGaa   >  2:2583505/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAAACGGTTAATGAACCAGACTCGCCGACAACCTTGCGGTACGTTTTTAAACGCCCCCTGAGGGTGGCGAACCAGCGCAACGAGATCGCTCAGTTGCAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCCATCTGCTTCAATGAGAAT  >  CP009273/3007771‑3008065

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