Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 183,669 C→A 5.2% R490L (CGC→CTC)  glnD ← uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273183,6690CA 5.2% 986.4 / 20.2 250R490L (CGC→CTC) glnDuridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (117/120);  new base A (10/3);  total (127/123)
Fisher's exact test for biased strand distribution p-value = 8.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 271 total aligned reads displayed.

GTGGAGTTCGGCAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGG  >  CP009273/183521‑183818
                                                                                                                                                    |                                                                                                                                                     
gTGGAGTTCGGCAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGAtggcg                                                                                                                                                      <  2:3183409/150‑1 (MQ=255)
 tGGAGTTCGGCAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGAtggcgc                                                                                                                                                     >  2:2329639/1‑150 (MQ=255)
 tGGAGTTCGGCAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGAtggcgc                                                                                                                                                     >  1:4165090/1‑150 (MQ=255)
 tGGAGTTCGGCAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGAtggcgc                                                                                                                                                     >  1:1712469/1‑150 (MQ=255)
  ggAGTTCGGCAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGAtggcgct                                                                                                                                                    <  1:2657836/150‑1 (MQ=255)
   gAGTTCGGCAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGAtggcgctg                                                                                                                                                   >  1:1912446/1‑150 (MQ=255)
          gCAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGtt                                                                                                                                            <  2:3141974/150‑1 (MQ=255)
           cAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGttt                                                                                                                                           >  1:1796138/1‑150 (MQ=255)
            aaaaTGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTttc                                                                                                                                          <  2:185462/150‑1 (MQ=255)
             aaaTGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTttct                                                                                                                                         >  2:2724677/1‑150 (MQ=255)
              aaTGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTttctt                                                                                                                                        <  2:1213166/150‑1 (MQ=255)
                tGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGAGCTGGCGCGTTTCTTCa                                                                                                                                      <  1:1029041/150‑1 (MQ=255)
                 gCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCAc                                                                                                                                     >  1:2415507/1‑150 (MQ=255)
                  cACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACt                                                                                                                                    <  2:886974/150‑1 (MQ=255)
                     tACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGc                                                                                                                                 >  2:1331587/1‑150 (MQ=255)
                       cATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCaa                                                                                                                               >  1:201791/1‑150 (MQ=255)
                       cATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCaa                                                                                                                               >  1:1635684/1‑150 (MQ=255)
                       cATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCaa                                                                                                                               >  1:1104538/1‑150 (MQ=255)
                       cATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCaa                                                                                                                               >  1:1501739/1‑150 (MQ=255)
                        aTCCTGAGCACCGAGAATGTAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTTAAAAAGCGCGGCGATGAAAAAAAGATCAGTTGACGGCAGGCGCGGCCATACGTCCACACACAAAGGATGGCGCAGGCGCGTTTCTTCACTGGaaaa                                                                                                                              >  1:4066703/1‑150 (MQ=255)
                         tCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCaaaa                                                                                                                             <  1:2027549/150‑1 (MQ=255)
                               gCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTcc                                                                                                                       >  1:2532243/1‑150 (MQ=255)
                               gCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTcc                                                                                                                       >  1:2854034/1‑150 (MQ=255)
                                cACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCa                                                                                                                      >  2:1144765/1‑150 (MQ=255)
                                   cGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGtt                                                                                                                   <  1:3136808/150‑1 (MQ=255)
                                    gagaATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGttt                                                                                                                  >  2:2994119/1‑150 (MQ=255)
                                     agaATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTc                                                                                                                 >  1:4025881/1‑150 (MQ=255)
                                      gaATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCa                                                                                                                <  1:1834621/150‑1 (MQ=255)
                                        aTGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGc                                                                                                              >  1:3322872/1‑150 (MQ=255)
                                        aTGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGc                                                                                                              >  2:4155895/1‑150 (MQ=255)
                                        aTGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGc                                                                                                              <  2:2046058/150‑1 (MQ=255)
                                             gTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCAc                                                                                                         >  1:584709/1‑150 (MQ=255)
                                              tGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCAcg                                                                                                        <  1:2709548/150‑1 (MQ=255)
                                                gTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCAcgcg                                                                                                      <  1:2418220/150‑1 (MQ=255)
                                                        cgcgTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTAt                                                                                              <  1:2924653/150‑1 (MQ=255)
                                                        cgcgTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTAt                                                                                              <  1:2208594/150‑1 (MQ=255)
                                                        cgcgTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGAGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTAt                                                                                              >  2:1290190/1‑150 (MQ=255)
                                                          cgTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGt                                                                                            >  1:1067166/1‑150 (MQ=255)
                                                          cgTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGt                                                                                            >  1:1232002/1‑150 (MQ=255)
                                                          cgTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGt                                                                                            <  2:199948/150‑1 (MQ=255)
                                                           gTCCATTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGtt                                                                                           <  2:1237416/150‑1 (MQ=255)
                                                            tCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACACCGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTc                                                                                          <  2:1380380/150‑1 (MQ=255)
                                                               tttGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATc                                                                                       >  1:297717/1‑150 (MQ=255)
                                                                ttGGAGATATCGTGAAACAGCTCGGAGATGAAAATCAGCTCAGTTGAAGGCTGGCGCGGACAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTATCCAGTTTCAGCATCACGAGGATAGTATGTTCATcc                                                                                      <  2:354675/150‑1 (MQ=255)
                                                                  ggCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCAc                                                                                    >  1:1346422/1‑150 (MQ=255)
                                                                     gATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCgt                                                                                 >  2:3232167/1‑150 (MQ=255)
                                                                      atatCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCgtg                                                                                >  2:2890172/1‑150 (MQ=255)
                                                                       tatCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCgtgt                                                                               >  2:771214/1‑150 (MQ=255)
                                                                               aaCAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTgg                                                                       >  2:2349657/1‑150 (MQ=255)
                                                                               aaCAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTgg                                                                       <  2:3429575/150‑1 (MQ=255)
                                                                                aCAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGa                                                                      <  1:17768/150‑1 (MQ=255)
                                                                                      cggcgATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGAt                                                                <  2:3927367/150‑1 (MQ=255)
                                                                                        gcgATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCa                                                              <  2:941681/150‑1 (MQ=255)
                                                                                          gATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCaaa                                                            >  2:652799/1‑150 (MQ=255)
                                                                                           aTGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAAc                                                           <  2:2140616/150‑1 (MQ=255)
                                                                                           aTGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAAc                                                           <  2:958696/150‑1 (MQ=255)
                                                                                            tGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACt                                                          >  2:1483522/1‑150 (MQ=255)
                                                                                            tGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACt                                                          <  2:39664/150‑1 (MQ=255)
                                                                                            tGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCGGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACt                                                          <  2:2807979/150‑1 (MQ=255)
                                                                                              aaaaTCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTACTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGc                                                        >  1:3482608/1‑150 (MQ=255)
                                                                                                 aTCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATc                                                     <  2:2589698/150‑1 (MQ=255)
                                                                                                    aGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGc                                                  <  1:3987763/150‑1 (MQ=255)
                                                                                                    aGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGc                                                  <  2:3861643/150‑1 (MQ=255)
                                                                                                         aGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCcgac                                             >  1:3455646/1‑150 (MQ=255)
                                                                                                         aGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCcgac                                             >  2:2620718/1‑150 (MQ=255)
                                                                                                          gTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCcgacg                                            >  2:635303/1‑150 (MQ=255)
                                                                                                           ttGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCcgacga                                           <  2:2420706/150‑1 (MQ=255)
                                                                                                           ttGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCcgacga                                           >  1:3174010/1‑150 (MQ=255)
                                                                                                           ttGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCcgacga                                           >  1:3000027/1‑150 (MQ=255)
                                                                                                           ttGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCcgacga                                           <  2:379568/150‑1 (MQ=255)
                                                                                                            tGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGat                                          <  1:3603771/150‑1 (MQ=255)
                                                                                                            tGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGat                                          <  1:3708833/150‑1 (MQ=255)
                                                                                                            tGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGat                                          <  1:3973989/150‑1 (MQ=255)
                                                                                                            tGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGat                                          <  2:3999541/150‑1 (MQ=255)
                                                                                                              aCGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGatat                                        <  1:2794940/150‑1 (MQ=255)
                                                                                                              aCGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGatat                                        <  1:2460671/150‑1 (MQ=255)
                                                                                                               cGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATg                                       >  1:3462255/1‑150 (MQ=255)
                                                                                                                 gCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCATTTTCAGCATCAAGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCg                                     <  1:733232/150‑1 (MQ=255)
                                                                                                                 gCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCg                                     <  2:1029109/150‑1 (MQ=255)
                                                                                                                  cAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGa                                    >  2:2979705/1‑150 (MQ=255)
                                                                                                                  cAGGCGCGGCCAGACGTCCACACACAACGGATGGAGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGa                                    >  2:533042/1‑150 (MQ=255)
                                                                                                                   aGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGAc                                   >  1:3444756/1‑150 (MQ=255)
                                                                                                                    ggCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGAcc                                  <  2:2819384/150‑1 (MQ=255)
                                                                                                                    ggCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGAcc                                  <  1:905014/150‑1 (MQ=255)
                                                                                                                    ggCGCGGCCAGACGTCCACACACAACGGATGGAGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGAcc                                  >  2:2640598/1‑150 (MQ=255)
                                                                                                                      cgcgGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCAt                                >  2:3522564/1‑150 (MQ=255)
                                                                                                                      cgcgGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCAt                                >  2:3522606/1‑150 (MQ=255)
                                                                                                                       gcgGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCAtt                               <  2:432534/150‑1 (MQ=255)
                                                                                                                          gCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCg                            >  1:2205538/1‑150 (MQ=255)
                                                                                                                             aGACGTCCACACACATCGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCa                         <  2:1161835/150‑1 (MQ=255)
                                                                                                                              gACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCat                        >  1:1893564/1‑150 (MQ=255)
                                                                                                                                 gTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCAtata                     >  1:1741768/1‑150 (MQ=255)
                                                                                                                                 gTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCAtata                     >  1:2336472/1‑150 (MQ=255)
                                                                                                                                 gTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCAtata                     >  2:2282788/1‑150 (MQ=255)
                                                                                                                                  tCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATAc                    >  2:3669963/1‑150 (MQ=255)
                                                                                                                                  tCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATAc                    >  1:2120324/1‑150 (MQ=255)
                                                                                                                                      cacacaACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACgcgc                >  1:3917503/1‑150 (MQ=255)
                                                                                                                                        cacaACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCg              <  2:1705389/150‑1 (MQ=255)
                                                                                                                                            aCGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCa          <  1:3572727/150‑1 (MQ=255)
                                                                                                                                                  ggcgctggcgcGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTAt    >  2:3035752/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGAGTTCGGCAAAATGCACTACATCCTGAGCACCGAGAATGGAGTGGTCGCCGCCGCGTCCTTTGGCGATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGG  >  CP009273/183521‑183818

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: