Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 800,779 G→T 100% E607D (GAG→GAT)  ybhJ → putative hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273800,7790GT95.9% 974.1 / 27.0 296E607D (GAG→GAT) ybhJputative hydratase
Reads supporting (aligned to +/- strand):  ref base G (3/9);  new base T (147/137);  total (150/146)
Fisher's exact test for biased strand distribution p-value = 8.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 299 total aligned reads displayed.

CGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTG  >  CP009273/800632‑800926
                                                                                                                                                   |                                                                                                                                                   
cGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgt                                                                                                                                                   <  2:3032261/150‑1 (MQ=255)
    gAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGttt                                                                                                                                               <  2:3762412/150‑1 (MQ=255)
     aaaCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGACTGTCAGCGAGCTGACAGATGTGTTTg                                                                                                                                              >  2:2572284/1‑150 (MQ=255)
       aCCTCTTCTTTTCGTACAAATCCGTTTGGTGTGGGGTAGTTTACCCTGTCTCGCCGCGATCCCGGTTTTGTTAGCTGATTTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgcg                                                                                                                                            <  1:1353234/150‑1 (MQ=255)
         ctctTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgcgcg                                                                                                                                          >  2:205942/1‑150 (MQ=255)
         ctctTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgcgcg                                                                                                                                          >  1:2800062/1‑150 (MQ=255)
         ctctTCTTATCGTTAAAATCGGTTTGGTGTGGCGGATTTTACCCGGTCGCGCCGCGATCCGGTTTATGTTAGAATAATTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgcgcg                                                                                                                                          <  1:927862/150‑1 (MQ=255)
           cttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCa                                                                                                                                        >  1:675964/1‑150 (MQ=255)
              cttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTa                                                                                                                                     >  1:4685928/1‑150 (MQ=255)
                  tCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCa                                                                                                                                 >  1:1035158/1‑150 (MQ=255)
                    gTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGa                                                                                                                               >  2:2571068/1‑150 (MQ=255)
                      tCAAATCCGTTTGGTGTGGGGGAGTTTACCCTGTCTCGCCGCGTTCCGGGTTTTGTTCGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGAtt                                                                                                                             <  1:3481243/150‑1 (MQ=255)
                         aaTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCt                                                                                                                          >  2:2538990/1‑150 (MQ=255)
                            ccgttgggTGTGGGGGTGTTTACCCTGTCTCGCCGCGATCCCGTTTTTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGt                                                                                                                       <  2:3790483/144‑1 (MQ=255)
                             cGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTc                                                                                                                      >  2:1702740/1‑150 (MQ=255)
                              gATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                                                                                     >  1:4078127/1‑150 (MQ=255)
                              gATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                                                                                     >  1:3632421/1‑150 (MQ=255)
                              gATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGGTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                                                                                     >  1:315319/1‑150 (MQ=255)
                               agtGGTCGGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGTTTATGTTAGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAg                                                                                                                    <  2:4603113/148‑1 (MQ=255)
                                       ggCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCatat                                                                                                            >  2:422121/1‑150 (MQ=255)
                                           ttgTTTACCCTGTCTCGCCGCGATCCCGGTTTTTTTAGCAGAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                                        <  1:3466125/148‑1 (MQ=255)
                                           gAGTTTACCCGGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                                        <  2:3905614/150‑1 (MQ=255)
                                            aGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                                                                                       >  2:1559601/1‑150 (MQ=255)
                                              tttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCg                                                                                                     >  1:789668/1‑150 (MQ=255)
                                              tttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCg                                                                                                     >  1:2134429/1‑150 (MQ=255)
                                               ttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGc                                                                                                    >  2:3133998/1‑150 (MQ=255)
                                                tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCt                                                                                                   >  2:691525/1‑150 (MQ=255)
                                                tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCt                                                                                                   >  1:3778715/1‑150 (MQ=255)
                                                 aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTg                                                                                                  >  1:2305884/1‑150 (MQ=255)
                                                 aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTg                                                                                                  >  1:3205298/1‑150 (MQ=255)
                                                   ccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCa                                                                                                >  1:1201860/1‑150 (MQ=255)
                                                    cTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCaa                                                                                               <  2:3247736/150‑1 (MQ=255)
                                                       tctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACt                                                                                            >  2:429774/1‑150 (MQ=255)
                                                        ctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTg                                                                                           >  1:3381230/1‑150 (MQ=255)
                                                         tcGCCGCGATCCCGGTTATGTTCGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGa                                                                                          <  2:495554/150‑1 (MQ=255)
                                                         tcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGa                                                                                          >  1:3950645/1‑150 (MQ=255)
                                                             cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                                      >  1:4334249/1‑150 (MQ=255)
                                                             cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                                      >  1:2457501/1‑150 (MQ=255)
                                                             cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                                      >  1:64011/1‑150 (MQ=255)
                                                                gATCCCGGTTATGTTCGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgc                                                                                   <  2:838535/150‑1 (MQ=255)
                                                                     ggtttttGTTAGCATAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                                                                              <  1:1409633/144‑1 (MQ=255)
                                                                     gggTTTTGTTAGCCTAATTACAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                                                                              <  1:2007569/149‑1 (MQ=255)
                                                                      ggTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATg                                                                             >  2:197445/1‑150 (MQ=255)
                                                                         tATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                                          >  2:3511781/1‑150 (MQ=255)
                                                                          aTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCt                                                                         >  1:1435715/1‑150 (MQ=255)
                                                                           tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa                                                                        >  2:1805347/1‑150 (MQ=255)
                                                                               tgCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCg                                                                    <  2:1648567/149‑1 (MQ=255)
                                                                               aGCAGAATTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCg                                                                    <  1:153160/150‑1 (MQ=255)
                                                                               aGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCg                                                                    <  2:4556372/150‑1 (MQ=255)
                                                                               aGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCg                                                                    <  2:329701/150‑1 (MQ=255)
                                                                                 cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGt                                                                  >  2:2978006/1‑150 (MQ=255)
                                                                                 cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGt                                                                  <  2:2140918/150‑1 (MQ=255)
                                                                                  aGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg                                                                 <  2:1044543/150‑1 (MQ=255)
                                                                                   gAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGTAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGa                                                                >  1:597046/1‑150 (MQ=255)
                                                                                       tAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAcc                                                            >  2:2806769/1‑150 (MQ=255)
                                                                                        aaaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                           >  1:1364303/1‑150 (MQ=255)
                                                                                             gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGa                                                      >  1:2858235/1‑150 (MQ=255)
                                                                                             gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGa                                                      >  1:4109653/1‑150 (MQ=255)
                                                                                               ctgctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATg                                                    <  1:682158/150‑1 (MQ=255)
                                                                                                  ctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGtt                                                 >  2:948761/1‑150 (MQ=255)
                                                                                                  ctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGtt                                                 <  1:1714989/150‑1 (MQ=255)
                                                                                                  ctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGtt                                                 >  1:4276410/1‑150 (MQ=255)
                                                                                                    gAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCt                                               <  1:1115193/150‑1 (MQ=255)
                                                                                                    gAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCt                                               <  1:2060583/150‑1 (MQ=255)
                                                                                                     aGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTg                                              >  2:4721553/1‑150 (MQ=255)
                                                                                                     aGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTg                                              >  2:4004844/1‑150 (MQ=255)
                                                                                                      gCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgc                                             <  1:1014757/150‑1 (MQ=255)
                                                                                                      gCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgc                                             <  1:2548287/150‑1 (MQ=255)
                                                                                                       cTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgcg                                            >  2:891745/1‑150 (MQ=255)
                                                                                                        tGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgcgc                                           >  1:2940702/1‑150 (MQ=255)
                                                                                                           aaaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTg                                        <  2:2259740/150‑1 (MQ=255)
                                                                                                           aaaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTg                                        <  1:422121/150‑1 (MQ=255)
                                                                                                            aaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGa                                       <  2:2800062/150‑1 (MQ=255)
                                                                                                            aaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGa                                       >  2:615793/1‑150 (MQ=255)
                                                                                                            aaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGa                                       >  1:2670725/1‑150 (MQ=255)
                                                                                                             aaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTACGCAGAGTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGaa                                      <  1:3900/150‑1 (MQ=255)
                                                                                                             aaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGaa                                      >  1:835004/1‑150 (MQ=255)
                                                                                                                cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAg                                   >  1:4529983/1‑150 (MQ=255)
                                                                                                                cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAg                                   >  1:1872679/1‑150 (MQ=255)
                                                                                                                cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAg                                   >  2:2700066/1‑150 (MQ=255)
                                                                                                                 aGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAgg                                  >  1:522609/1‑150 (MQ=255)
                                                                                                                    gTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAggcgg                               >  2:624575/1‑150 (MQ=255)
                                                                                                                          cGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCt                         >  2:1497644/1‑150 (MQ=255)
                                                                                                                           ggggTATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg                        <  1:1624380/150‑1 (MQ=255)
                                                                                                                              gAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCa                     <  1:3039334/150‑1 (MQ=255)
                                                                                                                               aaTGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAg                    >  1:4439007/1‑150 (MQ=255)
                                                                                                                                  gTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt                 >  2:2883347/1‑150 (MQ=255)
                                                                                                                                     aGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGtgtg              >  2:1614498/1‑150 (MQ=255)
                                                                                                                                      gCGAGCTGACAGATGTGTGTGCGCGCATTAAGCAGATTGCGGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGa             <  1:4545874/150‑1 (MQ=255)
                                                                                                                                        gAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAtt           <  2:3950645/150‑1 (MQ=255)
                                                                                                                                         aGCTGACAGATGTGTTTGCGCGCAGTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTg          <  2:1200582/150‑1 (MQ=255)
                                                                                                                                          gCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTgg         <  2:3988473/150‑1 (MQ=255)
                                                                                                                                          gCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTgg         <  1:4000514/150‑1 (MQ=255)
                                                                                                                                            tGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCg       <  2:789668/150‑1 (MQ=255)
                                                                                                                                             gACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg      >  2:4400615/1‑150 (MQ=255)
                                                                                                                                             gACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg      >  1:217407/1‑150 (MQ=255)
                                                                                                                                             gACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg      <  2:1287040/150‑1 (MQ=255)
                                                                                                                                              aCAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGt     >  2:3404935/1‑150 (MQ=255)
                                                                                                                                              aCAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGt     <  2:1843792/150‑1 (MQ=255)
                                                                                                                                                 gagGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTg  <  1:2571068/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTG  >  CP009273/800632‑800926

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: ‑