Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,234 C→T 41.2% E27K (GAA→AAA)  citC ← [citrate [pro‑3S]‑lyase] ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273647,2340CT41.2% 207.6 / inf 287E27K (GAA→AAA) citC[citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base C (71/97);  new base T (53/65);  total (125/162)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 288 total aligned reads displayed.

TCAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGAT  >  CP009273/647087‑647383
                                                                                                                                                   |                                                                                                                                                     
tCAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTAt                                                                                                                                                     <  2:958639/150‑1 (MQ=255)
  aGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGc                                                                                                                                                   <  1:3652775/150‑1 (MQ=255)
         gTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAAtt                                                                                                                                            <  1:4944574/150‑1 (MQ=255)
           ccTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTgg                                                                                                                                          <  2:335146/150‑1 (MQ=255)
            cTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTggg                                                                                                                                         >  2:2243948/1‑150 (MQ=255)
              tCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCg                                                                                                                                       >  2:2736014/1‑150 (MQ=255)
               cACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTATCATGCAGGAATTGGGCGa                                                                                                                                      >  1:2441989/1‑150 (MQ=255)
                aCCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGAt                                                                                                                                     <  2:1672895/150‑1 (MQ=255)
                   gcgGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTc                                                                                                                                  <  2:1765887/150‑1 (MQ=255)
                    cggacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTcc                                                                                                                                 <  2:1315969/150‑1 (MQ=255)
                     ggacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCg                                                                                                                                <  1:3337388/150‑1 (MQ=255)
                        cggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCa                                                                                                                             <  1:623211/150‑1 (MQ=255)
                           aTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCAttt                                                                                                                          <  1:3507079/150‑1 (MQ=255)
                           aTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCAACGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAttt                                                                                                                          <  2:4916432/150‑1 (MQ=255)
                            ttCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtttt                                                                                                                         <  2:4972988/150‑1 (MQ=255)
                             tCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCAttttt                                                                                                                        >  2:4348285/1‑150 (MQ=255)
                               aCTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAttttttt                                                                                                                      <  2:404189/150‑1 (MQ=255)
                                cTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTa                                                                                                                     >  2:900416/1‑150 (MQ=255)
                                cTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTa                                                                                                                     >  1:4551304/1‑150 (MQ=255)
                                 tGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAt                                                                                                                    >  2:2343502/1‑150 (MQ=255)
                                  gATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtt                                                                                                                   >  2:1743047/1‑150 (MQ=255)
                                  gATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtt                                                                                                                   <  2:2389409/150‑1 (MQ=255)
                                  gATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtt                                                                                                                   <  1:5102425/150‑1 (MQ=255)
                                  gATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtt                                                                                                                   <  1:3561598/150‑1 (MQ=255)
                                   aTAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAttt                                                                                                                  >  1:2322152/1‑150 (MQ=255)
                                     aGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTc                                                                                                                >  2:147952/1‑150 (MQ=255)
                                         aCGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGTAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGaa                                                                                                            >  1:771893/1‑150 (MQ=255)
                                             aTTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                        <  2:4509683/150‑1 (MQ=255)
                                               ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtttt                                                                                                      <  2:2746939/150‑1 (MQ=255)
                                                    aataTTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTAcgcg                                                                                                 >  2:4864812/1‑150 (MQ=255)
                                                      tatTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCggg                                                                                               <  2:3770063/150‑1 (MQ=255)
                                                      tatTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCggg                                                                                               >  1:423051/1‑150 (MQ=255)
                                                            cGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGaaaa                                                                                         >  1:4420979/1‑150 (MQ=255)
                                                            cGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTGTTTTATTTTCTGAACGTTTTACGCGGGTGaaaa                                                                                         >  1:3110211/1‑150 (MQ=255)
                                                             ggCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAat                                                                                        >  2:1981259/1‑150 (MQ=255)
                                                             ggCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAat                                                                                        >  1:1610179/1‑150 (MQ=255)
                                                               cAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAatat                                                                                      >  2:2044184/1‑150 (MQ=255)
                                                               cAATTCCACAGCACGCGATAAGCTTTTCATCGCGGGTTACGGTTATAATTACTTCGAGTGTGGTGTCAACGCTCATATCATTTTTATGCAGGAATTTGTCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCTGGTGAAAatat                                                                                      <  2:4948293/150‑1 (MQ=255)
                                                                aaTTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATc                                                                                     >  2:1581883/1‑150 (MQ=255)
                                                                aaTTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATc                                                                                     >  2:2546230/1‑150 (MQ=255)
                                                                   tCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCAtt                                                                                  >  1:4464161/1‑150 (MQ=255)
                                                                   tCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCAtt                                                                                  <  1:261636/150‑1 (MQ=255)
                                                                    ccACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTg                                                                                 <  2:4241744/150‑1 (MQ=255)
                                                                    ccACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTg                                                                                 <  1:1514314/150‑1 (MQ=255)
                                                                    ccACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTg                                                                                 >  1:818688/1‑150 (MQ=255)
                                                                    ccACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTg                                                                                 >  1:2817608/1‑150 (MQ=255)
                                                                     cACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGc                                                                                >  1:386314/1‑150 (MQ=255)
                                                                     cACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGc                                                                                >  2:3545004/1‑150 (MQ=255)
                                                                      aCCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGcc                                                                               <  2:3607783/150‑1 (MQ=255)
                                                                      aCCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGcc                                                                               >  2:2267151/1‑150 (MQ=255)
                                                                       ccGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACTCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCg                                                                              <  2:2213169/150‑1 (MQ=255)
                                                                       ccGCACGCGATAAGCTTTTCATCGAGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCg                                                                              <  1:5050713/150‑1 (MQ=255)
                                                                             gcgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACata                                                                        <  1:4323607/150‑1 (MQ=255)
                                                                              cgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataa                                                                       <  1:3590710/150‑1 (MQ=255)
                                                                               gATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataat                                                                      >  1:4469074/1‑150 (MQ=255)
                                                                                aTAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataata                                                                     <  1:1413917/150‑1 (MQ=255)
                                                                                    gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaata                                                                 >  2:3512225/1‑150 (MQ=255)
                                                                                    gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaata                                                                 >  2:1075949/1‑150 (MQ=255)
                                                                                    gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGACTTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaata                                                                 >  2:355979/1‑150 (MQ=255)
                                                                                    gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaata                                                                 <  1:472165/150‑1 (MQ=255)
                                                                                      ttttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGt                                                               >  2:2879713/1‑150 (MQ=255)
                                                                                      ttttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGt                                                               >  1:1148964/1‑150 (MQ=255)
                                                                                       tttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTa                                                              >  1:2025085/1‑150 (MQ=255)
                                                                                       tttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTa                                                              >  1:2024965/1‑150 (MQ=255)
                                                                                           aTCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                          <  2:2834388/150‑1 (MQ=255)
                                                                                            tCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTg                                                         <  1:4840579/150‑1 (MQ=255)
                                                                                                  gTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTg                                                   <  1:1560771/150‑1 (MQ=255)
                                                                                                  gTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTg                                                   <  2:1903576/150‑1 (MQ=255)
                                                                                                    tACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCa                                                 <  2:3923174/150‑1 (MQ=255)
                                                                                                     aCGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAt                                                <  2:4227517/150‑1 (MQ=255)
                                                                                                     aCGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAt                                                <  1:2063943/150‑1 (MQ=255)
                                                                                                         taataaATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttg                                            >  2:4620247/1‑150 (MQ=255)
                                                                                                            taaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTa                                         <  1:4541006/150‑1 (MQ=255)
                                                                                                               ataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGGGCATGTTGTTATCg                                      >  1:2759463/1‑150 (MQ=255)
                                                                                                               ataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCg                                      >  2:3309105/1‑150 (MQ=255)
                                                                                                                taCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGa                                     >  1:3621943/1‑150 (MQ=255)
                                                                                                                 aCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAt                                    <  2:3705544/150‑1 (MQ=255)
                                                                                                                  cTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAtt                                   >  2:345200/1‑150 (MQ=255)
                                                                                                                   ttCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAttt                                  <  1:4801279/150‑1 (MQ=255)
                                                                                                                    tCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTg                                 <  1:212045/150‑1 (MQ=255)
                                                                                                                    tCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTg                                 <  2:4553326/150‑1 (MQ=255)
                                                                                                                    tCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTg                                 <  2:3218325/150‑1 (MQ=255)
                                                                                                                      gACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCa                               <  1:227939/150‑1 (MQ=255)
                                                                                                                       aCTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCaa                              <  2:539155/150‑1 (MQ=255)
                                                                                                                        cTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCTTTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAAc                             <  2:861388/150‑1 (MQ=255)
                                                                                                                           tggtgTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGaa                          >  1:1262232/1‑150 (MQ=255)
                                                                                                                              tgtCAACGCTCAAATCATTGTCATGCAGGAATTGGGCGCTTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgt                       <  2:3605214/150‑1 (MQ=255)
                                                                                                                               gtCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgtt                      <  1:3710092/150‑1 (MQ=255)
                                                                                                                                 cAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgttgt                    >  1:2962720/1‑150 (MQ=255)
                                                                                                                                   aCGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTc                  >  1:1891346/1‑150 (MQ=255)
                                                                                                                                    cGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCa                 <  1:2905079/150‑1 (MQ=255)
                                                                                                                                     gCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCaa                <  1:3651553/150‑1 (MQ=255)
                                                                                                                                       tCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATg              >  2:3737286/1‑150 (MQ=255)
                                                                                                                                         aaaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGtt            <  2:3821937/150‑1 (MQ=255)
                                                                                                                                         aaaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGtt            >  1:1079751/1‑150 (MQ=255)
                                                                                                                                         aaaTCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGtt            <  1:4170545/150‑1 (MQ=255)
                                                                                                                                          aaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTCCGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTg           <  2:4261494/150‑1 (MQ=255)
                                                                                                                                           aTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGc          >  2:691408/1‑150 (MQ=255)
                                                                                                                                              aTTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCaaa       <  2:4744071/150‑1 (MQ=255)
                                                                                                                                              aTTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCaaa       >  1:591553/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGAT  >  CP009273/647087‑647383

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: