Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,185,155 G→T 25.5% I291I (ATC→ATA thiC ← phosphomethylpyrimidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,185,1550GT25.5% 561.4 / 152.5 270I291I (ATC→ATAthiCphosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base G (120/80);  new base T (35/34);  total (155/115)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 271 total aligned reads displayed.

ACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGAAGAGGTGACCGCCGA  >  CP009273/4185006‑4185301
                                                                                                                                                     |                                                                                                                                                  
    aGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCg                                                                                                                                                <  1:1253284/150‑1 (MQ=255)
    aGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCg                                                                                                                                                <  1:4510623/150‑1 (MQ=255)
     gcgcAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCgg                                                                                                                                               >  1:3093967/1‑150 (MQ=255)
      cgcAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGc                                                                                                                                              <  2:5110068/150‑1 (MQ=255)
      cgcAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGc                                                                                                                                              >  2:3710783/1‑150 (MQ=255)
      cgcAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGAGCGGCACTGTACCGATCGGc                                                                                                                                              >  2:3992639/1‑150 (MQ=255)
       gcagcaGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCa                                                                                                                                             <  2:5311706/150‑1 (MQ=255)
       gcagcaGTACGCCCGCATGGAAAGAGAAGTAATCCACACCTTGCACGGCCTGTTCCAGCTGAGTGTCGCGGAACGCTTCCAAGGTAAGATCTTCGGCGAACCCGTTAACCATATCCAGCGCCAGGTAGATCGGCACTGAACCTATCGGCa                                                                                                                                             >  2:2631364/1‑150 (MQ=255)
         agcagTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCAcc                                                                                                                                           >  1:4759540/1‑150 (MQ=255)
           cagTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCgg                                                                                                                                         >  2:1226776/1‑150 (MQ=255)
            agTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCggg                                                                                                                                        >  2:586940/1‑150 (MQ=255)
              tACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCt                                                                                                                                      <  1:2702419/150‑1 (MQ=255)
                 gCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGtt                                                                                                                                   <  2:1220736/150‑1 (MQ=255)
                  cccGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTa                                                                                                                                  <  1:3599007/150‑1 (MQ=255)
                   ccGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTAc                                                                                                                                 >  2:4185795/1‑150 (MQ=255)
                    cGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACg                                                                                                                                >  1:3258330/1‑150 (MQ=255)
                        tGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCaaa                                                                                                                            <  1:1673656/150‑1 (MQ=255)
                          gATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAAt                                                                                                                          >  2:69247/1‑150 (MQ=255)
                           aTAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATc                                                                                                                         >  2:1719205/1‑150 (MQ=255)
                            tAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATcc                                                                                                                        >  1:2749390/1‑150 (MQ=255)
                            tAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATcc                                                                                                                        <  2:3522201/150‑1 (MQ=255)
                             aGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCa                                                                                                                       <  2:2169130/150‑1 (MQ=255)
                              gTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATATTCGGAGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCAc                                                                                                                      >  2:1665977/1‑150 (MQ=255)
                               tGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACt                                                                                                                     <  2:1865653/150‑1 (MQ=255)
                                 aaGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTcg                                                                                                                   >  2:488546/1‑150 (MQ=255)
                                 aaGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTcg                                                                                                                   >  2:2835195/1‑150 (MQ=255)
                                   gTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTcgcg                                                                                                                 >  2:7695/1‑150 (MQ=255)
                                     aaTCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCggg                                                                                                               <  1:4235110/150‑1 (MQ=255)
                                     aaTCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCggg                                                                                                               >  2:4836685/1‑150 (MQ=255)
                                        ccACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGttt                                                                                                            <  1:2153595/150‑1 (MQ=255)
                                        ccACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGttt                                                                                                            >  2:2246647/1‑150 (MQ=255)
                                         cacaCCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTc                                                                                                           >  1:229376/1‑150 (MQ=255)
                                         cacaCCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTc                                                                                                           <  1:397579/150‑1 (MQ=255)
                                         cacaCCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTc                                                                                                           >  1:2933383/1‑150 (MQ=255)
                                         cacaCCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTc                                                                                                           >  1:4959019/1‑150 (MQ=255)
                                             ccTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGa                                                                                                       <  1:2106017/150‑1 (MQ=255)
                                                tGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAata                                                                                                    >  2:2036419/1‑150 (MQ=255)
                                                 gCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAatat                                                                                                   <  2:2603460/150‑1 (MQ=255)
                                                  cTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAAtata                                                                                                  <  1:5328689/150‑1 (MQ=255)
                                                   tCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAg                                                                                                 <  2:1495440/150‑1 (MQ=255)
                                                    cGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCNAAATCCACTCGCGGGTTTCGTGAATATAGc                                                                                                >  2:2523802/1‑150 (MQ=255)
                                                    cGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGc                                                                                                <  1:3675786/150‑1 (MQ=255)
                                                      gCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGa                                                                                              >  1:1886637/1‑150 (MQ=255)
                                                       ccTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGAc                                                                                             <  1:3386378/150‑1 (MQ=255)
                                                        cTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGAcc                                                                                            <  1:3535768/150‑1 (MQ=255)
                                                        cTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCCCTCGCGGGTTTCGTGAATATAGCGAcc                                                                                            <  1:909134/150‑1 (MQ=255)
                                                        cTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGAcc                                                                                            >  2:879443/1‑150 (MQ=255)
                                                        cTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGAcc                                                                                            <  1:1407941/150‑1 (MQ=255)
                                                         tGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCg                                                                                           >  2:3941823/1‑150 (MQ=255)
                                                                 cagcGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGgagag                                                                                   >  2:64665/1‑150 (MQ=255)
                                                                  agcGTGGCGCGGAACGCTTCCCAGGTAAGTTCTTCGGCGATCCCGTTACCCTTCTCCAGCGCCGGGTAGATCGGCACTGTACCGTTCGGCACCGGGCTGTTACGCAAAATCCACTCGGGGGTTTCGTGAATATAGCGACCGGTGgagagt                                                                                  <  2:2823164/150‑2 (MQ=255)
                                                                     gtgtCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATcc                                                                               <  2:4373127/150‑1 (MQ=255)
                                                                      tgtCGCGGAACGCTTCCCCGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCa                                                                              <  2:3365655/150‑1 (MQ=255)
                                                                          gcgGAACGCTTCCCTGTTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTATATCGGCACTGTACCGATCGGCACCGGGCTGTTACGAAAAATCCACTCGCGGGTTTCGTGAATATCGCGACCGGTGGAGAGATCCATCAc                                                                          <  2:1100494/150‑1 (MQ=255)
                                                                            ggAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCg                                                                        >  1:4321131/1‑150 (MQ=255)
                                                                            ggAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCg                                                                        >  2:3247626/1‑150 (MQ=255)
                                                                                cGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATc                                                                    >  1:818278/1‑150 (MQ=255)
                                                                                 gCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATcc                                                                   >  1:75457/1‑150 (MQ=255)
                                                                                 gCTTCCCAGGTAAGATCGTCGGCGATCCCGTTCCCCTTCTCCAGGGCCGGGTAGATCGGCACTGTACCTAGGGGCACGGGGCTGTTAGGCAATCTCCACTCGGGGGTTTCGTGATTTTAGCGACCGGTGGAGAGATCCATCACCGTATcc                                                                   <  2:2952100/150‑1 (MQ=255)
                                                                                      ccAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTcc                                                              >  2:4969741/1‑150 (MQ=255)
                                                                                        aGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTcccc                                                            <  2:5158844/150‑1 (MQ=255)
                                                                                         ggTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCa                                                           >  2:299520/1‑150 (MQ=255)
                                                                                            aaGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAgcg                                                        >  1:5222203/1‑150 (MQ=255)
                                                                                             aGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAgcgc                                                       >  2:3623325/1‑150 (MQ=255)
                                                                                                  ttCGGCGATCCCGTTAACCTACCCCCGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGa                                                  <  2:4980873/150‑1 (MQ=255)
                                                                                                    cGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGAcc                                                >  2:347793/1‑150 (MQ=255)
                                                                                                      gCGATCCCGTTACCCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCAt                                              <  1:204287/150‑1 (MQ=255)
                                                                                                        gATCCCTTTAACCTTCTCCAGCGCCTGGTCGATCGGCACTGTCCCGATCGGCACGGGGCTGTTACGCAAAAGCCACTCGCGGGTTTCGTGAATATCGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATAc                                            <  1:1987608/150‑1 (MQ=255)
                                                                                                        gATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATAc                                            >  1:2416361/1‑150 (MQ=255)
                                                                                                        gATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATAc                                            >  1:5125963/1‑150 (MQ=255)
                                                                                                        gATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATAc                                            >  1:158950/1‑150 (MQ=255)
                                                                                                           cccGTTAACCTTCTCCCGCGCCTGGTAGATCGGCACTGTACCGATCGGCACAGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGAGAATATCGCGACCGGTGGAGAGATCCATCACCGAATCCGCTCCCCAGCGCGTGGACCATACCAg                                         >  2:2344537/1‑150 (MQ=255)
                                                                                                           cccGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAg                                         <  2:4737139/150‑1 (MQ=255)
                                                                                                              gTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCtt                                      >  2:317477/1‑150 (MQ=255)
                                                                                                               ttAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCttt                                     >  1:1716016/1‑150 (MQ=255)
                                                                                                                 aaCCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTc                                   >  1:2810354/1‑150 (MQ=255)
                                                                                                                 aaCCTTCCCCAGCGCCTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTCCGAAAAATCCCCTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTc                                   <  1:3765252/150‑1 (MQ=255)
                                                                                                                  aCCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTcc                                  >  1:2896717/1‑150 (MQ=255)
                                                                                                                      tctcCCGCGCCTGGTTGATCGGCAGTGTACCGATCGGCACGGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCActt                              <  1:1076046/150‑1 (MQ=255)
                                                                                                                      tctcCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCActt                              >  2:3920426/1‑150 (MQ=255)
                                                                                                                       ctcCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCGGCTTTTCCActtc                             >  1:5115555/1‑150 (MQ=255)
                                                                                                                         ccAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCActtctt                           >  2:2082402/1‑150 (MQ=255)
                                                                                                                           aGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGTGCGTGGACCATACCAGCTTTTCCActtcttct                         >  1:2216861/1‑150 (MQ=255)
                                                                                                                           aGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCActtcttct                         >  2:5329831/1‑150 (MQ=255)
                                                                                                                            gcgcCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTACGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCActtcttctt                        <  2:3631659/150‑1 (MQ=255)
                                                                                                                               ccTGGTAGATCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGa                     >  2:2739542/1‑150 (MQ=255)
                                                                                                                                 tggaataTCGGAACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGCGTTTAGTGAATATAGCGACCGGTGGAGAGATCAAACACCGTATCCGCTCCACAGCGCGTGGACCAAACAAGCTTATCAACTTCTTCTTCGAtc                   >  2:5032413/7‑149 (MQ=255)
                                                                                                                                 tGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATg                   >  2:4041601/1‑150 (MQ=255)
                                                                                                                                 tGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATg                   >  1:411249/1‑150 (MQ=255)
                                                                                                                                  ggTTGATCGGCAGTTTACCGTTCGGCACGGGGTTTTTACGCAAAATCCACTGGGGGGTTTGGTTAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATgg                  <  1:2832123/150‑1 (MQ=255)
                                                                                                                                  ggTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATgg                  >  1:1172646/1‑150 (MQ=255)
                                                                                                                                   gTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGa                 >  2:205840/1‑150 (MQ=255)
                                                                                                                                    tAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGaa                >  1:97782/1‑150 (MQ=255)
                                                                                                                                    tAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGaa                >  2:2205446/1‑150 (MQ=255)
                                                                                                                                     aGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGAag               <  2:3084851/150‑1 (MQ=255)
                                                                                                                                       aTCGGCACTGTACCTATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGAagag             <  1:4876397/150‑1 (MQ=255)
                                                                                                                                       aTCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGAagag             >  1:800291/1‑150 (MQ=255)
                                                                                                                                           gCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGAAGAGGTGa         <  2:4544529/150‑1 (MQ=255)
                                                                                                                                                 tACCGATCGGCACGGGGCTGTTACGCAAAATCCCCTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGAAGAGGTGAccgccg   <  1:954902/150‑1 (MQ=255)
                                                                                                                                                  aCCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGAAGAGGTGACCGCCGa  <  1:3838949/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTAAGATCTTCGGCGATCCCGTTAACCTTCTCCAGCGCCTGGTAGATCGGCACTGTACCGATCGGCACCGGGCTGTTACGCAAAATCCACTCGCGGGTTTCGTGAATATAGCGACCGGTGGAGAGATCCATCACCGTATCCGCTCCCCAGCGCGTGGACCATACCAGCTTTTCCACTTCTTCTTCGATGGAAGAGGTGACCGCCGA  >  CP009273/4185006‑4185301

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: