Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,358 A→G 41.8% N41D (AAC→GAC) ‡ talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,3580AG41.8% 78.8 / 154.3 122N41D (AAC→GAC) ‡talBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base G (28/23);  total (63/59)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 123 total aligned reads displayed.

ACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCG  >  CP009273/8211‑8507
                                                                                                                                                   |                                                                                                                                                     
aCAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCTCAACCAACCCTTCTCTCATTCTTAAc                                                                                                                                                     <  1:905634/150‑1 (MQ=255)
  aaGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcgc                                                                                                                                                   <  2:947455/150‑1 (MQ=255)
         ttAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcgcagcg                                                                                                                                               >  1:442265/1‑147 (MQ=255)
            aaGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGAt                                                                                                                                          <  1:658031/149‑1 (MQ=255)
               agaAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGAGGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCg                                                                                                                                      >  1:1042122/1‑150 (MQ=255)
                  aaTACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGaa                                                                                                                                   <  2:833653/150‑1 (MQ=255)
                     aCTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATAc                                                                                                                                >  1:1052834/1‑150 (MQ=255)
                     aCTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATAc                                                                                                                                >  2:102263/1‑150 (MQ=255)
                      cTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATAcc                                                                                                                               <  1:750697/150‑1 (MQ=255)
                      cTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATAcc                                                                                                                               <  2:1028636/150‑1 (MQ=255)
                      acaaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTAcc                                                                                                                               >  1:315741/6‑150 (MQ=255)
                      acaaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTAcc                                                                                                                               >  1:428332/6‑150 (MQ=255)
                       caaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCg                                                                                                                              <  2:369418/146‑1 (MQ=255)
                       caaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCg                                                                                                                              <  1:468668/146‑1 (MQ=255)
                       caaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCg                                                                                                                              <  2:266040/146‑1 (MQ=255)
                        aaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                                                                                                                             <  1:388126/147‑1 (MQ=255)
                          cATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGt                                                                                                                             >  2:787570/1‑148 (MQ=255)
                           aTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAg                                                                                                                          >  1:258949/1‑150 (MQ=255)
                           aTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAg                                                                                                                          >  1:734995/1‑150 (MQ=255)
                             gACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGtt                                                                                                                        >  2:1126570/1‑150 (MQ=255)
                               cGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctga                                                                                                                      >  2:1038402/1‑150 (MQ=255)
                               cGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttga                                                                                                                      >  2:633300/1‑150 (MQ=255)
                                ggATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgat                                                                                                                     <  2:49848/150‑1 (MQ=255)
                                ggATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgat                                                                                                                     <  1:1075467/150‑1 (MQ=255)
                                 gACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttgatt                                                                                                                    >  2:1192639/1‑150 (MQ=255)
                                  ataaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgatt                                                                                                                    >  1:688317/3‑149 (MQ=255)
                                  aCAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttgatt                                                                                                                    >  1:817462/1‑149 (MQ=255)
                                   taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga                                                                                                                  >  2:1169343/2‑150 (MQ=255)
                                   taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga                                                                                                                  <  2:809958/149‑1 (MQ=255)
                                   taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga                                                                                                                  <  2:531907/149‑1 (MQ=255)
                                     aaTTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTtgattgat                                                                                                                 >  1:214340/1‑149 (MQ=255)
                                      aTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATtgacga                                                                                                               >  2:383958/1‑150 (MQ=255)
                                       ttGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTgatgat                                                                                                              <  2:632217/150‑1 (MQ=255)
                                          aCCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCt                                                                                                           <  2:868868/150‑1 (MQ=255)
                                                cTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGcc                                                                                                     <  2:421435/150‑1 (MQ=255)
                                                   cGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGcc                                                                                                     <  2:1129310/147‑1 (MQ=255)
                                                   cGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTgg                                                                                                  <  1:960182/150‑1 (MQ=255)
                                                   cGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTgg                                                                                                  >  2:487879/1‑150 (MQ=255)
                                                    gTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGAGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTggg                                                                                                 <  1:80216/150‑1 (MQ=255)
                                                    gTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTggg                                                                                                 >  1:1191317/1‑150 (MQ=255)
                                                    gTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTggg                                                                                                 >  1:564387/1‑150 (MQ=255)
                                                    gTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTgg                                                                                                  <  2:134679/149‑1 (MQ=255)
                                                      cAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGc                                                                                                >  2:409302/1‑149 (MQ=255)
                                                        gTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGctcg                                                                                             <  1:831030/150‑4 (MQ=255)
                                                        gTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGctcg                                                                                             <  2:30185/150‑4 (MQ=255)
                                                           caccacTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGctcgtag                                                                                          >  2:896458/1‑144 (MQ=255)
                                                           caccacTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGctcgtag                                                                                          >  2:615151/1‑144 (MQ=255)
                                                           caccacTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGctcgtag                                                                                          >  1:511431/1‑144 (MQ=255)
                                                           caccacTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGctcgtag                                                                                          >  2:875033/1‑144 (MQ=255)
                                                           caccacCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGaaa                                                                                            >  1:672247/1‑148 (MQ=255)
                                                           caccacCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAAc                                                                                           <  1:633300/149‑1 (MQ=255)
                                                            accaccGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAAcag                                                                                         <  2:1052834/150‑1 (MQ=255)
                                                             ccaccGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAAcagc                                                                                        <  2:785236/150‑1 (MQ=255)
                                                                ccGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCaga                                                                                     <  2:817462/150‑1 (MQ=255)
                                                                   tAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCa                                                                                  <  1:1026921/150‑1 (MQ=255)
                                                                           gACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATcgc                                                                          <  2:720596/150‑1 (MQ=255)
                                                                              acTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATcgcgcg                                                                       >  1:1142939/1‑150 (MQ=255)
                                                                                  gggACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCgcagc                                                                   <  1:725406/150‑1 (MQ=255)
                                                                                      tattgcgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCgcag                                                                    >  1:1132482/5‑145 (MQ=255)
                                                                                      tattgcgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCgcag                                                                    >  1:253692/5‑145 (MQ=255)
                                                                                      cATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGAt                                                               <  1:779517/150‑1 (MQ=255)
                                                                                       attgcgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATc                                                              >  2:955598/4‑150 (MQ=255)
                                                                                       aTCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATc                                                              <  1:918938/150‑1 (MQ=255)
                                                                                           cgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctcg                                                          <  2:315741/150‑5 (MQ=255)
                                                                                           cgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctc                                                           >  1:510553/1‑146 (MQ=255)
                                                                                             gCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGAc                                                        <  1:787570/150‑1 (MQ=255)
                                                                                              cAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctcgatg                                                       >  1:839056/1‑143 (MQ=255)
                                                                                              cAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctcgatg                                                       >  2:402266/1‑143 (MQ=255)
                                                                                                 tGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCgca                                                                     >  1:133896/1‑133 (MQ=255)
                                                                                                  gAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCgcag                                                                    <  2:133896/133‑1 (MQ=255)
                                                                                                  gAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGcgac                                                   >  2:17371/1‑150 (MQ=255)
                                                                                                   aaGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGcgacc                                                  <  2:5709/150‑1 (MQ=255)
                                                                                                         tATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACaaa                                            >  2:425519/1‑150 (MQ=255)
                                                                                                             aaCCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTg                                         <  1:1034559/149‑1 (MQ=255)
                                                                                                              gcctcAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGc                                       >  2:904138/5‑150 (MQ=255)
                                                                                                               cctcAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGcc                                      <  1:277255/147‑2 (MQ=255)
                                                                                                                  cAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTa                                   <  2:313671/150‑1 (MQ=255)
                                                                                                                  cAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTa                                   <  2:706123/150‑1 (MQ=255)
                                                                                                                    ggATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                 <  2:962763/150‑1 (MQ=255)
                                                                                                                    ggATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                 >  2:594738/1‑150 (MQ=255)
                                                                                                                     gATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                 >  1:967438/1‑149 (MQ=255)
                                                                                                                      aTGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTgaaca                               >  2:1109305/1‑148 (MQ=255)
                                                                                                                      aTGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTgaaca                               >  2:738123/1‑148 (MQ=255)
                                                                                                                      aTGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAata                               <  1:979128/150‑1 (MQ=255)
                                                                                                                       tGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAacat                              <  1:610249/150‑1 (MQ=255)
                                                                                                                         ccACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTg                            >  1:965700/1‑150 (MQ=255)
                                                                                                                          caccacTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCgg                           >  2:442801/1‑150 (MQ=255)
                                                                                                                          cacaACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTgg                           >  2:582833/1‑150 (MQ=255)
                                                                                                                           accactaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGt                          <  2:182023/144‑1 (MQ=255)
                                                                                                                           acaACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGt                          >  2:248138/1‑150 (MQ=255)
                                                                                                                              actaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCtc                       >  1:443385/4‑149 (MQ=255)
                                                                                                                              actaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCtc                       >  2:950846/4‑149 (MQ=255)
                                                                                                                                   cccTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAAt                  <  2:1137026/150‑1 (MQ=255)
                                                                                                                                   cccTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAAt                  <  1:796092/150‑1 (MQ=255)
                                                                                                                                    ccTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATc                 >  1:479670/1‑150 (MQ=255)
                                                                                                                                        ctctcATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTg              <  2:678858/149‑1 (MQ=255)
                                                                                                                                             aTTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAGCAGCAGAGCAACGATCGCGCGCAGCAGGTCGTGGGCGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAAc          <  2:442265/148‑1 (MQ=255)
                                                                                                                                               tCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGt      >  2:123565/1‑150 (MQ=255)
                                                                                                                                                cTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGt      <  2:1191317/149‑1 (MQ=255)
                                                                                                                                                   aaCGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCg  >  2:303597/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCG  >  CP009273/8211‑8507

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: