Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 233,829 C→A 19.9% S3* (TCA→TAA)  yafT → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273233,8290CA19.9% 783.4 / 129.4 296S3* (TCA→TAA) yafTlipoprotein
Reads supporting (aligned to +/- strand):  ref base C (137/100);  new base A (25/34);  total (162/134)
Fisher's exact test for biased strand distribution p-value = 4.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 300 total aligned reads displayed.

AGAATAAATTTATATTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTT  >  CP009273/233683‑233976
                                                                                                                                                  |                                                                                                                                                   
    tAAATTTATATTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTaaaaaa                                                                                                                                                >  1:2223263/1‑148 (MQ=255)
    tAAATTTATATTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCaaaaa                                                                                                                                                >  1:3091610/1‑148 (MQ=255)
    tAAATTTATATTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCaaaaa                                                                                                                                                >  1:1554760/1‑148 (MQ=255)
    tAAATTTATATTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCaaaaa                                                                                                                                                >  1:3358314/1‑148 (MQ=255)
       aTTTATATTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGc                                                                                                                                             >  1:1262742/1‑148 (MQ=255)
        tttATATTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCt                                                                                                                                            <  2:2265646/148‑1 (MQ=255)
         ttATATTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCtt                                                                                                                                           <  2:1969303/148‑1 (MQ=255)
            tatTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCACAAAAGCTttgt                                                                                                                                        >  2:3323086/1‑148 (MQ=255)
             atTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTttgtt                                                                                                                                       >  2:2125585/1‑148 (MQ=255)
               tGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGc                                                                                                                                     <  2:103032/148‑1 (MQ=255)
                 aTTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCat                                                                                                                                   >  1:2838249/1‑148 (MQ=255)
                    aaaTGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATAtg                                                                                                                                <  1:1669100/148‑1 (MQ=255)
                       tgaatgTATATTTCAAATTTATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATAtgtgt                                                                                                                             <  2:3476152/148‑1 (MQ=255)
                       tgaatgTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATAtgtgt                                                                                                                             <  2:1335958/148‑1 (MQ=255)
                        gaatgTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATAtgtgtg                                                                                                                            >  2:596153/1‑148 (MQ=255)
                        gaatgTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATAtgtgtg                                                                                                                            >  2:2841774/1‑148 (MQ=255)
                        gaatgTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATAtgtgtg                                                                                                                            >  2:2989047/1‑148 (MQ=255)
                        gaatgTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATAtgtgtg                                                                                                                            >  2:2959390/1‑148 (MQ=255)
                        gaaggTATATTTCAAATTGAATTTGTTTGTTATTATTTAAGTTAGGTATATAATTATAGTCCGTTATCAATGCTAAATATTATAATCATAATTACAGGAGATGGAGTGTCCAATTATGAATTAAACAAAGCTTTGATGCATAtgtgtg                                                                                                                            >  2:788566/1‑148 (MQ=255)
                          atgTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGtt                                                                                                                          >  1:2757254/1‑148 (MQ=255)
                           tgTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTa                                                                                                                         >  1:1410174/1‑148 (MQ=255)
                            gTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTAt                                                                                                                        >  2:64642/1‑148 (MQ=255)
                              atatTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATtc                                                                                                                      >  1:2096227/1‑148 (MQ=255)
                              atatTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATtc                                                                                                                      >  2:2495892/1‑148 (MQ=255)
                                atTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATtctc                                                                                                                    >  1:1643154/1‑148 (MQ=255)
                                  ttCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAATCTTTGTTGCATATGTGTGTTATTCTCgc                                                                                                                  >  1:2070617/1‑148 (MQ=255)
                                  ttCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCgc                                                                                                                  <  2:526894/148‑1 (MQ=255)
                                   tCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCgct                                                                                                                 >  2:1366358/1‑148 (MQ=255)
                                       aTTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCtt                                                                                                             >  1:250409/1‑148 (MQ=255)
                                       aTTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCtt                                                                                                             >  1:3489048/1‑148 (MQ=255)
                                       aTTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCtt                                                                                                             >  2:66531/1‑148 (MQ=255)
                                        ttGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTg                                                                                                            >  2:2363133/1‑148 (MQ=255)
                                         tGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGc                                                                                                           >  2:918171/1‑148 (MQ=255)
                                           aTTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAg                                                                                                         >  2:822278/1‑148 (MQ=255)
                                            ttttGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAgg                                                                                                        <  1:3331168/148‑1 (MQ=255)
                                            ttttGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAgg                                                                                                        >  1:1087073/1‑148 (MQ=255)
                                            ttttGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAgg                                                                                                        <  1:2255958/148‑1 (MQ=255)
                                                   tgttATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGcc                                                                                                 >  1:833284/1‑148 (MQ=255)
                                                    gttATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCct                                                                                                <  1:2859248/148‑1 (MQ=255)
                                                    gttATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCct                                                                                                <  2:3393138/148‑1 (MQ=255)
                                                    gttATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCct                                                                                                <  1:341922/148‑1 (MQ=255)
                                                       attattTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTg                                                                                             >  2:2891464/1‑148 (MQ=255)
                                                        ttattTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGa                                                                                            <  1:700014/148‑1 (MQ=255)
                                                          attTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAAt                                                                                          >  2:1251221/1‑148 (MQ=255)
                                                           tttAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAAtc                                                                                         >  1:2614146/1‑148 (MQ=255)
                                                            ttAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAAtct                                                                                        >  1:1147400/1‑148 (MQ=255)
                                                               aGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAAtcttct                                                                                     >  2:392724/1‑148 (MQ=255)
                                                               aGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAAtcttct                                                                                     <  2:3706146/148‑1 (MQ=255)
                                                                gTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTa                                                                                    <  2:309560/148‑1 (MQ=255)
                                                                 tGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATGTTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTAt                                                                                   <  1:2365921/148‑1 (MQ=255)
                                                                   aGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtg                                                                                 <  1:2989047/148‑1 (MQ=255)
                                                                    ggTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtga                                                                                <  1:791521/148‑1 (MQ=255)
                                                                     gTATATAATTAGAGTCCGTTACCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgat                                                                               >  2:217253/1‑148 (MQ=255)
                                                                         ataATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGaaa                                                                           <  2:3358314/148‑1 (MQ=255)
                                                                         ataATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGaaa                                                                           >  2:2536895/1‑148 (MQ=255)
                                                                            aTTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAaaga                                                                        >  2:2235477/1‑148 (MQ=255)
                                                                                  gTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGaaaa                                                                  <  1:1341959/148‑1 (MQ=255)
                                                                                  gTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGaaaa                                                                  <  2:3532147/148‑1 (MQ=255)
                                                                                   tCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGaaaaa                                                                 <  2:1452570/148‑1 (MQ=255)
                                                                                         aTCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCa                                                           <  1:181233/148‑1 (MQ=255)
                                                                                           cAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAg                                                         >  2:2610093/1‑148 (MQ=255)
                                                                                           cAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAg                                                         >  2:3727877/1‑148 (MQ=255)
                                                                                           cAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGaaaaa                                                                 <  1:1177095/140‑1 (MQ=255)
                                                                                           aaattatAAATATTATAATATTTTTGACAGGCGAGGTAGTGTCAAATTATGAATTAAAAAAAGATTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAg                                                         <  2:3791045/142‑1 (MQ=255)
                                                                                              tGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTca                                                      <  2:1578567/148‑1 (MQ=255)
                                                                                              tGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTca                                                      <  2:1932989/148‑1 (MQ=255)
                                                                                               gCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgatg                                                                              >  1:3163031/1‑123 (MQ=255)
                                                                                               gCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGACAAGAAGAAAAAAGCTCAAGTcac                                                     >  2:403786/1‑148 (MQ=255)
                                                                                               gCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTcac                                                     >  1:1216796/1‑148 (MQ=255)
                                                                                                 tAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTcacac                                                   >  1:558599/1‑148 (MQ=255)
                                                                                                  aaaTATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTcacaca                                                  <  1:2891234/148‑1 (MQ=255)
                                                                                                   aaTATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACaa                                                 <  2:2738236/148‑1 (MQ=255)
                                                                                                       ttCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTa                                             <  1:3620960/148‑1 (MQ=255)
                                                                                                       ttCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTa                                             >  1:3793605/1‑148 (MQ=255)
                                                                                                        tCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTaa                                            <  2:1227515/148‑1 (MQ=255)
                                                                                                        tCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTaa                                            <  2:1227167/148‑1 (MQ=255)
                                                                                                        tCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTATGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTaa                                            <  1:2551346/148‑1 (MQ=255)
                                                                                                          tAATCATTATGACAGGCGAGGGAGTGTCCAATTACGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAata                                          <  2:1475484/148‑1 (MQ=255)
                                                                                                           aaTCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAatat                                         >  1:1033248/1‑148 (MQ=255)
                                                                                                             tCATTATGACAGGCGAGGTAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAATAATATta                                       <  2:3573644/148‑1 (MQ=255)
                                                                                                                ttATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaata                                    <  1:1586748/148‑1 (MQ=255)
                                                                                                                ttATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaata                                    >  1:2673324/1‑148 (MQ=255)
                                                                                                                 tATGACAGGCGAGGGAGTGTCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaataa                                   >  1:2303397/1‑148 (MQ=255)
                                                                                                                   tGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATaaaa                                 <  1:1200013/148‑1 (MQ=255)
                                                                                                                     aCAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAac                               >  2:1669215/1‑148 (MQ=255)
                                                                                                                      cAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAaca                              <  1:918171/148‑1 (MQ=255)
                                                                                                                       aGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAacac                             <  1:3068825/148‑1 (MQ=255)
                                                                                                                           gAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTccc                         <  2:3231331/148‑1 (MQ=255)
                                                                                                                              ggAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAg                      <  1:378968/148‑1 (MQ=255)
                                                                                                                              ggAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAg                      >  1:1930049/1‑148 (MQ=255)
                                                                                                                                    tCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTg                >  2:1802403/1‑148 (MQ=255)
                                                                                                                                    tCCAATTATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTg                <  1:43050/148‑1 (MQ=255)
                                                                                                                                          tATGAATTAAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTgacagac          >  1:1743078/1‑148 (MQ=255)
                                                                                                                                            tGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACaa        >  1:965730/1‑148 (MQ=255)
                                                                                                                                            tGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACaa        >  2:1168581/1‑148 (MQ=255)
                                                                                                                                              aaTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAg      >  2:590598/1‑148 (MQ=255)
                                                                                                                                              aaTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAg      <  1:1441167/148‑1 (MQ=255)
                                                                                                                                               aTTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGa     <  2:1087027/148‑1 (MQ=255)
                                                                                                                                                 tCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAtt   >  2:2052329/1‑148 (MQ=255)
                                                                                                                                                  cAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGAttt  <  1:66531/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGAATAAATTTATATTGATTAAATGAATGTATATTTCAAATTGATTTTGTTTGTTATTATTTAAGTGAGGTATATAATTAGAGTCCGTTATCAATGCTAAATATTCTAATCATTATGACAGGCGAGGGAGTGTCCAATTATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACTGACAGACAAAGATTT  >  CP009273/233683‑233976

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: