Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,892,304 C→A 90.6% P150P (CCG→CCT yieK ← putative 6‑phosphogluconolactonase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,892,3040CA90.6% 863.1 / 79.1 318P150P (CCG→CCTyieKputative 6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand):  ref base C (20/10);  new base A (149/139);  total (169/149)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 318 total aligned reads displayed.

GCTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAG  >  CP009273/3892158‑3892449
                                                                                                                                                  |                                                                                                                                                 
gcTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCgg                                                                                                                                                 >  2:1815361/1‑149 (MQ=255)
 cTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAggg                                                                                                                                                >  2:724782/1‑149 (MQ=255)
 cTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAggg                                                                                                                                                <  2:3297949/149‑1 (MQ=255)
 cTTTACCCGCGCCGCTGACAATTATCAGCAAGTGTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAggg                                                                                                                                                <  2:3917964/149‑1 (MQ=255)
     aCCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACt                                                                                                                                            <  1:1067602/149‑1 (MQ=255)
        cgcgCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCa                                                                                                                                         <  2:577470/149‑1 (MQ=255)
           gccgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACgg                                                                                                                                      >  1:4379987/1‑149 (MQ=255)
            ccgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGt                                                                                                                                     >  2:3588242/1‑149 (MQ=255)
            ccgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTAGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGt                                                                                                                                     >  2:2930133/1‑149 (MQ=255)
             cgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGtt                                                                                                                                    <  1:3611304/149‑1 (MQ=255)
              gcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGttt                                                                                                                                   <  2:1571809/149‑1 (MQ=255)
              gcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGttt                                                                                                                                   >  1:3000276/1‑149 (MQ=255)
                tGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGc                                                                                                                                 >  2:1034629/1‑149 (MQ=255)
                tGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGc                                                                                                                                 >  1:3804083/1‑149 (MQ=255)
                 taCAATTATCAGCAAGTGTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCt                                                                                                                                <  2:281792/148‑1 (MQ=255)
                 gACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCt                                                                                                                                >  1:4122543/1‑149 (MQ=255)
                 gACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCt                                                                                                                                <  1:3404108/149‑1 (MQ=255)
                  aCAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTc                                                                                                                               >  1:434007/1‑149 (MQ=255)
                   cAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCa                                                                                                                              >  2:1876013/1‑149 (MQ=255)
                     aTTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATg                                                                                                                            <  2:4738427/149‑1 (MQ=255)
                     aTTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATg                                                                                                                            >  2:2873104/1‑149 (MQ=255)
                       tATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGa                                                                                                                          >  1:4267692/1‑149 (MQ=255)
                       tATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGa                                                                                                                          >  1:3893380/1‑149 (MQ=255)
                        aTCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGac                                                                                                                         >  1:4627987/1‑148 (MQ=255)
                        aTCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGaa                                                                                                                         >  1:2590581/1‑149 (MQ=255)
                        aTCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGaa                                                                                                                         >  1:4052350/1‑149 (MQ=255)
                        aTCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGaa                                                                                                                         >  1:1357090/1‑149 (MQ=255)
                          cAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAAt                                                                                                                       <  2:114495/149‑1 (MQ=255)
                           aGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATg                                                                                                                      >  1:454473/1‑149 (MQ=255)
                           aGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATg                                                                                                                      >  1:1741960/1‑149 (MQ=255)
                           aGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATg                                                                                                                      >  1:948091/1‑149 (MQ=255)
                            gCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATgg                                                                                                                     >  2:2390516/1‑149 (MQ=255)
                            gCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATgg                                                                                                                     <  2:481611/149‑1 (MQ=255)
                              aaGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGt                                                                                                                   >  2:1317257/1‑149 (MQ=255)
                              aaGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGt                                                                                                                   <  1:1267129/149‑1 (MQ=255)
                                gTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCg                                                                                                                 <  2:2801323/149‑1 (MQ=255)
                                 tttttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGAGAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGt                                                                                                                >  2:2022248/1‑149 (MQ=255)
                                  ttttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTa                                                                                                               <  2:3241646/149‑1 (MQ=255)
                                  ttttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTa                                                                                                               <  2:1650188/149‑1 (MQ=255)
                                    ttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATCTCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTAtt                                                                                                             <  1:876977/149‑1 (MQ=255)
                                    ttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTAtt                                                                                                             <  2:407082/149‑1 (MQ=255)
                                    ttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTAtt                                                                                                             <  1:1319873/149‑1 (MQ=255)
                                       gccgccATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCgg                                                                                                          >  1:2953994/1‑149 (MQ=255)
                                        ccgccATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGc                                                                                                         <  2:1690184/149‑1 (MQ=255)
                                           ccATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCaaa                                                                                                      >  2:2136240/1‑149 (MQ=255)
                                           ccATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCaaa                                                                                                      <  2:2221620/149‑1 (MQ=255)
                                            cATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAAt                                                                                                     <  2:866579/149‑1 (MQ=255)
                                              tgatgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTg                                                                                                   >  1:282771/1‑149 (MQ=255)
                                               gatgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGc                                                                                                  >  2:2468927/1‑149 (MQ=255)
                                                        cGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAAt                                                                                         <  1:4747081/149‑1 (MQ=255)
                                                        cGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAAt                                                                                         <  2:3785961/149‑1 (MQ=255)
                                                           ccccATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGAc                                                                                      <  1:1801018/149‑1 (MQ=255)
                                                           ccccATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGAc                                                                                      <  1:4190938/149‑1 (MQ=255)
                                                            cccATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGAcc                                                                                     <  2:2413659/149‑1 (MQ=255)
                                                                tCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCt                                                                                 >  1:1129493/1‑149 (MQ=255)
                                                                tCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCt                                                                                 >  1:3483378/1‑149 (MQ=255)
                                                                 cGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTg                                                                                <  1:2232114/149‑1 (MQ=255)
                                                                 cGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTg                                                                                >  1:3453702/1‑149 (MQ=255)
                                                                   tCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCa                                                                              >  2:2104720/1‑149 (MQ=255)
                                                                     aCGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCCTCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCAcc                                                                            <  2:4151578/149‑1 (MQ=255)
                                                                      cGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCt                                                                           >  1:264521/1‑149 (MQ=255)
                                                                        tAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTaa                                                                         <  1:643337/149‑1 (MQ=255)
                                                                        tAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTaa                                                                         <  2:3200010/149‑1 (MQ=255)
                                                                        tAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTaa                                                                         <  1:2724806/149‑1 (MQ=255)
                                                                          gCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATc                                                                       >  2:442011/1‑149 (MQ=255)
                                                                           cTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATcc                                                                      <  2:3470568/149‑1 (MQ=255)
                                                                            tGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGACATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAATTGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATccc                                                                     <  1:1235151/149‑1 (MQ=255)
                                                                             gTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCa                                                                    >  1:198826/1‑149 (MQ=255)
                                                                               cAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGc                                                                  <  2:138349/149‑1 (MQ=255)
                                                                                aGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCCTCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGca                                                                 <  2:4257818/149‑1 (MQ=255)
                                                                                aGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGTGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGca                                                                 >  2:1010426/1‑149 (MQ=255)
                                                                                     cTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCaccacc                                                            >  1:4467303/1‑149 (MQ=255)
                                                                                      tAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCAccacca                                                           >  1:23737/1‑149 (MQ=255)
                                                                                               gTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTc                                                  <  2:1474660/149‑1 (MQ=255)
                                                                                               gTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTc                                                  <  2:4513001/149‑1 (MQ=255)
                                                                                                   ccgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGcc                                              <  1:1592270/149‑1 (MQ=255)
                                                                                                         aaCTCGCCATGCGCCCCAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcg                                        >  1:2539338/1‑149 (MQ=255)
                                                                                                         aaCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcg                                        >  2:2363514/1‑149 (MQ=255)
                                                                                                          aCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcgc                                       >  1:3653274/1‑149 (MQ=255)
                                                                                                          aCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcgc                                       >  1:2156377/1‑149 (MQ=255)
                                                                                                             cGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGcc                                    <  1:3710400/149‑1 (MQ=255)
                                                                                                             cGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGcc                                    <  2:2594082/149‑1 (MQ=255)
                                                                                                              gCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCa                                   >  2:537725/1‑149 (MQ=255)
                                                                                                                 aTGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGtt                                <  2:3504892/149‑1 (MQ=255)
                                                                                                                    cgcCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCt                             <  1:2355217/149‑1 (MQ=255)
                                                                                                                     gcCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTg                            <  1:1472971/149‑1 (MQ=255)
                                                                                                                       cacaATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGAt                          <  1:3123580/149‑1 (MQ=255)
                                                                                                                        acaATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATc                         >  2:3063319/1‑149 (MQ=255)
                                                                                                                        acaATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATc                         >  2:2682655/1‑149 (MQ=255)
                                                                                                                           atatCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATg                      <  1:578223/149‑1 (MQ=255)
                                                                                                                             atCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCt                    <  1:3326679/149‑1 (MQ=255)
                                                                                                                             atCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCt                    <  1:48031/149‑1 (MQ=255)
                                                                                                                             atCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAAGGGGTCGTATTCGGCAAATTGCCACAAACATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCt                    <  2:4036940/149‑1 (MQ=255)
                                                                                                                                 aCCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTcgcg                >  2:1940747/1‑149 (MQ=255)
                                                                                                                                       tctcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGtt          >  2:1700669/1‑149 (MQ=255)
                                                                                                                                          cccTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATc       <  1:2481852/149‑1 (MQ=255)
                                                                                                                                            cTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCaa     >  2:4781741/1‑149 (MQ=255)
                                                                                                                                            cTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCaa     <  2:2197653/149‑1 (MQ=255)
                                                                                                                                             tGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAAt    >  2:669296/1‑149 (MQ=255)
                                                                                                                                              ggATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATa   <  2:1752974/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAG  >  CP009273/3892158‑3892449

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: