Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 379,966 G→A 12.6% S36L (TCG→TTG)  yaiS ← putative PIG‑L family deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273379,9660GA12.6% 951.2 / 100.1 309S36L (TCG→TTG) yaiSputative PIG‑L family deacetylase
Reads supporting (aligned to +/- strand):  ref base G (108/162);  new base A (21/18);  total (129/180)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 313 total aligned reads displayed.

ATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTATACTACAAGCACAACG  >  CP009273/379818‑380114
                                                                                                                                                    |                                                                                                                                                    
aTTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCg                                                                                                                                                      >  1:3613734/1‑149 (MQ=255)
   aaGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCAATg                                                                                                                                                   >  2:2391063/1‑149 (MQ=255)
    aGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATgc                                                                                                                                                  <  1:2049270/149‑1 (MQ=255)
         aaTAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCcac                                                                                                                                             <  2:2008540/149‑1 (MQ=255)
          aTAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCcaca                                                                                                                                            >  1:3294088/1‑149 (MQ=255)
            aGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGc                                                                                                                                          >  1:4537227/1‑149 (MQ=255)
             gTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGcc                                                                                                                                         >  2:1068655/1‑149 (MQ=255)
                 ggtggCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAAt                                                                                                                                     <  2:1428086/149‑1 (MQ=255)
                   tggCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTc                                                                                                                                   <  1:573073/149‑1 (MQ=255)
                   tggCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTc                                                                                                                                   >  2:1103210/1‑149 (MQ=255)
                   tggCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTc                                                                                                                                   <  1:2385226/149‑1 (MQ=255)
                     gCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCta                                                                                                                                 <  2:2873122/149‑1 (MQ=255)
                       aCCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCtata                                                                                                                               >  2:839318/1‑149 (MQ=255)
                        cccTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTatat                                                                                                                              >  1:4797571/1‑149 (MQ=255)
                        cccTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTatat                                                                                                                              >  2:4572891/1‑149 (MQ=255)
                         ccTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATc                                                                                                                             <  2:2257614/149‑1 (MQ=255)
                          cTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCg                                                                                                                            >  2:4864977/1‑149 (MQ=255)
                            aaTATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACTGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAATCGATGCGCCACAGCCTAATTCTATATCGTc                                                                                                                          >  2:4112218/1‑149 (MQ=255)
                             atatCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCa                                                                                                                         <  1:816716/149‑1 (MQ=255)
                              tatCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAg                                                                                                                        <  2:2494985/149‑1 (MQ=255)
                              tatCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAg                                                                                                                        <  2:3982103/149‑1 (MQ=255)
                              tatCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAg                                                                                                                        <  1:3344984/149‑1 (MQ=255)
                               atCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAgg                                                                                                                       >  2:5197503/1‑149 (MQ=255)
                                tCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGa                                                                                                                      <  2:897399/149‑1 (MQ=255)
                                tCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGa                                                                                                                      <  1:335649/149‑1 (MQ=255)
                                 cTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATACATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGAt                                                                                                                     <  2:4683898/149‑1 (MQ=255)
                                  tttAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCTTTTTTGCGCAAGACGAGCAAGCAATGCGCCACAGCCTAATTCTATATCGTCAGGAtg                                                                                                                    >  2:1713943/1‑149 (MQ=255)
                                  tttAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAAGTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGAtg                                                                                                                    <  2:2374574/149‑1 (MQ=255)
                                             gcgcGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAg                                                                                                         <  2:3435122/149‑1 (MQ=255)
                                               gcgATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCa                                                                                                       <  2:4433439/149‑1 (MQ=255)
                                                 gATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAg                                                                                                     <  2:4921545/149‑1 (MQ=255)
                                                 gATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAg                                                                                                     <  1:302570/149‑1 (MQ=255)
                                                   ttcttcATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGa                                                                                                   >  2:3974875/1‑149 (MQ=255)
                                                   ttcttcATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCAGGTAGTCATCACCACGGCTGCTATATAAATTCCTTTTTGAGCAAGCCGTGCAAGCGATGCGACACATCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGa                                                                                                   <  1:1024485/149‑1 (MQ=255)
                                                    tcttcATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGAt                                                                                                  >  2:2167094/1‑149 (MQ=255)
                                                    tcttcATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGAt                                                                                                  >  2:2845012/1‑149 (MQ=255)
                                                       tcATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATAcc                                                                                               <  1:2734163/149‑1 (MQ=255)
                                                       tcATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATAcc                                                                                               >  1:918025/1‑149 (MQ=255)
                                                       tcATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATAcc                                                                                               >  1:692093/1‑149 (MQ=255)
                                                         aTGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCct                                                                                             <  2:175677/149‑1 (MQ=255)
                                                           ggCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCcttt                                                                                           <  1:1634226/149‑1 (MQ=255)
                                                            gCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCctttc                                                                                          <  1:3176332/149‑1 (MQ=255)
                                                              gatcgatTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCctttctt                                                                                        >  2:1622004/1‑149 (MQ=255)
                                                               atcgatTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCctttcttt                                                                                       >  1:4711764/1‑149 (MQ=255)
                                                                  gatTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCctttcatccgg                                                                                    >  2:1844501/1‑143 (MQ=255)
                                                                  gatTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTAt                                                                                    <  1:3943065/149‑1 (MQ=255)
                                                                  gatTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTAt                                                                                    <  1:3484863/149‑1 (MQ=255)
                                                                   attattCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTAtt                                                                                   <  2:3133406/149‑1 (MQ=255)
                                                                    ttattCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCTTTTTTGCGCAAGACGAGCAAGCAATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTAttt                                                                                  >  2:2535387/1‑149 (MQ=255)
                                                                    ttattCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTAttt                                                                                  >  1:4183032/1‑149 (MQ=255)
                                                                      attCCATCTGTGACATAGTTGCCGGTAGTCATCAACAAGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCTTCAGGATGTGCGACAATAGCAAGGATACCCTTTCTTTTATTTGc                                                                                <  1:3752247/149‑1 (MQ=255)
                                                                       ttCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCGAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGcc                                                                               >  1:4115647/1‑149 (MQ=255)
                                                                          cATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCgaa                                                                            <  1:1739315/149‑1 (MQ=255)
                                                                            tCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCgaaga                                                                          >  2:46597/1‑149 (MQ=255)
                                                                            tCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCgaaga                                                                          >  1:5015643/1‑149 (MQ=255)
                                                                              tgtgCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGaaa                                                                        >  2:2750319/1‑149 (MQ=255)
                                                                                  ccAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAaggag                                                                    >  1:4294417/1‑149 (MQ=255)
                                                                                  ccAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAaggag                                                                    <  1:1269923/149‑1 (MQ=255)
                                                                                     gagTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGc                                                                 <  1:2891015/149‑1 (MQ=255)
                                                                                          gCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAAt                                                            >  1:1138521/1‑149 (MQ=255)
                                                                                           ccGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATc                                                           <  1:1370285/149‑1 (MQ=255)
                                                                                            cGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCt                                                          >  2:5163379/1‑149 (MQ=255)
                                                                                                aGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGAAGGGCTGAATCTaaaa                                                      <  1:4011869/149‑1 (MQ=255)
                                                                                                   catcaCCACGGCTGCGATATAAATTCTTTTTTGCGCAAGACGAGCAAGCAATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCt                                                   >  1:4308896/1‑149 (MQ=255)
                                                                                                      caccacGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTGTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTAt                                                >  2:1308936/1‑149 (MQ=255)
                                                                                                         cacGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCa                                             >  1:427407/1‑149 (MQ=255)
                                                                                                           cGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACt                                           <  1:659273/149‑1 (MQ=255)
                                                                                                                gCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATg                                      <  2:2605029/149‑1 (MQ=255)
                                                                                                                 cGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGa                                     <  1:1891170/149‑1 (MQ=255)
                                                                                                                   atatAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGaca                                   >  1:3820565/1‑149 (MQ=255)
                                                                                                                   atatAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGaca                                   <  2:3299996/149‑1 (MQ=255)
                                                                                                                   atatAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGaca                                   >  1:976407/1‑149 (MQ=255)
                                                                                                                     ataAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACt                                 >  1:1885691/1‑149 (MQ=255)
                                                                                                                      taAATTCCTTTTTGCGCAAGACGAGCAAGCAATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTc                                <  2:2081653/149‑1 (MQ=255)
                                                                                                                       aaaTTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTcc                               <  1:1858719/149‑1 (MQ=255)
                                                                                                                        aaTTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCa                              <  2:4241907/149‑1 (MQ=255)
                                                                                                                        aaTTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCa                              <  2:3013647/149‑1 (MQ=255)
                                                                                                                         aTTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCAt                             >  2:2663437/1‑149 (MQ=255)
                                                                                                                         aTTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCAt                             <  1:4930515/149‑1 (MQ=255)
                                                                                                                         aTTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCAt                             <  1:1676580/149‑1 (MQ=255)
                                                                                                                          ttCTTTTTTGCGCAAGACGAGCAAGCAATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCAtt                            >  1:4435312/1‑149 (MQ=255)
                                                                                                                          ttCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCAtt                            <  1:5121882/149‑1 (MQ=255)
                                                                                                                           tcttttttGCGCAAGACGAGCAAGCAATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCAttt                           >  1:2163444/4‑149 (MQ=255)
                                                                                                                             cTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATttta                         <  2:477755/149‑1 (MQ=255)
                                                                                                                             cTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATttta                         >  1:3562191/1‑149 (MQ=255)
                                                                                                                             cTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATttta                         <  2:3531247/149‑1 (MQ=255)
                                                                                                                              tttttGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATtttat                        <  2:2909649/149‑1 (MQ=255)
                                                                                                                               ttttGCGCAAGACGAGCAAGCAATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATtttatt                       <  2:3794446/149‑1 (MQ=255)
                                                                                                                                tttGCGCAAGACGAGCAAGCAATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATtttattt                      <  1:4934581/149‑1 (MQ=255)
                                                                                                                                 ttGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATtttattta                     <  2:133025/149‑1 (MQ=255)
                                                                                                                                  tGCGCAAGACGAGCAAGCAATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTttatttat                    <  2:1479181/149‑1 (MQ=255)
                                                                                                                                     gcAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTtattat                 <  2:3719664/149‑1 (MQ=255)
                                                                                                                                        aGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTAtact              >  1:2537719/1‑149 (MQ=255)
                                                                                                                                         gACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTAtacta             <  2:3184582/149‑1 (MQ=255)
                                                                                                                                            gAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTATACTACaa          >  2:598684/1‑149 (MQ=255)
                                                                                                                                             aGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTATACTACAAg         >  2:295959/1‑149 (MQ=255)
                                                                                                                                                aaGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTATACTACAAGcac      <  2:3594217/149‑1 (MQ=255)
                                                                                                                                                aaGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTATACTACAAGcac      <  1:1252241/149‑1 (MQ=255)
                                                                                                                                                 aGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTATACTACAAGcaca     >  1:3245846/1‑149 (MQ=255)
                                                                                                                                                    gATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTATACTACAAGCACAACg  >  1:2018095/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACGGCTGCGATATAAATTCCTTTTTGCGCAAGACGAGCAAGCGATGCGCCACAGCCTAATTCTATATCGTCAGGATGTGCGCCAATAGCAAGGATACCCTTTCTTTTATTTGCCGAAGAAAGGAGGGCTGAATCTAAAACCTTATCCACTTAATGACACTCCATTTTATTTATTATACTACAAGCACAACG  >  CP009273/379818‑380114

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: