Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,336,681 A→T 14.5% 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,6810AT14.5% 914.8 / 130.9 310Y155N (TAC→AAC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (152/113);  new base T (20/25);  total (172/138)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 315 total aligned reads displayed.

AGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGC  >  CP009273/4336532‑4336830
                                                                                                                                                     |                                                                                                                                                     
 gTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTAc                                                                                                                                                      >  1:1000180/1‑150 (MQ=255)
 gTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTAc                                                                                                                                                      >  2:156792/1‑150 (MQ=255)
  tCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACa                                                                                                                                                     <  1:4127365/150‑1 (MQ=255)
  tCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACa                                                                                                                                                     >  1:734503/1‑150 (MQ=255)
  tCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACa                                                                                                                                                     >  1:2041965/1‑150 (MQ=255)
    agcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATg                                                                                                                                                   <  1:3923348/150‑1 (MQ=255)
     gcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGt                                                                                                                                                  <  2:1110685/150‑1 (MQ=255)
         ggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAg                                                                                                                                              >  1:2436035/1‑150 (MQ=255)
            tttGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAgcgc                                                                                                                                           >  2:1053315/1‑150 (MQ=255)
              tGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGcgccg                                                                                                                                         >  2:1868734/1‑150 (MQ=255)
                gTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCa                                                                                                                                       <  1:1446363/150‑1 (MQ=255)
                 tttCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCAt                                                                                                                                      <  2:4311548/150‑1 (MQ=255)
                  ttCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCAtt                                                                                                                                     >  2:3704227/1‑150 (MQ=255)
                   tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTc                                                                                                                                    <  2:24425/150‑1 (MQ=255)
                      tGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCtgt                                                                                                                                 >  2:3787397/1‑150 (MQ=255)
                       ggTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCtgtg                                                                                                                                <  2:3495623/150‑1 (MQ=255)
                        gTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGa                                                                                                                               >  2:2377132/1‑150 (MQ=255)
                        gTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCGACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGa                                                                                                                               <  2:127927/150‑1 (MQ=255)
                          agagATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAAt                                                                                                                             <  2:1506088/150‑1 (MQ=255)
                          agagATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAAt                                                                                                                             <  2:1028354/150‑1 (MQ=255)
                              aTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgcg                                                                                                                         >  1:556389/1‑150 (MQ=255)
                               tACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgcgc                                                                                                                        >  2:2790747/1‑150 (MQ=255)
                               tACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgcgc                                                                                                                        >  1:2872834/1‑150 (MQ=255)
                                  gTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAgg                                                                                                                     <  2:3106960/150‑1 (MQ=255)
                                  gTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAgg                                                                                                                     <  2:3836682/150‑1 (MQ=255)
                                     ttgttgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTa                                                                                                                  >  2:3397552/1‑150 (MQ=255)
                                           gCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCg                                                                                                            <  1:1001431/150‑1 (MQ=255)
                                            cAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGa                                                                                                           >  2:1594923/1‑150 (MQ=255)
                                            cAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGa                                                                                                           >  2:916927/1‑150 (MQ=255)
                                             agagCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGat                                                                                                          >  1:1657443/1‑150 (MQ=255)
                                               agCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat                                                                                                        <  2:71567/150‑1 (MQ=255)
                                                gccgccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGAtata                                                                                                       <  1:683222/150‑1 (MQ=255)
                                                  cgccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAgg                                                                                                     >  1:3073322/1‑150 (MQ=255)
                                                   gccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgt                                                                                                    >  2:3576919/1‑150 (MQ=255)
                                                    ccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgtg                                                                                                   >  1:802967/1‑150 (MQ=255)
                                                      gccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGtt                                                                                                 <  1:3692436/150‑1 (MQ=255)
                                                        cTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtat                                                                                               <  2:2924752/150‑1 (MQ=255)
                                                        cTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtat                                                                                               >  1:3277870/1‑150 (MQ=255)
                                                         tttCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtata                                                                                              <  1:158169/150‑1 (MQ=255)
                                                              caacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCAGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGc                                                                                         >  1:1313336/1‑150 (MQ=255)
                                                               aacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                                                                        >  2:3983201/1‑150 (MQ=255)
                                                               aacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                                                                        >  1:1267831/1‑150 (MQ=255)
                                                                acgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCt                                                                                       <  1:1553764/150‑1 (MQ=255)
                                                                acgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCt                                                                                       <  2:4232892/150‑1 (MQ=255)
                                                                  gcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCttt                                                                                     <  2:3885901/150‑1 (MQ=255)
                                                                  gcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCttt                                                                                     <  1:355879/150‑1 (MQ=255)
                                                                  gcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCttt                                                                                     <  2:3884898/150‑1 (MQ=255)
                                                                   caaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTc                                                                                    <  1:1891728/150‑1 (MQ=255)
                                                                    aaaGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCg                                                                                   >  1:1705027/1‑150 (MQ=255)
                                                                      aGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGAc                                                                                 <  2:1747774/150‑1 (MQ=255)
                                                                        gAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAg                                                                               <  2:886642/150‑1 (MQ=255)
                                                                         aaTTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGc                                                                              <  1:348137/150‑1 (MQ=255)
                                                                            ttGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGtt                                                                           <  2:2550422/150‑1 (MQ=255)
                                                                               tACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttct                                                                        >  2:3291574/1‑150 (MQ=255)
                                                                               tACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttct                                                                        >  2:2625296/1‑150 (MQ=255)
                                                                               tACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttct                                                                        >  2:2095492/1‑150 (MQ=255)
                                                                                aCTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttctt                                                                       >  2:135630/1‑150 (MQ=255)
                                                                                 cTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttcttc                                                                      >  1:2933252/1‑150 (MQ=255)
                                                                                 cTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTttcttc                                                                      >  2:2125743/1‑150 (MQ=255)
                                                                                  tCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATCTAGGTGTTATAGACGCCTCTCGACAGCGTTTCTtca                                                                     >  2:3478462/1‑150 (MQ=255)
                                                                                  tCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtca                                                                     <  1:3838167/150‑1 (MQ=255)
                                                                                   cATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtcat                                                                    <  1:3603772/150‑1 (MQ=255)
                                                                                   cATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtcat                                                                    >  2:2250612/1‑150 (MQ=255)
                                                                                    aTCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtcatc                                                                   >  2:3171525/1‑150 (MQ=255)
                                                                                    aTCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtcatc                                                                   >  1:890023/1‑150 (MQ=255)
                                                                                     tCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtcatca                                                                  <  1:2373253/150‑1 (MQ=255)
                                                                                       gCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAcc                                                                <  2:672279/150‑1 (MQ=255)
                                                                                       gCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAcc                                                                <  2:718941/150‑1 (MQ=255)
                                                                                        ccATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccg                                                               >  1:2903547/1‑150 (MQ=255)
                                                                                          aTCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgcc                                                             <  1:4326087/150‑1 (MQ=255)
                                                                                                  cgTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTcc                                                     <  2:3666498/150‑1 (MQ=255)
                                                                                                  cgTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTcc                                                     >  2:474636/1‑150 (MQ=255)
                                                                                                  cgTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTcc                                                     <  2:3019464/150‑1 (MQ=255)
                                                                                                   gTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCca                                                    >  2:4093500/1‑150 (MQ=255)
                                                                                                      cAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCacac                                                 <  1:31724/150‑1 (MQ=255)
                                                                                                      cAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCacac                                                 <  2:776784/150‑1 (MQ=255)
                                                                                                       aGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACAcg                                                <  1:2807915/150‑1 (MQ=255)
                                                                                                         gTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACAcgcg                                              <  2:1675221/150‑1 (MQ=255)
                                                                                                         gTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACAcgcg                                              <  1:3105165/150‑1 (MQ=255)
                                                                                                           cGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCt                                            <  1:3783113/150‑1 (MQ=255)
                                                                                                                 gcgcAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtt                                      <  1:1356959/150‑1 (MQ=255)
                                                                                                                   gcaggcagGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGCGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtttt                                    >  2:2182667/1‑150 (MQ=255)
                                                                                                                    caggcaggTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttttt                                   >  1:92987/1‑150 (MQ=255)
                                                                                                                    caggcaggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttttt                                   <  2:1657443/150‑1 (MQ=255)
                                                                                                                     aggcaggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtttttt                                  >  1:362663/1‑150 (MQ=255)
                                                                                                                        caggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGTTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTAcc                               >  1:2311184/1‑150 (MQ=255)
                                                                                                                        caggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTAcc                               >  1:3386536/1‑150 (MQ=255)
                                                                                                                         aggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcg                              >  2:4196571/1‑150 (MQ=255)
                                                                                                                         aggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcg                              >  1:4038344/1‑150 (MQ=255)
                                                                                                                              aGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACgat                         <  2:1000180/150‑1 (MQ=255)
                                                                                                                                     gtgtTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATcgcg                  <  1:1392120/150‑1 (MQ=255)
                                                                                                                                      tgtTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCgg                 >  1:3331438/1‑150 (MQ=255)
                                                                                                                                       gtTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGt                >  2:4165175/1‑150 (MQ=255)
                                                                                                                                         tGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGc              >  2:2678824/1‑150 (MQ=255)
                                                                                                                                          gACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGcc             >  1:2584903/1‑150 (MQ=255)
                                                                                                                                            ccTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCgg           >  1:2728797/1‑150 (MQ=255)
                                                                                                                                                  tGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCt     >  1:4152872/1‑150 (MQ=255)
                                                                                                                                                  tGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCt     >  2:3236326/1‑150 (MQ=255)
                                                                                                                                                    tACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTgg   >  2:1191385/1‑150 (MQ=255)
                                                                                                                                                     aCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGc  <  1:2125386/150‑1 (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: