Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 110,883 A→T 8.7% Y318F (TAC→TTC)  guaC → GMP reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273110,8830AT 8.7% 735.4 / 48.6 231Y318F (TAC→TTC) guaCGMP reductase
Reads supporting (aligned to +/- strand):  ref base A (70/141);  new base T (11/9);  total (81/150)
Fisher's exact test for biased strand distribution p-value = 8.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 231 total aligned reads displayed.

ATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCG  >  CP009273/110735‑111028
                                                                                                                                                    |                                                                                                                                                 
aTGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATAc                                                                                                                                                  <  1:1109280/150‑1 (MQ=255)
    gCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTg                                                                                                                                              <  2:819358/150‑1 (MQ=255)
     cTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTgg                                                                                                                                             >  1:139628/1‑150 (MQ=255)
       ccGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTgggg                                                                                                                                           <  2:2455058/150‑1 (MQ=255)
          aGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCtt                                                                                                                                        >  2:1699493/1‑150 (MQ=255)
          aGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCtt                                                                                                                                        >  1:2857150/1‑150 (MQ=255)
              tGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACg                                                                                                                                    >  2:4153901/1‑150 (MQ=255)
              tGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACg                                                                                                                                    >  1:2212676/1‑150 (MQ=255)
              tGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACg                                                                                                                                    <  1:2624032/150‑1 (MQ=255)
                cGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGcc                                                                                                                                  >  1:830757/1‑150 (MQ=255)
                cGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGcc                                                                                                                                  <  1:2659547/150‑1 (MQ=255)
                  aTGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTg                                                                                                                                <  1:570527/150‑1 (MQ=255)
                      aaCGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAag                                                                                                                            >  2:3471542/1‑150 (MQ=255)
                      aaCGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAag                                                                                                                            <  2:1769213/150‑1 (MQ=255)
                          tCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCt                                                                                                                        <  2:1095866/150‑1 (MQ=255)
                          tCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCt                                                                                                                        <  1:3117373/150‑1 (MQ=255)
                          tCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCt                                                                                                                        <  2:3910534/150‑1 (MQ=255)
                          tCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCt                                                                                                                        <  1:868830/150‑1 (MQ=255)
                          tCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCt                                                                                                                        <  2:459388/150‑1 (MQ=255)
                           cACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTg                                                                                                                       <  2:1551759/150‑1 (MQ=255)
                            aCGTTGGCGGCGTTGCGGAATATCGAGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACAGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGa                                                                                                                      >  1:3157117/1‑150 (MQ=255)
                               ttGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTCAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCa                                                                                                                   <  1:3174350/150‑1 (MQ=255)
                               ttGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCa                                                                                                                   <  2:2684446/150‑1 (MQ=255)
                                 ggcggcGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAg                                                                                                                 >  2:3629530/1‑150 (MQ=255)
                                 ggcggcGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAg                                                                                                                 >  2:415775/1‑150 (MQ=255)
                                   cggcgTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATCCCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAgcg                                                                                                               >  1:76252/1‑150 (MQ=255)
                                     gcgTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGca                                                                                                             <  1:773572/150‑1 (MQ=255)
                                     gcgTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGca                                                                                                             <  2:559876/150‑1 (MQ=255)
                                      cgTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGcac                                                                                                            >  1:3679585/1‑150 (MQ=255)
                                      cgTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGcac                                                                                                            <  1:896212/150‑1 (MQ=255)
                                         tGCGGAATATCGCGCAGCAGAAGGTAACACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGcaccac                                                                                                         <  1:1567422/150‑1 (MQ=255)
                                         tGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGcaccac                                                                                                         <  2:457543/150‑1 (MQ=255)
                                            ggAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGtt                                                                                                      <  2:1335788/150‑1 (MQ=255)
                                            ggAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGtt                                                                                                      <  2:235328/150‑1 (MQ=255)
                                            ggAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGtt                                                                                                      <  1:28713/150‑1 (MQ=255)
                                               atatCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTAt                                                                                                   <  2:1663722/150‑1 (MQ=255)
                                               atatCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTAt                                                                                                   <  1:1733253/150‑1 (MQ=255)
                                               agatCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTAt                                                                                                   <  2:3863483/148‑1 (MQ=255)
                                                tatCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTAtt                                                                                                  <  1:1893370/150‑1 (MQ=255)
                                                 atCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTc                                                                                                 <  2:1569551/150‑1 (MQ=255)
                                                 atCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTc                                                                                                 <  1:1645276/150‑1 (MQ=255)
                                                  tCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCg                                                                                                <  1:3528895/150‑1 (MQ=255)
                                                    gcgcAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCgtg                                                                                              <  1:3130937/150‑1 (MQ=255)
                                                     cgcAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCgtgt                                                                                             <  1:487624/150‑1 (MQ=255)
                                                        agcagAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCa                                                                                          <  1:284717/150‑1 (MQ=255)
                                                              aGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACa                                                                                    <  2:1391311/150‑1 (MQ=255)
                                                              aGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACa                                                                                    >  2:3523896/1‑150 (MQ=255)
                                                               ggTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACaa                                                                                   >  1:240324/1‑150 (MQ=255)
                                                                  aaaaCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGaa                                                                                >  2:4174203/1‑150 (MQ=255)
                                                                    aaCCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGaaaa                                                                              <  1:636932/150‑1 (MQ=255)
                                                                       cGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCg                                                                           <  1:2003159/150‑1 (MQ=255)
                                                                        gTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGc                                                                          >  1:3047617/1‑150 (MQ=255)
                                                                          tAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCAt                                                                        >  1:955763/1‑150 (MQ=255)
                                                                          tAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCAt                                                                        <  1:3632398/150‑1 (MQ=255)
                                                                           aaGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATc                                                                       >  2:2919502/1‑150 (MQ=255)
                                                                           aaGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATc                                                                       <  1:430370/150‑1 (MQ=255)
                                                                              cTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTc                                                                    <  1:2696872/150‑1 (MQ=255)
                                                                               tGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTca                                                                   <  2:2482558/150‑1 (MQ=255)
                                                                                gCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTcaa                                                                  >  2:2568547/1‑150 (MQ=255)
                                                                                gCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTcaa                                                                  <  2:3548848/150‑1 (MQ=255)
                                                                                 ccGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTcaac                                                                 <  2:3154669/150‑1 (MQ=255)
                                                                                 ccGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTcaac                                                                 <  1:2888122/150‑1 (MQ=255)
                                                                                      gCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCt                                                            <  1:2785496/150‑1 (MQ=255)
                                                                                      gCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCt                                                            <  2:1390950/150‑1 (MQ=255)
                                                                                       cGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTg                                                           <  2:3144098/150‑1 (MQ=255)
                                                                                        gAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGt                                                          >  2:3703239/1‑150 (MQ=255)
                                                                                           gCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAAt                                                       <  1:205440/150‑1 (MQ=255)
                                                                                            cccGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAAtc                                                      <  1:2123792/150‑1 (MQ=255)
                                                                                                gTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTccc                                                  >  2:3993775/1‑150 (MQ=255)
                                                                                                 ttGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCa                                                 <  2:2973942/150‑1 (MQ=255)
                                                                                                  tGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCaa                                                >  1:2759565/1‑150 (MQ=255)
                                                                                                  tGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCCGTAATCTCCCaa                                                >  1:1431661/1‑150 (MQ=255)
                                                                                                        tACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTgg                                          <  2:2648319/150‑1 (MQ=255)
                                                                                                        tACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTgg                                          <  2:171747/150‑1 (MQ=255)
                                                                                                         aCCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGc                                         >  2:515687/1‑150 (MQ=255)
                                                                                                           cgcgcgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGt                                       <  2:1112898/150‑1 (MQ=255)
                                                                                                           cgcgcgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGt                                       <  2:3747798/150‑1 (MQ=255)
                                                                                                               cgAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAg                                   >  2:2285423/1‑150 (MQ=255)
                                                                                                                gagaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGc                                  <  2:1936070/150‑1 (MQ=255)
                                                                                                                    tatTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAaca                              <  2:3184564/150‑1 (MQ=255)
                                                                                                                     atTTTGGGCGGCCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAacac                             >  2:3124733/1‑150 (MQ=255)
                                                                                                                     atTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAacac                             <  2:3547023/150‑1 (MQ=255)
                                                                                                                        ttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGcc                          <  2:2677157/150‑1 (MQ=255)
                                                                                                                        ttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGcc                          >  2:1184948/1‑150 (MQ=255)
                                                                                                                        ttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGcc                          <  2:2781309/150‑1 (MQ=255)
                                                                                                                         tGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCa                         <  1:1993822/150‑1 (MQ=255)
                                                                                                                             cggcCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGt                     >  2:186342/1‑150 (MQ=255)
                                                                                                                               gcCTGCGTTCAGCTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTa                   >  1:3829569/1‑150 (MQ=255)
                                                                                                                                ccTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTat                  >  1:345836/1‑150 (MQ=255)
                                                                                                                                ccTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTat                  <  2:2175398/150‑1 (MQ=255)
                                                                                                                                 cTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatc                 >  2:2320935/1‑150 (MQ=255)
                                                                                                                                  tGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcc                >  1:2981667/1‑150 (MQ=255)
                                                                                                                                     gTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccat             >  2:1426215/1‑150 (MQ=255)
                                                                                                                                     gTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccat             <  1:3004772/150‑1 (MQ=255)
                                                                                                                                     gTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccat             <  1:1780956/150‑1 (MQ=255)
                                                                                                                                     gTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGAAGGAACAAGAAAACCGCTTCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGAGTTatcccat             <  2:4067746/150‑1 (MQ=255)
                                                                                                                                          gCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCAc        >  2:3939216/1‑150 (MQ=255)
                                                                                                                                           cTTGTACATTCGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACt       >  2:3046965/1‑150 (MQ=255)
                                                                                                                                           cTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACt       <  2:719887/150‑1 (MQ=255)
                                                                                                                                               tACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatc   <  1:404737/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATGAGCTCCGAGTCTGCGATGAAACGTCACGTTGGCGGCGTTGCGGAATATCGCGCAGCAGAAGGTAAAACCGTTAAGCTGCCGCTGCGAGGCCCGGTTGAAAATACCGCGCGAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCG  >  CP009273/110735‑111028

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