Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,422,841 C→A 62.3% D293Y (GAT→TAT)  purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,422,8410CA62.3% 164.6 / inf 236D293Y (GAT→TAT) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (52/37);  new base A (68/79);  total (120/116)
Fisher's exact test for biased strand distribution p-value = 8.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 242 total aligned reads displayed.

GCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCACATT  >  CP009273/2422694‑2422990
                                                                                                                                                   |                                                                                                                                                     
gctgctGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATAtt                                                                                                                                                     <  2:277450/150‑1 (MQ=255)
  tgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTcc                                                                                                                                                   <  2:2533985/150‑1 (MQ=255)
  tgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTcc                                                                                                                                                   <  1:1624349/150‑1 (MQ=255)
   gctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTccc                                                                                                                                                  >  2:2296355/1‑150 (MQ=255)
   gctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTccc                                                                                                                                                  <  2:677929/150‑1 (MQ=255)
    ctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCa                                                                                                                                                 <  1:441981/150‑1 (MQ=255)
     tgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCAt                                                                                                                                                >  2:4340415/1‑150 (MQ=255)
      ggCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCAtt                                                                                                                                               <  1:3818413/150‑1 (MQ=255)
           ggCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCgg                                                                                                                                          >  2:57918/1‑150 (MQ=255)
            gCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCggg                                                                                                                                         <  2:209509/150‑1 (MQ=255)
            gCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCggg                                                                                                                                         <  2:209852/150‑1 (MQ=255)
               tGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCaa                                                                                                                                      <  1:3261129/150‑1 (MQ=255)
                gATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAAt                                                                                                                                     <  1:2015438/150‑1 (MQ=255)
                 aTAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAAtt                                                                                                                                    >  2:4354588/1‑150 (MQ=255)
                    aaGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtc                                                                                                                                 >  2:2243371/1‑150 (MQ=255)
                       gTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCg                                                                                                                              <  2:3152007/150‑1 (MQ=255)
                       gTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCg                                                                                                                              <  2:721064/150‑1 (MQ=255)
                        tGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGc                                                                                                                             >  2:884264/1‑150 (MQ=255)
                           ggCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCa                                                                                                                          >  1:2498719/1‑150 (MQ=255)
                               aaCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGttt                                                                                                                      >  2:2474466/1‑150 (MQ=255)
                                       ggTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCa                                                                                                              >  2:2540133/1‑150 (MQ=255)
                                       ggTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCa                                                                                                              >  2:144988/1‑150 (MQ=255)
                                         tttttAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtc                                                                                                                                 <  1:4829271/129‑1 (MQ=255)
                                         tttttAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtc                                                                                                                                 >  2:4829271/1‑129 (MQ=255)
                                             tAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTca                                                                                                        <  1:2574162/150‑1 (MQ=255)
                                              aaCGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTcac                                                                                                       <  2:4241153/150‑1 (MQ=255)
                                              aaCGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGGTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTcac                                                                                                       <  2:4632303/150‑1 (MQ=255)
                                                cGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCacac                                                                                                     <  1:2490264/150‑1 (MQ=255)
                                                cGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCacac                                                                                                     <  1:2367990/150‑1 (MQ=255)
                                                   aGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcgcg                                                                                                  >  1:355256/1‑150 (MQ=255)
                                                       cTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCt                                                                                              <  1:1108249/150‑1 (MQ=255)
                                                        tGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTg                                                                                             <  1:3612143/150‑1 (MQ=255)
                                                        tGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTg                                                                                             <  1:1127366/150‑1 (MQ=255)
                                                        tGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTg                                                                                             <  2:3013014/150‑1 (MQ=255)
                                                        aggCGTTACGGTTTGCACAGAAAACGATCAATTTCCAGCGCGATATCACACGAGTGTTCTGGGATAGGGATCACCACGTCGATATCAAGATAGTCCAATTCGCGGGCAAATATCTCGACCAGATTCGAGCCCATATTCACCCGCGCGCTg                                                                                             >  2:2582339/2‑150 (MQ=255)
                                                         ggCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGt                                                                                            >  1:2558758/1‑150 (MQ=255)
                                                         ggCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGt                                                                                            <  1:178697/150‑1 (MQ=255)
                                                            ggTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTaaa                                                                                         <  2:2415077/150‑1 (MQ=255)
                                                             gTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAAc                                                                                        >  2:387621/1‑150 (MQ=255)
                                                             gTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAAc                                                                                        >  2:186220/1‑150 (MQ=255)
                                                               aCGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACgg                                                                                      >  2:2039459/1‑150 (MQ=255)
                                                                cGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGa                                                                                     >  1:466753/1‑150 (MQ=255)
                                                                        ccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGt                                                                             >  2:1780576/1‑150 (MQ=255)
                                                                         cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTc                                                                            <  1:1022573/150‑1 (MQ=255)
                                                                         cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTc                                                                            <  1:1305370/150‑1 (MQ=255)
                                                                         cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTc                                                                            <  1:1884028/150‑1 (MQ=255)
                                                                         cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTc                                                                            <  1:3140931/150‑1 (MQ=255)
                                                                         cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTc                                                                            <  1:3468547/150‑1 (MQ=255)
                                                                          aGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCg                                                                           >  2:3692108/1‑150 (MQ=255)
                                                                          aGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCg                                                                           <  1:2718826/150‑1 (MQ=255)
                                                                          aGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATACTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCg                                                                           >  1:3698195/1‑150 (MQ=255)
                                                                           gAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGa                                                                          <  2:4401665/150‑1 (MQ=255)
                                                                           gAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGa                                                                          <  2:1606111/150‑1 (MQ=255)
                                                                           gAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGa                                                                          >  1:4116966/1‑150 (MQ=255)
                                                                            aagaCGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGAt                                                                         <  2:3329835/147‑1 (MQ=255)
                                                                            aaTACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGAt                                                                         <  2:3329884/150‑1 (MQ=255)
                                                                             aTACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATa                                                                        <  1:4397432/150‑1 (MQ=255)
                                                                                cGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAAc                                                                     <  2:557938/150‑1 (MQ=255)
                                                                                  aGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGa                                                                   <  2:3566698/150‑1 (MQ=255)
                                                                                   gCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAg                                                                  <  2:2356810/150‑1 (MQ=255)
                                                                                     aaTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTc                                                                >  1:1462157/1‑150 (MQ=255)
                                                                                     aaTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTc                                                                >  1:675744/1‑150 (MQ=255)
                                                                                     aaTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTc                                                                <  1:208734/150‑1 (MQ=255)
                                                                                      aTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCATTTTTCTCGCCCAGTTTCGGGCCGATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGTTAAACGAGTcc                                                               <  2:490485/150‑1 (MQ=255)
                                                                                      aTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTcc                                                               <  1:2322295/150‑1 (MQ=255)
                                                                                         tCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcggg                                                            >  2:3538053/1‑150 (MQ=255)
                                                                                              cgcgATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCgggcgggc                                                       <  1:3398209/150‑1 (MQ=255)
                                                                                                 gATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCaaa                                                    >  1:981380/1‑150 (MQ=255)
                                                                                                  atatCACCCGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAg                                                   <  1:2060944/150‑1 (MQ=255)
                                                                                                    atCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTCAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGta                                                 >  1:4539884/1‑150 (MQ=255)
                                                                                                      cacaCGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtata                                               <  2:1641692/150‑1 (MQ=255)
                                                                                                      cacaCGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtata                                               >  1:2783058/1‑150 (MQ=255)
                                                                                                      cacaCGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtata                                               <  1:2884699/150‑1 (MQ=255)
                                                                                                       acacGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatac                                              <  1:1938199/150‑1 (MQ=255)
                                                                                                             ggTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTc                                        >  2:3196230/1‑150 (MQ=255)
                                                                                                              gTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCa                                       <  1:2296355/150‑1 (MQ=255)
                                                                                                                 tCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAAc                                    >  1:3479703/1‑150 (MQ=255)
                                                                                                                  cTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACa                                   <  2:4841058/150‑1 (MQ=255)
                                                                                                                    gggatcgggatcACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAgg                                 >  2:4865144/1‑150 (MQ=255)
                                                                                                                     ggatcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGc                                <  1:3856469/150‑1 (MQ=255)
                                                                                                                        tcgggatcACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACg                             <  1:3192340/150‑1 (MQ=255)
                                                                                                                          gggatcACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCCAAGTATACATACTCAAACAGGCACGGa                           >  2:1672558/1‑150 (MQ=255)
                                                                                                                          gggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGa                           >  1:786447/1‑150 (MQ=255)
                                                                                                                              tcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGc                       >  2:2203068/1‑150 (MQ=255)
                                                                                                                                accacGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTg                     <  2:2931745/150‑1 (MQ=255)
                                                                                                                                  cacGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGAc                   <  2:1437098/150‑1 (MQ=255)
                                                                                                                                   acGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGAcc                  >  2:364236/1‑150 (MQ=255)
                                                                                                                                     gTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCgg                >  2:1319356/1‑150 (MQ=255)
                                                                                                                                      tagATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGTTTGCTGACCGGa               <  1:2979199/148‑1 (MQ=255)
                                                                                                                                      tCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGa               >  2:1810155/1‑150 (MQ=255)
                                                                                                                                      tCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGa               >  1:2129158/1‑150 (MQ=255)
                                                                                                                                       cGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAt              <  2:4260786/150‑1 (MQ=255)
                                                                                                                                       cGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAt              >  1:1436834/1‑150 (MQ=255)
                                                                                                                                         atatCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTg            >  2:2311821/1‑150 (MQ=255)
                                                                                                                                            tCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCGGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCa         <  1:4537257/150‑1 (MQ=255)
                                                                                                                                             ccAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAg        >  1:522500/1‑150 (MQ=255)
                                                                                                                                             ccAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAg        >  1:522444/1‑150 (MQ=255)
                                                                                                                                             ccAGATATTCCCATTCGCGGGCAATGTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAg        <  1:4775497/150‑1 (MQ=255)
                                                                                                                                                   cTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCACAtt  >  2:4494518/1‑150 (MQ=255)
                                                                                                                                                   cTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCACAtt  >  2:1644701/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCACATT  >  CP009273/2422694‑2422990

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: