Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 658,096 C→A 100% M1M (ATG→ATT) † ybeD ← UPF0250 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273658,0960CA96.3% 619.2 / 12.2 204M1M (ATG→ATT) †ybeDUPF0250 family protein
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (81/115);  total (85/119)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 204 total aligned reads displayed.

TTGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAAT  >  CP009273/657951‑658240
                                                                                                                                                 |                                                                                                                                                
ttGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAAt                                                                                                                                                <  2:791526/148‑1 (MQ=255)
ttGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAAt                                                                                                                                                <  1:1237730/148‑1 (MQ=255)
ttGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAAt                                                                                                                                                >  2:122042/1‑148 (MQ=255)
 tGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGGTTGGTTTTAATg                                                                                                                                               >  2:1642861/1‑148 (MQ=255)
   gTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTc                                                                                                                                             >  1:3444689/1‑148 (MQ=255)
    ttttACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCa                                                                                                                                            <  2:1645227/148‑1 (MQ=255)
    ttttACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCa                                                                                                                                            <  1:1382313/148‑1 (MQ=255)
     tttACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAg                                                                                                                                           <  1:218832/148‑1 (MQ=255)
        aCCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTc                                                                                                                                        >  1:2601788/1‑148 (MQ=255)
        aCCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTc                                                                                                                                        >  1:382894/1‑148 (MQ=255)
          cGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCg                                                                                                                                      <  1:220685/148‑1 (MQ=255)
               gggTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCTCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGcgta                                                                                                                                 <  1:2467912/148‑1 (MQ=255)
               gggTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGcgta                                                                                                                                 >  2:2514991/1‑148 (MQ=255)
                   gtAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGcgtaacgt                                                                                                                             <  2:2869266/148‑1 (MQ=255)
                     aGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGcgtaacgtaa                                                                                                                           <  1:404879/148‑1 (MQ=255)
                       tCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAAtt                                                                                                                         >  1:2902144/1‑148 (MQ=255)
                          ccTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGGAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTaaa                                                                                                                      >  2:3029299/1‑148 (MQ=255)
                           cTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAAt                                                                                                                     <  2:1530588/148‑1 (MQ=255)
                             ggCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAg                                                                                                                   >  2:2259196/1‑148 (MQ=255)
                                gcATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCaa                                                                                                                <  1:32517/148‑1 (MQ=255)
                                 cATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAAc                                                                                                               >  2:2362646/1‑148 (MQ=255)
                                  aTGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACt                                                                                                              >  1:2665240/1‑148 (MQ=255)
                                   tGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTc                                                                                                             >  2:323493/1‑148 (MQ=255)
                                   tGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTc                                                                                                             <  2:25776/148‑1 (MQ=255)
                                   tGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTc                                                                                                             <  1:2022258/148‑1 (MQ=255)
                                    ggCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTcc                                                                                                            <  1:76993/148‑1 (MQ=255)
                                    ggCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTcc                                                                                                            <  2:761852/148‑1 (MQ=255)
                                       gcTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTATGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGc                                                                                                         <  2:3324456/148‑1 (MQ=255)
                                       gcTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGc                                                                                                         <  2:3185407/148‑1 (MQ=255)
                                           tACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGa                                                                                                     >  2:1182902/1‑148 (MQ=255)
                                           tACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGa                                                                                                     >  2:522990/1‑148 (MQ=255)
                                             ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAg                                                                                                   <  2:2833597/148‑1 (MQ=255)
                                             ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAg                                                                                                   <  2:3284390/148‑1 (MQ=255)
                                                 ttCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCgg                                                                                               >  1:1096837/1‑148 (MQ=255)
                                                     accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGt                                                                                           <  1:2060820/148‑1 (MQ=255)
                                                        accTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttt                                                                                        >  1:2262668/1‑148 (MQ=255)
                                                        accTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttt                                                                                        <  1:2818275/148‑1 (MQ=255)
                                                           tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                                                                                     <  1:122042/148‑1 (MQ=255)
                                                            gATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAc                                                                                    <  2:1518745/148‑1 (MQ=255)
                                                                aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                                                                                <  1:1370342/148‑1 (MQ=255)
                                                                 aCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                                                                               >  2:2258627/1‑148 (MQ=255)
                                                                 aCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                                                                               >  1:272665/1‑148 (MQ=255)
                                                                  ccAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGc                                                                              >  1:75219/1‑148 (MQ=255)
                                                                   cAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCt                                                                             <  2:2139129/148‑1 (MQ=255)
                                                                    aGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCtt                                                                            >  2:3466894/1‑148 (MQ=255)
                                                                     gCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTa                                                                           <  2:1745021/148‑1 (MQ=255)
                                                                     gCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTa                                                                           <  1:3036281/148‑1 (MQ=255)
                                                                       tCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                                                                         <  1:211473/148‑1 (MQ=255)
                                                                        cAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                                                                        >  1:2391450/1‑148 (MQ=255)
                                                                         aGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGtat                                                                       <  2:589578/148‑1 (MQ=255)
                                                                         aGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGtat                                                                       <  1:1308940/148‑1 (MQ=255)
                                                                           gTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTatat                                                                     <  1:1013791/148‑1 (MQ=255)
                                                                           gTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTatat                                                                     >  2:340175/1‑148 (MQ=255)
                                                                               cGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgggg                                                                 >  1:988141/1‑148 (MQ=255)
                                                                                gCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGa                                                                <  1:1114270/148‑1 (MQ=255)
                                                                                 ccTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc                                                               >  2:2625852/1‑148 (MQ=255)
                                                                                   tGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACgg                                                             <  1:3314124/148‑1 (MQ=255)
                                                                                      cccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGaaa                                                          <  2:2946751/148‑1 (MQ=255)
                                                                                      cccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGaaa                                                          <  2:2244368/148‑1 (MQ=255)
                                                                                      cccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTATGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGaaa                                                          <  2:2069727/148‑1 (MQ=255)
                                                                                       ccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAt                                                         <  1:3123241/148‑1 (MQ=255)
                                                                                        cATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAtt                                                        >  2:2956970/1‑148 (MQ=255)
                                                                                        cATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAtt                                                        >  1:810608/1‑148 (MQ=255)
                                                                                         aTAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTa                                                       >  1:1610132/1‑148 (MQ=255)
                                                                                           aaCTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTaca                                                     <  1:2362646/148‑1 (MQ=255)
                                                                                             cTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACt                                                   >  1:2528924/1‑148 (MQ=255)
                                                                                                tGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTc                                                <  1:1695270/148‑1 (MQ=255)
                                                                                                 gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCa                                               >  2:39583/1‑148 (MQ=255)
                                                                                                  taagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCaa                                              <  2:2336502/148‑1 (MQ=255)
                                                                                                     gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtg                                           >  2:2662042/1‑148 (MQ=255)
                                                                                                     gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtg                                           >  1:1174852/1‑148 (MQ=255)
                                                                                                      taaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtgt                                          >  1:2501817/1‑148 (MQ=255)
                                                                                                      taaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAAGGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtgt                                          <  2:3109305/148‑1 (MQ=255)
                                                                                                        aaaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGttt                                        <  2:527805/148‑1 (MQ=255)
                                                                                                        aaaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAATGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGttt                                        <  1:1485358/148‑1 (MQ=255)
                                                                                                           ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAt                                     >  1:181751/1‑148 (MQ=255)
                                                                                                            gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAtt                                    >  1:234839/1‑148 (MQ=255)
                                                                                                               tAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACttt                                                 >  1:2632288/1‑132 (MQ=255)
                                                                                                                aGGGAATTCACGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCGCCGGCGTAACGTAATTAAAGAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTa                                <  2:2810314/148‑1 (MQ=255)
                                                                                                                aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTa                                <  1:2505023/148‑1 (MQ=255)
                                                                                                                  ggAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaac                              <  2:1832487/148‑1 (MQ=255)
                                                                                                                  ggAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaac                              >  1:618010/1‑148 (MQ=255)
                                                                                                                      ttCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTATTTTTTaaccaaa                          >  1:2138353/1‑148 (MQ=255)
                                                                                                                      ttCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaa                          <  2:542074/148‑1 (MQ=255)
                                                                                                                      ttCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaa                          <  2:1201058/148‑1 (MQ=255)
                                                                                                                           gCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGt                     >  2:3264776/1‑148 (MQ=255)
                                                                                                                           gCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGt                     >  1:492627/1‑148 (MQ=255)
                                                                                                                            cAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTg                    >  2:2087025/1‑148 (MQ=255)
                                                                                                                               ttCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATg                 >  2:244662/1‑148 (MQ=255)
                                                                                                                               ttCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATg                 >  2:1331135/1‑148 (MQ=255)
                                                                                                                                 cGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGa               <  2:2918604/148‑1 (MQ=255)
                                                                                                                                   ttAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAAc             <  1:953212/148‑1 (MQ=255)
                                                                                                                                    tAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACa            >  2:3089084/1‑148 (MQ=255)
                                                                                                                                    tAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACa            <  1:2258627/148‑1 (MQ=255)
                                                                                                                                      aGTTTGGTTTTAATGTCATCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAtt          <  2:3266256/148‑1 (MQ=255)
                                                                                                                                        tttGGTTTTAATGTCAGCTCCGGCGTAACTTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaa        <  1:2956970/148‑1 (MQ=255)
                                                                                                                                           ggTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaattt     <  1:323493/148‑1 (MQ=255)
                                                                                                                                           ggTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaattt     >  2:2986881/1‑148 (MQ=255)
                                                                                                                                            gTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaattta    <  1:1653544/148‑1 (MQ=255)
                                                                                                                                            gTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaattta    <  1:3003781/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAAT  >  CP009273/657951‑658240

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: