Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,297,485 G→C 100% G44G (GGC→GGG alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,297,4850GC97.0% 996.8 / 16.3 251G44G (GGC→GGGalsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base C (126/117);  total (130/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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 252 total aligned reads displayed.

GCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCT  >  CP009273/4297339‑4297628
                                                                                                                                                  |                                                                                                                                               
gccgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGccc                                                                                                                                               >  1:2602693/1‑149 (MQ=255)
  cgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCgg                                                                                                                                             <  1:584031/149‑1 (MQ=255)
  cgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCgg                                                                                                                                             <  1:2869681/149‑1 (MQ=255)
  cgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCgg                                                                                                                                             <  1:2380564/149‑1 (MQ=255)
   gcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCggg                                                                                                                                            >  1:1033118/1‑149 (MQ=255)
   gcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCggg                                                                                                                                            >  1:1033070/1‑149 (MQ=255)
      gTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGc                                                                                                                                         >  2:3430114/1‑149 (MQ=255)
          aCGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATg                                                                                                                                     >  2:299048/1‑149 (MQ=255)
          aCGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATg                                                                                                                                     >  2:4203274/1‑149 (MQ=255)
           cGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGa                                                                                                                                    <  2:1657119/149‑1 (MQ=255)
              cAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACtt                                                                                                                                 <  1:2931601/149‑1 (MQ=255)
                  ttAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGc                                                                                                                             <  2:3680396/149‑1 (MQ=255)
                     ggCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGt                                                                                                                          >  2:1927392/1‑149 (MQ=255)
                     ggCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGt                                                                                                                          >  1:3934086/1‑149 (MQ=255)
                     ggCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGt                                                                                                                          >  2:1105313/1‑149 (MQ=255)
                     ggCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGt                                                                                                                          >  2:4125469/1‑149 (MQ=255)
                      gCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTc                                                                                                                         <  1:559002/149‑1 (MQ=255)
                       cGTAGAAAGAATGGTGCGTTGATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAACGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTcc                                                                                                                        <  1:4134167/149‑1 (MQ=255)
                         tAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGc                                                                                                                      >  1:2023000/1‑149 (MQ=255)
                         tAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCAGc                                                                                                                      >  1:3830181/1‑149 (MQ=255)
                             aataatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCtttt                                                                                                                  >  2:1054908/1‑149 (MQ=255)
                              ataatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCttttt                                                                                                                 <  1:3627776/149‑1 (MQ=255)
                                aatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCttttttt                                                                                                               <  1:4677608/149‑1 (MQ=255)
                                   ggTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGc                                                                                                            <  2:3553983/149‑1 (MQ=255)
                                        gTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATGATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAgtgt                                                                                                       >  1:3165273/1‑149 (MQ=255)
                                        gTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAgtgt                                                                                                       >  2:4633122/1‑149 (MQ=255)
                                         tttATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTa                                                                                                      <  1:4031464/149‑1 (MQ=255)
                                          ttATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAg                                                                                                     <  1:4173226/149‑1 (MQ=255)
                                            aTCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCg                                                                                                   >  2:862447/1‑149 (MQ=255)
                                             tCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATCCCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGt                                                                                                  >  2:3333608/1‑149 (MQ=255)
                                                ttACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTc                                                                                               >  2:4358227/1‑149 (MQ=255)
                                                     gACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTCGCGTTTCACCtt                                                                                          >  2:3060487/1‑149 (MQ=255)
                                                     gACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCtt                                                                                          >  1:1574886/1‑149 (MQ=255)
                                                      aCCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTc                                                                                         <  2:4061839/149‑1 (MQ=255)
                                                      aCCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTc                                                                                         <  1:3171172/149‑1 (MQ=255)
                                                       ccAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTct                                                                                        <  1:3551805/149‑1 (MQ=255)
                                                       ccAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTAct                                                                                        <  2:1531108/149‑1 (MQ=255)
                                                         aGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTctgc                                                                                      >  2:2734810/1‑149 (MQ=255)
                                                              cGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTcc                                                                                 <  1:906108/149‑1 (MQ=255)
                                                                       cATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACaaa                                                                        <  1:2014217/149‑1 (MQ=255)
                                                                        aTCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACaaaa                                                                       <  2:3447015/149‑1 (MQ=255)
                                                                         tCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAg                                                                      <  2:2023019/149‑1 (MQ=255)
                                                                         tCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAg                                                                      >  2:4098848/1‑149 (MQ=255)
                                                                         tCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAg                                                                      <  2:318946/149‑1 (MQ=255)
                                                                             cAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGa                                                                  >  1:407226/1‑149 (MQ=255)
                                                                                   aTGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCg                                                            <  1:2444064/149‑1 (MQ=255)
                                                                                     gACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGtt                                                          >  2:151680/1‑149 (MQ=255)
                                                                                      aCAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTg                                                         <  2:4528817/149‑1 (MQ=255)
                                                                                       cAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCGTCTTCTGTCCGCAGAAAAAAAAGGAAATGCGTTGc                                                        >  1:821471/1‑149 (MQ=255)
                                                                                        agcgagcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATCTGCGTTGcc                                                       >  2:3832930/1‑149 (MQ=255)
                                                                                            agcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCa                                                   <  2:634519/149‑1 (MQ=255)
                                                                                             gcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTATCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCat                                                  >  1:432050/1‑149 (MQ=255)
                                                                                             gcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCat                                                  >  1:1890729/1‑149 (MQ=255)
                                                                                             gcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCat                                                  <  1:2209307/149‑1 (MQ=255)
                                                                                               ggtaccgcgcCTGAGTTGCTCGTCAATCATTTCGCCGATCCCCGACACCAGGCCGGGAGCAATGACTGCGGCGGTCCGCTGTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGGTGCCCCCagat                                                <  2:4519425/143‑1 (MQ=255)
                                                                                                ttAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATc                                               >  1:329399/1‑149 (MQ=255)
                                                                                                ttAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATc                                               >  1:507300/1‑149 (MQ=255)
                                                                                                  aaaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCa                                             >  2:4715866/1‑149 (MQ=255)
                                                                                                   aaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCAc                                            <  2:2247483/149‑1 (MQ=255)
                                                                                                   aaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCAc                                            >  2:1630010/1‑149 (MQ=255)
                                                                                                   aaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCAc                                            <  1:2255365/149‑1 (MQ=255)
                                                                                                   aaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCAc                                            >  1:3244541/1‑149 (MQ=255)
                                                                                                      cgcCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGcc                                         <  2:1756811/149‑1 (MQ=255)
                                                                                                         cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACCAAAGCGGATATGCGTTGCCCCCCTATCCACGCCCGc                                      >  2:3463618/1‑149 (MQ=255)
                                                                                                         cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGc                                      <  2:2692417/149‑1 (MQ=255)
                                                                                                         cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGc                                      <  2:4021550/149‑1 (MQ=255)
                                                                                                         cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGc                                      >  2:1908312/1‑149 (MQ=255)
                                                                                                         cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGc                                      >  1:4762816/1‑149 (MQ=255)
                                                                                                         cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGc                                      <  2:8054/149‑1 (MQ=255)
                                                                                                          tGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCt                                     <  1:2317560/149‑1 (MQ=255)
                                                                                                          tGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCt                                     <  2:3035046/149‑1 (MQ=255)
                                                                                                             gTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacg                                  >  2:258024/1‑149 (MQ=255)
                                                                                                                gCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacgacg                               >  1:2181692/1‑149 (MQ=255)
                                                                                                                  tCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGtt                             <  1:2705148/149‑1 (MQ=255)
                                                                                                                    gTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTAt                           <  1:3381753/149‑1 (MQ=255)
                                                                                                                    gTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTAt                           >  1:3973088/1‑149 (MQ=255)
                                                                                                                     tCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATg                          >  2:4105152/1‑149 (MQ=255)
                                                                                                                     tCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATg                          >  2:3349752/1‑149 (MQ=255)
                                                                                                                      cAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGc                         <  1:1083665/149‑1 (MQ=255)
                                                                                                                      cAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGc                         <  2:208623/149‑1 (MQ=255)
                                                                                                                        aTCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTg                       >  2:182279/1‑149 (MQ=255)
                                                                                                                         tCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGt                      <  1:745659/149‑1 (MQ=255)
                                                                                                                            tttCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGtttt                   <  2:3859471/149‑1 (MQ=255)
                                                                                                                            gttCGCCGAGACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGtttt                   <  2:2219738/148‑1 (MQ=255)
                                                                                                                             ttCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGttttg                  <  1:3584906/149‑2 (MQ=255)
                                                                                                                               cGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGc                <  1:3758269/149‑1 (MQ=255)
                                                                                                                               cGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGc                <  2:2438298/149‑1 (MQ=255)
                                                                                                                                gCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCa               <  2:2770495/149‑1 (MQ=255)
                                                                                                                                gCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCa               >  1:4369136/1‑149 (MQ=255)
                                                                                                                                gCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCa               <  1:4298413/149‑1 (MQ=255)
                                                                                                                                 ccGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCAt              <  2:3327266/149‑1 (MQ=255)
                                                                                                                                  cGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTGCATGa            <  1:3901442/149‑1 (MQ=255)
                                                                                                                                    aTACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAg           <  2:428844/149‑1 (MQ=255)
                                                                                                                                    aTACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAg           >  2:836201/1‑149 (MQ=255)
                                                                                                                                    aTACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAg           <  2:714695/149‑1 (MQ=255)
                                                                                                                                     tACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAgg          >  1:3001121/1‑149 (MQ=255)
                                                                                                                                      aCCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGc         <  2:3109965/149‑1 (MQ=255)
                                                                                                                                       cccGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCt        <  2:2452189/149‑1 (MQ=255)
                                                                                                                                            cacCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTAcc   >  1:3011752/1‑149 (MQ=255)
                                                                                                                                             acCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCt  >  1:1453398/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCT  >  CP009273/4297339‑4297628

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: