Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,336,835 C→A 5.8% P103P (CCG→CCT 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,8350CA 5.8% 1151.5 / 29.2 304P103P (CCG→CCTfumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base C (154/132);  new base A (5/13);  total (159/145)
Fisher's exact test for biased strand distribution p-value = 4.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

TCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAG  >  CP009273/4336686‑4336982
                                                                                                                                                     |                                                                                                                                                   
 ccAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGTCAGGTCg                                                                                                                                                    >  2:3635676/1‑150 (MQ=255)
 ccAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCg                                                                                                                                                    <  2:4526335/150‑1 (MQ=255)
   aGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGgc                                                                                                                                                  >  2:1410197/1‑150 (MQ=255)
      gccgcATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGgcagc                                                                                                                                               >  2:858209/1‑150 (MQ=255)
             tCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCtt                                                                                                                                        >  1:3639181/1‑150 (MQ=255)
             tCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCtt                                                                                                                                        >  1:4339112/1‑150 (MQ=255)
             tCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCtt                                                                                                                                        <  1:1664220/150‑1 (MQ=255)
              cTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCttt                                                                                                                                       <  1:4539200/150‑1 (MQ=255)
              cTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCttt                                                                                                                                       <  1:2127538/150‑1 (MQ=255)
                 tgAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                                                                                                                                    >  1:922499/1‑150 (MQ=255)
                  gAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcg                                                                                                                                   >  2:1120412/1‑150 (MQ=255)
                  gAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcg                                                                                                                                   >  1:2793421/1‑150 (MQ=255)
                  gAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcg                                                                                                                                   >  1:2793905/1‑150 (MQ=255)
                   aaTAGCGCAGGTTATATTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcgg                                                                                                                                  >  2:2024094/1‑150 (MQ=255)
                      aGCGCAGGTTATCTTCGATAGAGGTGTTAGAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGa                                                                                                                               <  2:3961386/150‑1 (MQ=255)
                        cgcAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGAtt                                                                                                                             <  1:99242/150‑1 (MQ=255)
                        cgcAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGAtt                                                                                                                             <  1:2500702/150‑1 (MQ=255)
                        cgcAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATAGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGAtt                                                                                                                             <  1:3701596/150‑1 (MQ=255)
                         gcAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGAttt                                                                                                                            >  1:1540375/1‑150 (MQ=255)
                         gcAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGAttt                                                                                                                            >  1:3994624/1‑150 (MQ=255)
                          cAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTc                                                                                                                           >  2:878319/1‑150 (MQ=255)
                            ggTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCgg                                                                                                                         >  1:2368347/1‑150 (MQ=255)
                              ttATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAg                                                                                                                       >  2:2554430/1‑150 (MQ=255)
                               tATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACTCGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGTTTTCGGAGt                                                                                                                      <  1:2975637/150‑1 (MQ=255)
                                 tCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGttt                                                                                                                    <  2:3996658/150‑1 (MQ=255)
                                    tCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCtt                                                                                                                 <  2:3539678/150‑1 (MQ=255)
                                    tCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCtt                                                                                                                 <  2:1148132/150‑1 (MQ=255)
                                     cGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTa                                                                                                                >  2:1524942/1‑150 (MQ=255)
                                      gATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTaa                                                                                                               <  2:3414567/150‑1 (MQ=255)
                                       atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAg                                                                                                              <  1:2316370/150‑1 (MQ=255)
                                       atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAg                                                                                                              >  1:165229/1‑150 (MQ=255)
                                       atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAg                                                                                                              >  2:1892335/1‑150 (MQ=255)
                                         ataGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGaa                                                                                                            >  2:59621/1‑150 (MQ=255)
                                          taGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAAt                                                                                                           <  1:3260719/150‑1 (MQ=255)
                                           aGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAAtt                                                                                                          >  1:2060394/1‑150 (MQ=255)
                                           aGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAAtt                                                                                                          <  2:1378197/150‑1 (MQ=255)
                                             gtgtTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGc                                                                                                        >  2:956559/1‑150 (MQ=255)
                                               gtTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAg                                                                                                      <  1:625463/150‑1 (MQ=255)
                                                ttATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAgc                                                                                                     <  2:2375757/150‑1 (MQ=255)
                                                ttATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTAGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAgc                                                                                                     <  1:2432386/150‑1 (MQ=255)
                                                        gCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTAGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTa                                                                                             <  2:2336297/150‑1 (MQ=255)
                                                         ccTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTAc                                                                                            >  1:1426417/1‑150 (MQ=255)
                                                         ccTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTAc                                                                                            <  1:1541366/150‑1 (MQ=255)
                                                         ccTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTAc                                                                                            <  1:1683921/150‑1 (MQ=255)
                                                          cTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACt                                                                                           <  1:4123027/150‑1 (MQ=255)
                                                          cTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACt                                                                                           <  1:2754270/150‑1 (MQ=255)
                                                           tttCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACtt                                                                                          <  2:2773217/150‑1 (MQ=255)
                                                              cGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTc                                                                                       >  1:378777/1‑150 (MQ=255)
                                                                  aGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGttt                                                                                   <  1:2535678/150‑1 (MQ=255)
                                                                  aGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGttt                                                                                   <  2:2974042/150‑1 (MQ=255)
                                                                     gTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg                                                                                >  2:686685/1‑150 (MQ=255)
                                                                     gTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg                                                                                >  1:2708815/1‑150 (MQ=255)
                                                                      tttCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGc                                                                               >  1:3560769/1‑150 (MQ=255)
                                                                        tcttcATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTg                                                                             <  2:721669/150‑1 (MQ=255)
                                                                         cttcATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTgg                                                                            >  2:1533801/1‑150 (MQ=255)
                                                                          ttcATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGc                                                                           <  2:2060394/150‑1 (MQ=255)
                                                                          ttcATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTAGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGc                                                                           >  2:2859620/1‑150 (MQ=255)
                                                                            catcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCtt                                                                         >  1:1154967/1‑150 (MQ=255)
                                                                              tcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCt                                                                       >  1:1124423/1‑150 (MQ=255)
                                                                               caCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTg                                                                      <  1:1340059/150‑1 (MQ=255)
                                                                                    ccgccgGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCg                                                                 <  1:3012364/150‑1 (MQ=255)
                                                                                     cgccgGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGt                                                                >  1:995526/1‑150 (MQ=255)
                                                                                          gTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTgaaga                                                           >  1:1442055/1‑150 (MQ=255)
                                                                                             caccCGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTTAAGAATa                                                        <  2:2909882/150‑1 (MQ=255)
                                                                                             cacaCGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATa                                                        <  1:3773488/150‑1 (MQ=255)
                                                                                             cacaCGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGCAGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATa                                                        >  1:4312986/1‑150 (MQ=255)
                                                                                              acacGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAg                                                       >  2:149828/1‑150 (MQ=255)
                                                                                              acacGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAg                                                       >  1:1802833/1‑150 (MQ=255)
                                                                                                acGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGcc                                                     >  1:3599113/1‑150 (MQ=255)
                                                                                                   cgcTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCa                                                  >  2:2194772/1‑150 (MQ=255)
                                                                                                    gcTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCaa                                                 >  2:3039860/1‑150 (MQ=255)
                                                                                                       ggCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCt                                              <  1:1120341/150‑1 (MQ=255)
                                                                                                          cTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGtt                                           <  2:1787109/150‑1 (MQ=255)
                                                                                                             ttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCt                                        <  1:4395843/150‑1 (MQ=255)
                                                                                                             ttttAACGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCt                                        >  2:3216299/1‑150 (MQ=255)
                                                                                                              tttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTg                                       <  2:921076/150‑1 (MQ=255)
                                                                                                                tACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTAGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGgt                                     <  2:2754248/150‑1 (MQ=255)
                                                                                                                 aCCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGgtg                                    >  1:406293/1‑150 (MQ=255)
                                                                                                                    gacgaTGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGc                                 >  1:2409846/1‑150 (MQ=255)
                                                                                                                     acgaTGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGcc                                <  1:493524/150‑1 (MQ=255)
                                                                                                                     acgaTGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGcc                                <  1:858228/150‑1 (MQ=255)
                                                                                                                      cgaTGATCGCGGTGCCGGTATCCTGGCAGGTAGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCg                               <  2:2664135/150‑1 (MQ=255)
                                                                                                                       gatgatCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCgg                              >  2:1601550/1‑150 (MQ=255)
                                                                                                                        atgatCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCggg                             >  2:4227005/1‑150 (MQ=255)
                                                                                                                          gatCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCg                           >  2:2435571/1‑150 (MQ=255)
                                                                                                                           atCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCgg                          <  1:1840745/150‑1 (MQ=255)
                                                                                                                             cgcgGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAg                        <  2:172985/150‑1 (MQ=255)
                                                                                                                                 gTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATa                    >  1:3872748/1‑150 (MQ=255)
                                                                                                                                 gTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATa                    >  1:3872519/1‑150 (MQ=255)
                                                                                                                                    ccGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATaaaa                 <  1:2246347/150‑1 (MQ=255)
                                                                                                                                           gcTGGCAGGTCGGCAGCACGCCTTTGGCGGCGAGTTCGGAGTTTCTGAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGcgt          >  1:1039120/2‑150 (MQ=255)
                                                                                                                                             tGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGTTCGTGAAGAATAGCCGCAACCAGTGTCTGGTGTGCCGGGCGGAGCATAAAAGAAGcgtcg        >  1:4161632/1‑150 (MQ=255)
                                                                                                                                             tGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGcgtcg        >  1:3143678/1‑150 (MQ=255)
                                                                                                                                              ggCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGcgtcgt       >  2:621868/1‑150 (MQ=255)
                                                                                                                                                cAGGTCGGCAGCTCGCGTTTGGCGGCGATTTCGGAGTTTCTTCAGAAGTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGcgtcgtta     <  2:2737873/150‑3 (MQ=255)
                                                                                                                                                cAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGa     >  1:2625863/1‑150 (MQ=255)
                                                                                                                                                cAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGa     <  2:2738030/150‑1 (MQ=255)
                                                                                                                                                 aGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGaa    >  2:3773758/1‑150 (MQ=255)
                                                                                                                                                 aGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGaa    >  2:2615712/1‑150 (MQ=255)
                                                                                                                                                   gTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAg  <  2:1802793/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGTGTGCCGGGCGGAGCATAAAAGAAGCGTCGTGAAAG  >  CP009273/4336686‑4336982

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: