Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,336,681 A→T 40.3% 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,6810AT40.3% 195.1 / inf 256Y155N (TAC→AAC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (80/73);  new base T (52/51);  total (132/124)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 256 total aligned reads displayed.

TCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGG  >  CP009273/4336534‑4336829
                                                                                                                                                   |                                                                                                                                                    
tCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCa                                                                                                                                                    <  1:2606380/150‑1 (MQ=255)
tCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACa                                                                                                                                                    <  2:269011/150‑1 (MQ=255)
 cagcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCAt                                                                                                                                                   <  1:1199842/150‑1 (MQ=255)
  agcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATg                                                                                                                                                  >  2:3611617/1‑150 (MQ=255)
  agcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATg                                                                                                                                                  >  1:171741/1‑150 (MQ=255)
     agGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTcc                                                                                                                                               <  2:3655434/150‑1 (MQ=255)
          tttGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAgcgc                                                                                                                                          <  1:3059671/150‑1 (MQ=255)
            tGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGcgccg                                                                                                                                        >  2:3350192/1‑150 (MQ=255)
             ggTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGcgccgc                                                                                                                                       >  1:2527681/1‑150 (MQ=255)
              gTTTCCTGGGAGAGAAACGTTTTGTTGGCAGAGCCTCCGCCTTTCGCAACGCAAAGGAATTTGTAATCAACGCCATCAACCGCGTAAAGGTCGATTTGCGAAGGCAGGGTAGTGCCTGTGTTGACCTATTTGGTCATGTCCAGCGCCGCa                                                                                                                                      >  2:1244851/1‑150 (MQ=255)
                 tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTc                                                                                                                                   >  1:2431369/1‑150 (MQ=255)
                  ccTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCt                                                                                                                                  >  1:538880/1‑150 (MQ=255)
                  ccTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCt                                                                                                                                  >  2:4065534/1‑150 (MQ=255)
                   cTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGTCGCATTCtg                                                                                                                                 >  2:4463/1‑150 (MQ=255)
                   cTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGTCGCATTCtg                                                                                                                                 >  2:4377/1‑150 (MQ=255)
                   cTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCtg                                                                                                                                 <  1:1520125/150‑1 (MQ=255)
                    tGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCtgt                                                                                                                                <  1:636131/150‑1 (MQ=255)
                     ggTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCtgtg                                                                                                                               <  2:1209548/150‑1 (MQ=255)
                     ggTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCtgtg                                                                                                                               <  2:1209502/150‑1 (MQ=255)
                        agagATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAAt                                                                                                                            >  2:500398/1‑150 (MQ=255)
                        agagATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAAt                                                                                                                            <  1:2198094/150‑1 (MQ=255)
                           gaTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAgc                                                                                                                         >  2:1960253/1‑150 (MQ=255)
                           gaTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgc                                                                                                                         >  2:703604/1‑150 (MQ=255)
                           gaTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTACCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgc                                                                                                                         >  1:2662859/1‑150 (MQ=255)
                            aTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgcg                                                                                                                        <  1:1989190/150‑1 (MQ=255)
                               cGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAg                                                                                                                     <  2:677524/150‑1 (MQ=255)
                                gTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAgg                                                                                                                    <  1:3997528/150‑1 (MQ=255)
                                  tttGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGtt                                                                                                                  <  1:1933867/150‑1 (MQ=255)
                                   ttgttgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTa                                                                                                                 <  1:296719/150‑1 (MQ=255)
                                     gttgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATc                                                                                                               >  2:3726384/1‑150 (MQ=255)
                                      ttgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCt                                                                                                              <  1:1046272/150‑1 (MQ=255)
                                         gCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCg                                                                                                           <  2:3769297/150‑1 (MQ=255)
                                           agagCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGat                                                                                                         >  1:185746/1‑150 (MQ=255)
                                                cgccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAgg                                                                                                    >  1:1947850/1‑150 (MQ=255)
                                                cgccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAgg                                                                                                    <  1:4038306/150‑1 (MQ=255)
                                                cgccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAgg                                                                                                    <  2:394332/150‑1 (MQ=255)
                                                 gccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgt                                                                                                   <  1:1678595/150‑1 (MQ=255)
                                                 gccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgt                                                                                                   <  1:1101834/150‑1 (MQ=255)
                                                   cgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgtgt                                                                                                 <  2:418923/150‑1 (MQ=255)
                                                    gccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGtt                                                                                                <  1:1103396/150‑1 (MQ=255)
                                                    gccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTGTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGtt                                                                                                >  1:3722698/1‑150 (MQ=255)
                                                     ccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTta                                                                                               >  1:2780116/1‑150 (MQ=255)
                                                      cTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtat                                                                                              <  2:2366548/150‑1 (MQ=255)
                                                       tttCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtata                                                                                             <  1:403501/150‑1 (MQ=255)
                                                       tttCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGCCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtata                                                                                             >  2:107869/1‑150 (MQ=255)
                                                        ttCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAg                                                                                            <  1:1120253/150‑1 (MQ=255)
                                                          cgcaacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGAc                                                                                          >  1:2709133/1‑150 (MQ=255)
                                                             aacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                                                                       >  2:950078/1‑150 (MQ=255)
                                                             aacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                                                                       >  1:1403575/1‑150 (MQ=255)
                                                             aacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                                                                       >  1:10291/1‑150 (MQ=255)
                                                             aacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTGTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                                                                       >  2:1664192/1‑150 (MQ=255)
                                                               cgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCtt                                                                                     >  2:898488/1‑150 (MQ=255)
                                                                   aaGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGa                                                                                 <  2:1535003/150‑1 (MQ=255)
                                                                    aGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGAc                                                                                >  1:2556779/1‑150 (MQ=255)
                                                                     ggAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCATGTTAGTGCCGGTGTTGACCTCTTTTTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACa                                                                               <  1:732505/150‑1 (MQ=255)
                                                                     ggAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATCGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACa                                                                               >  2:1338789/1‑150 (MQ=255)
                                                                      gAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAg                                                                              >  1:58064/1‑150 (MQ=255)
                                                                      gAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAg                                                                              >  2:274893/1‑150 (MQ=255)
                                                                       aaTTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGc                                                                             >  2:569067/1‑150 (MQ=255)
                                                                        aTTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGAAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGAGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCg                                                                            >  2:1872879/1‑150 (MQ=255)
                                                                         tttGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGt                                                                           >  1:3106332/1‑150 (MQ=255)
                                                                           tGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGttt                                                                         <  2:2714522/150‑1 (MQ=255)
                                                                            gTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttc                                                                        >  2:3596062/1‑150 (MQ=255)
                                                                             tACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttct                                                                       <  2:2203924/150‑1 (MQ=255)
                                                                              aCTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttctt                                                                      >  2:3570904/1‑150 (MQ=255)
                                                                               cTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttcttc                                                                     <  1:2348810/150‑1 (MQ=255)
                                                                                tCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtca                                                                    <  1:476422/150‑1 (MQ=255)
                                                                                      ccATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccg                                                              >  1:395711/1‑150 (MQ=255)
                                                                                        aTCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgcc                                                            <  1:346485/150‑1 (MQ=255)
                                                                                        aTCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgcc                                                            <  1:2871819/150‑1 (MQ=255)
                                                                                           tACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgccgcc                                                         >  1:2120452/1‑150 (MQ=255)
                                                                                           tACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgccgcc                                                         >  2:1871080/1‑150 (MQ=255)
                                                                                            aCCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTTACCTCTATGTACATGTCCAGCGACGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGAGATTCTTCATCAccgccgccg                                                        >  2:1148944/1‑150 (MQ=255)
                                                                                             ccGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCgg                                                       >  2:2801864/1‑150 (MQ=255)
                                                                                              cgcgTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGt                                                      <  1:3829794/150‑1 (MQ=255)
                                                                                                 gTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCca                                                   >  2:1265211/1‑150 (MQ=255)
                                                                                                 gTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCca                                                   >  1:2380133/1‑150 (MQ=255)
                                                                                                    cAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCacac                                                <  2:1348057/150‑1 (MQ=255)
                                                                                                    cAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCacac                                                <  2:1618856/150‑1 (MQ=255)
                                                                                                     aGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACAcg                                               <  1:2332821/150‑1 (MQ=255)
                                                                                                       gTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACAcgcg                                             >  2:2040954/1‑150 (MQ=255)
                                                                                                          gATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTg                                          >  1:1305619/1‑150 (MQ=255)
                                                                                                                cgcAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                                    >  2:3661470/1‑150 (MQ=255)
                                                                                                                   aggcaggTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtttttt                                 <  1:3649658/150‑1 (MQ=255)
                                                                                                                    ggcaggTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTa                                <  1:522868/150‑1 (MQ=255)
                                                                                                                       aggTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcg                             >  2:3164104/1‑150 (MQ=255)
                                                                                                                       aggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcg                             <  2:2523468/150‑1 (MQ=255)
                                                                                                                           tAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcgacga                         >  2:2862320/1‑150 (MQ=255)
                                                                                                                            aGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACgat                        <  2:3905458/150‑1 (MQ=255)
                                                                                                                              tGCCGGTGTTGACCTCTTTGTACATGACCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCAACACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACgatga                      >  1:807491/1‑150 (MQ=255)
                                                                                                                                   gtgtTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATcgcg                 >  1:4171990/1‑150 (MQ=255)
                                                                                                                                    tgtTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCgg                >  1:506810/1‑150 (MQ=255)
                                                                                                                                       tGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGc             <  1:2075811/150‑1 (MQ=255)
                                                                                                                                            tctTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTa        >  2:3526849/1‑150 (MQ=255)
                                                                                                                                             ctTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTAt       <  1:3541422/150‑1 (MQ=255)
                                                                                                                                              tttGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATc      <  1:3277056/150‑1 (MQ=255)
                                                                                                                                              tttGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATc      >  1:3362252/1‑150 (MQ=255)
                                                                                                                                               ttGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATcc     <  1:2671023/150‑1 (MQ=255)
                                                                                                                                                tGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCt    <  2:3109147/150‑1 (MQ=255)
                                                                                                                                                 gTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTg   <  1:3169150/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGG  >  CP009273/4336534‑4336829

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: