Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 188,863 T→A 22.4% D174E (GAT→GAA)  pyrH → uridylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273188,8630TA22.4% 743.3 / 223.9 313D174E (GAT→GAA) pyrHuridylate kinase
Reads supporting (aligned to +/- strand):  ref base T (114/129);  new base A (35/35);  total (149/164)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 315 total aligned reads displayed.

CTGTTGCGCAACAACCGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGG  >  CP009273/188714‑189011
                                                                                                                                                     |                                                                                                                                                    
 tGTTGCGCAACAACCGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATc                                                                                                                                                     >  1:4606382/1‑150 (MQ=255)
   ttGCGCAACAACCGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCg                                                                                                                                                   <  1:3951690/150‑1 (MQ=255)
       gcAACAACCGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGa                                                                                                                                               <  2:4624649/150‑1 (MQ=255)
       gcAACAACCGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGa                                                                                                                                               <  1:3754185/150‑1 (MQ=255)
         aacaacCGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGaaa                                                                                                                                             >  1:2006019/1‑150 (MQ=255)
            aacCGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGAt                                                                                                                                          >  1:4351636/1‑150 (MQ=255)
            aacCGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGAt                                                                                                                                          >  1:1966238/1‑150 (MQ=255)
              ccGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATcc                                                                                                                                        <  2:3054710/150‑1 (MQ=255)
              ccGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATcc                                                                                                                                        <  2:2629433/150‑1 (MQ=255)
              ccGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATcc                                                                                                                                        >  1:3473948/1‑150 (MQ=255)
                     gtgATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCAAACCGCa                                                                                                                                 >  1:57693/1‑150 (MQ=255)
                     gtgATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCa                                                                                                                                 >  1:360071/1‑150 (MQ=255)
                       gATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAAc                                                                                                                               <  2:524686/150‑1 (MQ=255)
                        aTCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAAcc                                                                                                                              >  2:3931058/1‑150 (MQ=255)
                           ctctCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATg                                                                                                                           >  1:4226660/1‑150 (MQ=255)
                           ctctCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATg                                                                                                                           >  2:3654662/1‑150 (MQ=255)
                           ctctCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATg                                                                                                                           >  2:1876296/1‑150 (MQ=255)
                            tctcCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGt                                                                                                                          >  1:4109265/1‑150 (MQ=255)
                             ctcCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTa                                                                                                                         <  2:1796495/150‑1 (MQ=255)
                              tcCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTAc                                                                                                                        >  1:2205808/1‑150 (MQ=255)
                                 gccgGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAg                                                                                                                     <  1:649127/150‑1 (MQ=255)
                                   cgGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCa                                                                                                                   <  2:2282313/150‑1 (MQ=255)
                                      tACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACt                                                                                                                <  1:3917458/150‑1 (MQ=255)
                                      tACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACt                                                                                                                <  1:1758252/150‑1 (MQ=255)
                                        cAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGa                                                                                                              >  2:1716007/1‑150 (MQ=255)
                                         aGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGAc                                                                                                             >  2:4317716/1‑150 (MQ=255)
                                           gTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACtt                                                                                                           >  1:3197102/1‑150 (MQ=255)
                                            tAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTa                                                                                                          <  1:2399993/150‑1 (MQ=255)
                                             aaCCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTAc                                                                                                         <  1:943821/150‑1 (MQ=255)
                                              aCCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACa                                                                                                        <  2:391565/150‑1 (MQ=255)
                                               cccGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAg                                                                                                       >  2:4129117/1‑150 (MQ=255)
                                               cccGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAg                                                                                                       <  2:2563229/150‑1 (MQ=255)
                                                ccGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGc                                                                                                      <  1:4249070/150‑1 (MQ=255)
                                                 cGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCg                                                                                                     <  1:4447691/150‑1 (MQ=255)
                                                 cGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCg                                                                                                     >  2:3353692/1‑150 (MQ=255)
                                                     cTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGt                                                                                                 <  2:1029854/150‑1 (MQ=255)
                                                      tttACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTg                                                                                                >  2:625713/1‑150 (MQ=255)
                                                       ttACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGc                                                                                               >  2:4719784/1‑150 (MQ=255)
                                                       ttACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGc                                                                                               <  1:1313572/150‑1 (MQ=255)
                                                        tACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCt                                                                                              >  1:2266724/1‑150 (MQ=255)
                                                          ccaccGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTgg                                                                                            <  2:1213754/150‑1 (MQ=255)
                                                              cGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGaaaa                                                                                        >  2:3451033/1‑150 (MQ=255)
                                                               gACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGaaaaa                                                                                       <  1:125172/150‑1 (MQ=255)
                                                                    agcagcTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCt                                                                                  <  1:2247581/150‑1 (MQ=255)
                                                                    agcagcTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCt                                                                                  <  1:1580553/150‑1 (MQ=255)
                                                                     gcagcTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTg                                                                                 <  2:4473484/150‑1 (MQ=255)
                                                                      cagcTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGa                                                                                <  1:3432447/150‑1 (MQ=255)
                                                                      cagcTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGa                                                                                >  1:2782761/1‑150 (MQ=255)
                                                                       agcTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGaa                                                                               >  2:4095329/1‑150 (MQ=255)
                                                                          ttGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGt                                                                            <  1:1190851/150‑1 (MQ=255)
                                                                           tGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTc                                                                           >  2:4352275/1‑150 (MQ=255)
                                                                            gCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCa                                                                          <  2:1910197/150‑1 (MQ=255)
                                                                            gCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCa                                                                          <  2:1749882/150‑1 (MQ=255)
                                                                            gCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCa                                                                          >  2:1949648/1‑150 (MQ=255)
                                                                               tGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATgg                                                                       <  1:2935538/150‑1 (MQ=255)
                                                                               tGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATgg                                                                       <  1:2995944/150‑1 (MQ=255)
                                                                               tGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATgg                                                                       <  1:4872655/150‑1 (MQ=255)
                                                                                 cGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGAc                                                                     >  2:2697839/1‑150 (MQ=255)
                                                                                 cGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGAc                                                                     >  2:4269278/1‑150 (MQ=255)
                                                                                        tCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTggcgg                                                              <  2:674616/150‑1 (MQ=255)
                                                                                            aaTTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCtt                                                          >  2:671866/1‑150 (MQ=255)
                                                                                              ttGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCa                                                        >  1:1422346/1‑150 (MQ=255)
                                                                                                gAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACg                                                      >  1:3251582/1‑150 (MQ=255)
                                                                                                 aaGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGc                                                     >  1:499331/1‑150 (MQ=255)
                                                                                                   gCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTg                                                   >  1:2723053/1‑150 (MQ=255)
                                                                                                   gCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTg                                                   >  1:3829547/1‑150 (MQ=255)
                                                                                                     cGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGc                                                 >  1:3671356/1‑150 (MQ=255)
                                                                                                     cGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGc                                                 >  2:2511537/1‑150 (MQ=255)
                                                                                                       aTGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTc                                               <  1:3090914/150‑1 (MQ=255)
                                                                                                       aTGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTc                                               <  1:3349658/150‑1 (MQ=255)
                                                                                                       aTGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTc                                               <  2:4196238/150‑1 (MQ=255)
                                                                                                       aTGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTc                                               <  1:1849086/150‑1 (MQ=255)
                                                                                                          tggtgCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTg                                            >  2:2387991/1‑150 (MQ=255)
                                                                                                           ggtgCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGa                                           >  1:3705401/1‑150 (MQ=255)
                                                                                                           ggtgCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGa                                           <  2:2202289/150‑1 (MQ=255)
                                                                                                           ggtgCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGa                                           >  1:105617/1‑150 (MQ=255)
                                                                                                            gtgCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGAc                                          <  2:3600446/150‑1 (MQ=255)
                                                                                                              gCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCa                                        >  1:2816054/1‑150 (MQ=255)
                                                                                                               cTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCAt                                       >  2:4760779/1‑150 (MQ=255)
                                                                                                               cTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCAt                                       >  1:964115/1‑150 (MQ=255)
                                                                                                               cTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCAt                                       >  1:964283/1‑150 (MQ=255)
                                                                                                               cTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCAt                                       <  2:2006019/150‑1 (MQ=255)
                                                                                                                     gCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTa                                 <  2:2887844/150‑1 (MQ=255)
                                                                                                                      cAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTAc                                <  2:3233781/150‑1 (MQ=255)
                                                                                                                      cAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTAc                                <  2:1383430/150‑1 (MQ=255)
                                                                                                                      cAACCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTAc                                <  1:2889911/150‑1 (MQ=255)
                                                                                                                       aaCCAAAGTTGACGGCGTGTTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTAcc                               >  1:92013/1‑150 (MQ=255)
                                                                                                                         ccAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGa                             >  2:2173117/1‑150 (MQ=255)
                                                                                                                              gTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCAGAAATTACCGATTCgt                        >  1:777471/1‑150 (MQ=255)
                                                                                                                                tGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCgtgt                      <  1:595577/150‑1 (MQ=255)
                                                                                                                                   cGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGtttt                   <  1:4664291/150‑1 (MQ=255)
                                                                                                                                       gtgtTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAat               <  1:2518128/150‑1 (MQ=255)
                                                                                                                                         gtTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATACGTGTTTTCAatat             >  1:3647536/1‑150 (MQ=255)
                                                                                                                                         gtTTACCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAatat             <  1:2656523/150‑1 (MQ=255)
                                                                                                                                            tACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGaa          <  2:1391582/150‑1 (MQ=255)
                                                                                                                                             aCCGCTGAACCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAAc         <  2:1596376/150‑1 (MQ=255)
                                                                                                                                                gCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACaaa      <  2:3278415/150‑1 (MQ=255)
                                                                                                                                                 cTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAAc     <  2:2760755/150‑1 (MQ=255)
                                                                                                                                                   gATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCg   >  2:2999331/1‑150 (MQ=255)
                                                                                                                                                    aTCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCgg  >  2:595392/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGTTGCGCAACAACCGTGTGGTGATCCTCTCCGCCGGTACAGGTAACCCGTTCTTTACCACCGACTCAGCAGCTTGCCTGCGTGGTATCGAAATTGAAGCCGATGTGGTGCTGAAAGCAACCAAAGTTGACGGCGTGTTTACCGCTGATCCGGCGAAAGATCCAACCGCAACCATGTACGAGCAACTGACTTACAGCGAAGTGCTGGAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGG  >  CP009273/188714‑189011

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: ‑