Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,336,681 A→T 33.6% Y155N (TAC→AAC)  fumB ← anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,336,6810AT33.6% 427.2 / inf 315Y155N (TAC→AAC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (97/112);  new base T (54/51);  total (152/163)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 319 total aligned reads displayed.

AGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGC  >  CP009273/4336532‑4336830
                                                                                                                                                     |                                                                                                                                                     
 gTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTAc                                                                                                                                                      >  1:3572520/1‑150 (MQ=255)
   cagcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCAt                                                                                                                                                    <  2:409143/150‑1 (MQ=255)
     gcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCACCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGt                                                                                                                                                  <  1:2791854/150‑1 (MQ=255)
      cagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTc                                                                                                                                                 >  2:3190090/1‑150 (MQ=255)
        gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCa                                                                                                                                               <  2:2903348/150‑1 (MQ=255)
        gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCa                                                                                                                                               <  2:428542/150‑1 (MQ=255)
          gCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAgc                                                                                                                                             <  1:3746842/150‑1 (MQ=255)
           cTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAgcg                                                                                                                                            <  2:1545477/150‑1 (MQ=255)
            tttGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAgcgc                                                                                                                                           >  1:1245850/1‑150 (MQ=255)
             ttGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGcgcc                                                                                                                                          >  2:824820/1‑150 (MQ=255)
             ttGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGcgcc                                                                                                                                          >  2:2356655/1‑150 (MQ=255)
             ttGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGcgcc                                                                                                                                          <  1:1370837/150‑1 (MQ=255)
              tGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGcgccg                                                                                                                                         <  2:4399902/150‑1 (MQ=255)
                gTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCa                                                                                                                                       >  1:1142716/1‑150 (MQ=255)
                  ttCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCAtt                                                                                                                                     <  1:3913041/150‑1 (MQ=255)
                    ccTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCt                                                                                                                                   <  1:4120476/150‑1 (MQ=255)
                     cTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCtg                                                                                                                                  >  2:2509047/1‑150 (MQ=255)
                      tGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCtgt                                                                                                                                 >  1:1740459/1‑150 (MQ=255)
                      tGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCtgt                                                                                                                                 >  2:4549318/1‑150 (MQ=255)
                        gTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGa                                                                                                                               <  2:1719010/150‑1 (MQ=255)
                        gTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGa                                                                                                                               <  2:256738/150‑1 (MQ=255)
                        gTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGa                                                                                                                               >  2:3625416/1‑150 (MQ=255)
                         tAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGaa                                                                                                                              <  2:190120/150‑1 (MQ=255)
                          agagATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAAt                                                                                                                             <  1:1799378/150‑1 (MQ=255)
                          agagATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAAt                                                                                                                             <  2:4362467/150‑1 (MQ=255)
                           gagaTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATa                                                                                                                            <  2:4766943/150‑1 (MQ=255)
                           gagaTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATa                                                                                                                            >  1:4203925/1‑150 (MQ=255)
                            agaTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAg                                                                                                                           <  2:3845005/150‑1 (MQ=255)
                             gaTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgc                                                                                                                          >  2:4628216/1‑150 (MQ=255)
                             gaTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgc                                                                                                                          <  1:2047608/150‑1 (MQ=255)
                                 cGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAg                                                                                                                      <  2:1207811/150‑1 (MQ=255)
                                  gTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAgg                                                                                                                     >  1:1317848/1‑150 (MQ=255)
                                   ttttGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGt                                                                                                                    <  1:3357887/150‑1 (MQ=255)
                                    tttGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGtt                                                                                                                   <  2:1037964/150‑1 (MQ=255)
                                      tgttgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTAt                                                                                                                 <  1:717122/150‑1 (MQ=255)
                                         tgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGTCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCAAGCGCCGCATTCTGTGAATAGCGCAGGTTATCtt                                                                                                              >  1:3047682/1‑150 (MQ=255)
                                         tgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCtt                                                                                                              >  2:2954006/1‑150 (MQ=255)
                                          ggCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTc                                                                                                             <  2:2451514/150‑1 (MQ=255)
                                          ggCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTc                                                                                                             >  2:1800336/1‑150 (MQ=255)
                                          ggCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTc                                                                                                             <  1:2512316/150‑1 (MQ=255)
                                           gCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCg                                                                                                            <  1:665601/150‑1 (MQ=255)
                                             agagCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGat                                                                                                          >  2:1136158/1‑150 (MQ=255)
                                               agCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat                                                                                                        >  1:4346501/1‑150 (MQ=255)
                                                        cTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtat                                                                                               <  1:4518339/150‑1 (MQ=255)
                                                         tttCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtata                                                                                              <  2:4224668/150‑1 (MQ=255)
                                                             gcaacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACg                                                                                          >  1:4778133/1‑150 (MQ=255)
                                                             gcaacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACg                                                                                          <  1:53914/150‑1 (MQ=255)
                                                               aacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                                                                        <  2:1094258/150‑1 (MQ=255)
                                                               aacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                                                                        >  2:795414/1‑150 (MQ=255)
                                                                 cgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCtt                                                                                      <  1:4627847/150‑1 (MQ=255)
                                                                 cgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCtt                                                                                      <  1:1973959/150‑1 (MQ=255)
                                                                 cgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCtt                                                                                      <  1:4225082/150‑1 (MQ=255)
                                                                  gcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCttt                                                                                     >  2:640037/1‑150 (MQ=255)
                                                                  gcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCttt                                                                                     <  2:35067/150‑1 (MQ=255)
                                                                   caaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTc                                                                                    >  2:4587512/1‑150 (MQ=255)
                                                                    aaaGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCg                                                                                   <  2:3454239/150‑1 (MQ=255)
                                                                    aaaGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCg                                                                                   >  1:3144621/1‑150 (MQ=255)
                                                                    aaaGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCg                                                                                   >  1:4017908/1‑150 (MQ=255)
                                                                      aGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGAc                                                                                 >  1:2289134/1‑150 (MQ=255)
                                                                      aGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGAc                                                                                 >  2:1982297/1‑150 (MQ=255)
                                                                       ggAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACa                                                                                >  2:3914627/1‑150 (MQ=255)
                                                                         aaTTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGc                                                                              <  1:1792051/150‑1 (MQ=255)
                                                                               tACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttct                                                                        <  1:1095787/150‑1 (MQ=255)
                                                                               tACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGATATAGACGCCTTTCGACAGCGTttct                                                                        >  2:387965/1‑150 (MQ=255)
                                                                                 cTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttcttc                                                                      <  1:2684249/150‑1 (MQ=255)
                                                                                  tCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtca                                                                     >  1:119412/1‑150 (MQ=255)
                                                                                   cATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtcat                                                                    >  1:1856198/1‑150 (MQ=255)
                                                                                          aTCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgcc                                                             >  1:1027338/1‑150 (MQ=255)
                                                                                             tACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgccgcc                                                          >  2:3548566/1‑150 (MQ=255)
                                                                                              aCCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgccgccg                                                         >  2:2882581/1‑150 (MQ=255)
                                                                                                  cgTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTcc                                                     >  2:4153552/1‑150 (MQ=255)
                                                                                                  cgTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATAGAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTcc                                                     >  2:1645073/1‑150 (MQ=255)
                                                                                                   gTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCca                                                    <  1:2903302/150‑1 (MQ=255)
                                                                                                     aCAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCcaca                                                  >  2:1388622/1‑150 (MQ=255)
                                                                                                      cAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCacac                                                 <  2:3302032/150‑1 (MQ=255)
                                                                                                        ggTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACAcgc                                               >  1:3985557/1‑150 (MQ=255)
                                                                                                            gATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTg                                           >  2:4478274/1‑150 (MQ=255)
                                                                                                                 gcgcAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtt                                      <  1:2394492/150‑1 (MQ=255)
                                                                                                                   gcaggcagGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtttt                                    <  2:2818937/150‑1 (MQ=255)
                                                                                                                   gcaggcagGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtttt                                    <  2:2979396/150‑1 (MQ=255)
                                                                                                                   gcaggcagGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtttt                                    <  1:3210108/150‑1 (MQ=255)
                                                                                                                      ggcaggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTa                                 >  2:1775912/1‑150 (MQ=255)
                                                                                                                       gcaggTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGGTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTAc                                >  1:4718897/1‑150 (MQ=255)
                                                                                                                         aggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcg                              <  1:824820/150‑1 (MQ=255)
                                                                                                                         aggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcg                              <  1:3535975/150‑1 (MQ=255)
                                                                                                                           gTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcgac                            >  1:594686/1‑150 (MQ=255)
                                                                                                                           gTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcgac                            <  2:2775671/150‑1 (MQ=255)
                                                                                                                           gTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcgac                            <  1:1300838/150‑1 (MQ=255)
                                                                                                                              aGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACgat                         <  1:2753863/150‑1 (MQ=255)
                                                                                                                               gTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACgatg                        <  2:485206/150‑1 (MQ=255)
                                                                                                                                  ccGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATc                     <  2:2144246/150‑1 (MQ=255)
                                                                                                                                     gtgtTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATcgcg                  <  1:1513353/150‑1 (MQ=255)
                                                                                                                                        ttGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTg               <  1:317578/150‑1 (MQ=255)
                                                                                                                                         tGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGc              >  2:3124753/1‑150 (MQ=255)
                                                                                                                                             ctctTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGt          <  1:1644297/150‑1 (MQ=255)
                                                                                                                                                 ttGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACAATGATCGCGGTGCCGGTATcc      <  2:225390/150‑1 (MQ=255)
                                                                                                                                                  tGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCt     >  1:4496217/1‑150 (MQ=255)
                                                                                                                                                  tGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCt     <  2:2026517/150‑1 (MQ=255)
                                                                                                                                                    tACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTgg   <  2:3235786/150‑1 (MQ=255)
                                                                                                                                                     aCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGc  >  2:2679283/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGC  >  CP009273/4336532‑4336830

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: