Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,995,585 T→C 6.4% L100P (CTG→CCG) ‡ bkdAB → branched‑chain alpha‑keto acid dehydrogenase complex subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,995,5850TC 6.4% 529.6 / 16.6 170L100P (CTG→CCG) ‡bkdABbranched‑chain alpha‑keto acid dehydrogenase complex subunit beta
Reads supporting (aligned to +/- strand):  ref base T (83/76);  new base C (5/6);  total (88/82)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 170 total aligned reads displayed.

GGCAAATCGCGCGTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGAT  >  NC_002947/4995437‑4995732
                                                                                                                                                    |                                                                                                                                                   
 gCAAATCGCGCGTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTgg                                                                                                                                                   <  1:3603804/150‑1 (MQ=255)
 gCAAATCGCGCGTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTgg                                                                                                                                                   >  1:3088286/1‑150 (MQ=255)
 gCAAATCGCGCGTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTgg                                                                                                                                                   <  2:3329538/150‑1 (MQ=255)
   aaaTCGCGCGTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGcc                                                                                                                                                 >  1:4249204/1‑150 (MQ=255)
    aaTCGCGCGTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGccc                                                                                                                                                <  2:1674264/150‑1 (MQ=255)
       cgcgcgTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGcc                                                                                                                                             <  2:2198497/150‑1 (MQ=255)
       cgcgcgTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGcc                                                                                                                                             <  2:3969803/150‑1 (MQ=255)
        gcgcgTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCt                                                                                                                                            >  1:4101526/1‑150 (MQ=255)
          gcgGGTTCGACGCGCCCATCTCCGAGAGCGGCAGCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGCGCTCCCGGTCGCCGACGCCGGCTACCCGGCCGCCGACCAGGGCGGCGCCGAGCTGGCGCGCCGGc                                                                                                                                          >  2:3918147/1‑150 (MQ=255)
              gtTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTa                                                                                                                                      >  1:2017130/1‑150 (MQ=255)
                  gACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGCGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGt                                                                                                                                  >  1:890914/1‑150 (MQ=255)
                  gACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGt                                                                                                                                  >  1:3229579/1‑150 (MQ=255)
                     gcgcCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcg                                                                                                                               <  2:1217208/150‑1 (MQ=255)
                     gcgcCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcg                                                                                                                               >  1:4267346/1‑150 (MQ=255)
                       gcCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggc                                                                                                                             >  2:4192574/1‑150 (MQ=255)
                       gcCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggc                                                                                                                             >  2:1642412/1‑150 (MQ=255)
                          cATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccgg                                                                                                                          >  1:3683999/1‑150 (MQ=255)
                           aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                         >  1:1742710/1‑150 (MQ=255)
                           aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                         <  1:2759034/150‑1 (MQ=255)
                           aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                         >  1:973114/1‑150 (MQ=255)
                           aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                         <  2:2672572/150‑1 (MQ=255)
                           aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                         <  1:3242831/150‑1 (MQ=255)
                              tcCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAg                                                                                                                      >  1:2040095/1‑150 (MQ=255)
                              tcCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAg                                                                                                                      >  2:2948027/1‑150 (MQ=255)
                                cGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGtt                                                                                                                    <  1:1409203/150‑1 (MQ=255)
                                  agagCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCa                                                                                                                  <  2:2011210/150‑1 (MQ=255)
                                  agagCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCa                                                                                                                  <  2:3585321/150‑1 (MQ=255)
                                    agCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCAtt                                                                                                                <  1:3933201/150‑1 (MQ=255)
                                      cGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGc                                                                                                              <  1:4948026/150‑1 (MQ=255)
                                      cGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGc                                                                                                              <  1:966108/150‑1 (MQ=255)
                                        gCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGccc                                                                                                            >  2:3123017/1‑150 (MQ=255)
                                        gCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGccc                                                                                                            <  1:1706394/150‑1 (MQ=255)
                                         cATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGcccc                                                                                                           >  2:1237090/1‑150 (MQ=255)
                                          aTCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCg                                                                                                          >  1:2958879/1‑150 (MQ=255)
                                           tcgtcgGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGc                                                                                                         >  2:2741383/1‑150 (MQ=255)
                                             gtcgGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTg                                                                                                       >  2:2498929/1‑150 (MQ=255)
                                             gtcgGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGCCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTg                                                                                                       >  1:251984/1‑150 (MQ=255)
                                                    ccgccgTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTgc                                                                                                <  1:4970572/150‑1 (MQ=255)
                                                          tGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGc                                                                                          <  2:3088286/150‑1 (MQ=255)
                                                               aTGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTgc                                                                                     <  2:2443795/150‑1 (MQ=255)
                                                                 gggTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTgcgg                                                                                   >  2:214986/1‑150 (MQ=255)
                                                                       cTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCAGACTACTTCTACCCGGCATCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGcggcggcgg                                                                             >  1:3795084/1‑150 (MQ=255)
                                                                        tATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCggcggcggc                                                                            >  2:1939523/1‑150 (MQ=255)
                                                                        tATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCggcggcggc                                                                            >  2:1939577/1‑150 (MQ=255)
                                                                        tATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGACGACTACTTCTACCCGGCCTCCGACCAGATCGACTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCggcggcggc                                                                            >  2:4891548/1‑150 (MQ=255)
                                                                            gCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCt                                                                        >  1:253211/1‑150 (MQ=255)
                                                                             ccTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTa                                                                       <  1:2641026/150‑1 (MQ=255)
                                                                              cTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTAt                                                                      <  1:925484/150‑1 (MQ=255)
                                                                              cTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTAt                                                                      >  1:4548092/1‑150 (MQ=255)
                                                                                  gcCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcg                                                                  >  2:2092308/1‑150 (MQ=255)
                                                                                  gcCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcg                                                                  <  2:4106834/150‑1 (MQ=255)
                                                                                   cccGGTGGTGGAGATCCAGTTCGCCGACTAATTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcgg                                                                 <  1:4764443/150‑1 (MQ=255)
                                                                                    ccGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcggc                                                                <  2:4657977/150‑1 (MQ=255)
                                                                                    ccGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcggc                                                                >  1:3008840/1‑150 (MQ=255)
                                                                                      ggtggtGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCa                                                              <  2:738320/150‑1 (MQ=255)
                                                                                         ggtggAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGAc                                                           <  2:2784765/150‑1 (MQ=255)
                                                                                          gtggAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACt                                                          <  2:4500699/150‑1 (MQ=255)
                                                                                                gtCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGc                                                    <  1:843708/149‑1 (MQ=255)
                                                                                                 tgcAGTTTGCGGATTATGTCTATCCGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGcc                                                   >  2:1495558/3‑150 (MQ=255)
                                                                                                  gcAGTTTGCGGATTATGTCTATCCGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCa                                                  <  2:617276/149‑1 (MQ=255)
                                                                                                  ccAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCa                                                  <  2:4335693/150‑1 (MQ=255)
                                                                                                  ccAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCa                                                  <  1:3133067/150‑1 (MQ=255)
                                                                                                     gTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCagag                                               >  1:3550703/1‑150 (MQ=255)
                                                                                                           cGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGa                                         >  1:3919494/1‑150 (MQ=255)
                                                                                                            gACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGaa                                        <  1:214986/150‑1 (MQ=255)
                                                                                                             actactTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAg                                       <  1:4693182/150‑1 (MQ=255)
                                                                                                              ctactTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGc                                      <  2:4659136/150‑1 (MQ=255)
                                                                                                               tactTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCg                                     >  1:1014390/1‑150 (MQ=255)
                                                                                                               tactTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCg                                     >  1:339411/1‑150 (MQ=255)
                                                                                                                actTCTACCCTGCCTCCGAACAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCATTGCGGCGGCGGCATCTATGGCGGACAGACTCACAGCCAGAGCCCGGAAGCGa                                    <  1:3265287/150‑1 (MQ=255)
                                                                                                                actTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGa                                    <  1:4629170/150‑1 (MQ=255)
                                                                                                                 ctTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGAt                                   >  1:2045979/1‑150 (MQ=255)
                                                                                                                 ctTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGAt                                   <  1:2136662/150‑1 (MQ=255)
                                                                                                                 ctTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGAt                                   <  2:3919494/150‑1 (MQ=255)
                                                                                                                 ctTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGAt                                   <  2:3248705/150‑1 (MQ=255)
                                                                                                                  ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGACGCGATg                                  >  1:1780103/1‑150 (MQ=255)
                                                                                                                  ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATg                                  >  1:4252117/1‑150 (MQ=255)
                                                                                                                  ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATg                                  >  2:3727750/1‑150 (MQ=255)
                                                                                                                   tCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGt                                 <  2:1353909/150‑1 (MQ=255)
                                                                                                                    cTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGtt                                <  2:1648699/150‑1 (MQ=255)
                                                                                                                    cTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGtt                                <  2:2464567/150‑1 (MQ=255)
                                                                                                                     tatccGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTc                               >  2:3396033/4‑150 (MQ=255)
                                                                                                                      atccGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCa                              <  1:1016386/148‑1 (MQ=255)
                                                                                                                       cccGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCAc                             >  1:4670030/1‑150 (MQ=255)
                                                                                                                         cGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCAccc                           >  1:2017743/1‑150 (MQ=255)
                                                                                                                           gCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAg                         <  1:4293393/150‑1 (MQ=255)
                                                                                                                             cTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGTTGTTCACCCAGgt                       <  2:661444/150‑1 (MQ=255)
                                                                                                                             cTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGgt                       <  2:3567172/150‑1 (MQ=255)
                                                                                                                              tCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGAGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGgtg                      >  2:3072271/1‑150 (MQ=255)
                                                                                                                                cGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGtgtg                    <  1:280815/150‑1 (MQ=255)
                                                                                                                                cGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGtgtg                    >  1:776551/1‑150 (MQ=255)
                                                                                                                                cGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACACAGGtgtg                    >  2:1352952/1‑150 (MQ=255)
                                                                                                                                  aCCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCACCCATgtttcgg                  <  2:1172034/150‑6 (MQ=255)
                                                                                                                                  aCCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCACCCATgtttcgg                  <  2:5016162/150‑6 (MQ=255)
                                                                                                                                     aGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGcc               <  1:2377351/150‑1 (MQ=255)
                                                                                                                                             tcCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCAcc       >  2:1091820/1‑150 (MQ=255)
                                                                                                                                             tcCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCAcc       >  1:4434922/1‑150 (MQ=255)
                                                                                                                                             tcCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCAcc       >  1:599638/1‑150 (MQ=255)
                                                                                                                                                gAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCAGTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGCGCAGGCCGCACAcctgg    <  2:3431441/150‑4 (MQ=255)
                                                                                                                                                  gCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGAt  <  2:128869/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCAAATCGCGCGTGTTCGACGCGCCCATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGAT  >  NC_002947/4995437‑4995732

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: