Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,995,646 T→G 6.5% P120P (CCT→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,6460TG 6.5% 558.5 / 21.2 169P120P (CCT→CCGbkdABbranched‑chain alpha‑keto acid dehydrogenase complex subunit beta
Reads supporting (aligned to +/- strand):  ref base T (63/95);  new base G (5/6);  total (68/101)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 175 total aligned reads displayed.

ATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGAATGCGACGACCCGGT  >  NC_002947/4995500‑4995795
                                                                                                                                                  |                                                                                                                                                     
aTGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTgc                                                                                                                                                    <  2:2443795/150‑1 (MQ=255)
  gggTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTgcgg                                                                                                                                                  >  2:214986/1‑150 (MQ=255)
        cTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTAACCGGCATCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGcggcggcgg                                                                                                                                            >  2:4179382/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)
               cTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTAt                                                                                                                                     >  1:4548092/1‑150 (MQ=255)
               cTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTAt                                                                                                                                     <  2:2017130/150‑1 (MQ=255)
               cTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTAt                                                                                                                                     <  1:925484/150‑1 (MQ=255)
                   gcCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcg                                                                                                                                 <  2:4106834/150‑1 (MQ=255)
                   gcCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcg                                                                                                                                 >  1:3646295/1‑150 (MQ=255)
                    cccGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcgg                                                                                                                                >  1:3134409/1‑150 (MQ=255)
                    cccGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcgg                                                                                                                                <  1:4764364/150‑1 (MQ=255)
                     ccGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcggc                                                                                                                               >  1:3008840/1‑150 (MQ=255)
                     ccGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcggc                                                                                                                               >  1:3752820/1‑150 (MQ=255)
                       ggtggtGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCa                                                                                                                             <  2:738320/150‑1 (MQ=255)
                           gtggAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACt                                                                                                                         <  2:4500699/150‑1 (MQ=255)
                             ggAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTca                                                                                                                       <  1:1445347/150‑1 (MQ=255)
                                 gtCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGc                                                                                                                   <  1:3415590/149‑1 (MQ=255)
                                 gtCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGc                                                                                                                   <  1:843708/149‑1 (MQ=255)
                                  tgcAGTTTGCGGATTATGTCTATCCGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGcc                                                                                                                  >  1:3957254/3‑150 (MQ=255)
                                   gcAGTTTGCGGATTATGTCTATCCGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCa                                                                                                                 <  2:617276/149‑1 (MQ=255)
                                   ccAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCa                                                                                                                 <  2:4335693/150‑1 (MQ=255)
                                    cAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCag                                                                                                                <  1:3379097/150‑1 (MQ=255)
                                      gTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCagag                                                                                                              >  1:2462877/1‑150 (MQ=255)
                                              actactTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAg                                                                                                      <  1:723961/150‑1 (MQ=255)
                                              actactTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAg                                                                                                      >  2:890901/1‑150 (MQ=255)
                                               ctactTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGc                                                                                                     <  2:4659136/150‑1 (MQ=255)
                                                tactTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCg                                                                                                    >  1:1014390/1‑150 (MQ=255)
                                                 actTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGa                                                                                                   <  1:1766764/150‑1 (MQ=255)
                                                 actTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGa                                                                                                   <  1:4629170/150‑1 (MQ=255)
                                                  ctTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGAt                                                                                                  <  2:3248886/150‑1 (MQ=255)
                                                  ctTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGAt                                                                                                  <  1:2136662/150‑1 (MQ=255)
                                                   ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGACGCGATg                                                                                                 >  1:1780103/1‑150 (MQ=255)
                                                   ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATg                                                                                                 >  2:3422007/1‑150 (MQ=255)
                                                   ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATg                                                                                                 >  2:3727750/1‑150 (MQ=255)
                                                   ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATg                                                                                                 >  1:4252117/1‑150 (MQ=255)
                                                   ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATg                                                                                                 <  1:4262922/150‑1 (MQ=255)
                                                    tCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGt                                                                                                <  2:253211/150‑1 (MQ=255)
                                                    tCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGt                                                                                                <  2:1353789/150‑1 (MQ=255)
                                                     cTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGtt                                                                                               <  2:1648699/150‑1 (MQ=255)
                                                       atccGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCa                                                                                             <  1:1016386/148‑1 (MQ=255)
                                                        cccGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCAc                                                                                            >  1:4670030/1‑150 (MQ=255)
                                                          cGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCAccc                                                                                          >  1:2019034/1‑150 (MQ=255)
                                                            gcctatgACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCACCCat                                                                                        <  2:2700475/144‑2 (MQ=255)
                                                              cTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGTTGTTCACCCAGgt                                                                                      <  2:661444/150‑1 (MQ=255)
                                                              cTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGgt                                                                                      <  2:3567172/150‑1 (MQ=255)
                                                                   aCCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCACCCATgtttcgg                                                                                 <  2:5016162/150‑6 (MQ=255)
                                                                   aCCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCACCCATgtttcgg                                                                                 <  2:1172034/150‑6 (MQ=255)
                                                                       gATCGTCTCCGAGCTGGCCCGCCTGCGTGACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCt                                                                             <  2:3434783/150‑1 (MQ=255)
                                                                              tcCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCAcc                                                                      >  1:4434922/1‑150 (MQ=255)
                                                                              tcCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCAcc                                                                      >  1:2811567/1‑150 (MQ=255)
                                                                                   gCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGAt                                                                 <  2:128869/150‑1 (MQ=255)
                                                                                   gCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGAt                                                                 <  1:3875347/150‑1 (MQ=255)
                                                                                   gCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCCGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGAt                                                                 >  2:1479400/1‑150 (MQ=255)
                                                                                      ggCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGcc                                                              <  2:2040264/150‑1 (MQ=255)
                                                                                         ccGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTc                                                           >  2:1373562/1‑150 (MQ=255)
                                                                                         acGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTTACACGATGACCCTGCGCATGCCTTGAGGAGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGTTGTGCGGCCTGCGCACCGTGATGCCGTc                                                           >  2:3952258/2‑150 (MQ=255)
                                                                                          cGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTcc                                                          <  2:1101995/150‑1 (MQ=255)
                                                                                              tGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAAcc                                                      <  2:4249204/150‑1 (MQ=255)
                                                                                               gCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAAccc                                                     <  1:1993345/150‑1 (MQ=255)
                                                                                                 gTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCtt                                                   <  1:3123017/150‑1 (MQ=255)
                                                                                                 gTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCtt                                                   >  1:728447/1‑150 (MQ=255)
                                                                                                  ttACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACTGCCAGAGCGGGGAAGCGATGTTCACCGCGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTa                                                  <  2:251984/150‑1 (MQ=255)
                                                                                                  ttACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCTGTCCAACCCTTa                                                  >  2:3830602/1‑150 (MQ=255)
                                                                                                  ttACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTa                                                  >  1:32199/1‑150 (MQ=255)
                                                                                                   tACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTAt                                                 <  2:4267346/150‑1 (MQ=255)
                                                                                                   tACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTAt                                                 <  1:1284828/150‑1 (MQ=255)
                                                                                                   tACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTAt                                                 <  1:4179382/150‑1 (MQ=255)
                                                                                                   tACCGTTCGGCCGGCGAGTTCATTGCCACGCTGACCCTGATAATGCCTTGCGGCGGCTGCAACAATGGCGTCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACAAAGGTGTGCGGACAGCTCACCGAGAGGCCGACCAAACCATaa                                                 >  1:1395863/1‑149 (MQ=255)
                                                                                                      cGTTCGGCCGGCGAGTTCGTTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGAc                                              <  2:1742710/150‑1 (MQ=255)
                                                                                                            gcctgcGATTTCTTTGACCCTCTGACACTGCTCATGACTTGGGGCGGCGGCTTCTATGGCGGCCAGTCTATCAGCCAGCGCGCGGAAGCGTTGTTCACGCATGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCaaa                                        <  1:117516/150‑1 (MQ=255)
                                                                                                               ggcGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGc                                     >  1:2694364/1‑150 (MQ=255)
                                                                                                                 cGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCt                                   <  2:728447/150‑1 (MQ=255)
                                                                                                                  gAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACATCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTg                                  <  1:4891548/150‑1 (MQ=255)
                                                                                                                  gAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTg                                  <  2:3229579/150‑1 (MQ=255)
                                                                                                                   aGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGt                                 >  2:2377500/1‑150 (MQ=255)
                                                                                                                   aGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGt                                 <  1:4396490/150‑1 (MQ=255)
                                                                                                                   aGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGt                                 <  1:661287/150‑1 (MQ=255)
                                                                                                                   aGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGt                                 <  1:5043397/150‑1 (MQ=255)
                                                                                                                       cATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGAt                             <  1:1642412/150‑1 (MQ=255)
                                                                                                                       cATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGAt                             <  2:4101526/150‑1 (MQ=255)
                                                                                                                           gCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGcc                         >  1:347107/1‑150 (MQ=255)
                                                                                                                           gCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGcc                         <  1:2191798/150‑1 (MQ=255)
                                                                                                                            ccccGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCt                        <  1:2092308/150‑1 (MQ=255)
                                                                                                                              ccGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCtcg                      >  2:935203/1‑150 (MQ=255)
                                                                                                                               cGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCtcga                     >  2:2564255/1‑150 (MQ=255)
                                                                                                                               cGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCtcga                     >  2:2564247/1‑150 (MQ=255)
                                                                                                                                    aCCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGaa                <  1:2741383/150‑1 (MQ=255)
                                                                                                                                    aCCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGaa                >  2:2870958/1‑150 (MQ=255)
                                                                                                                                      ccTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGAATg              >  1:3251105/1‑150 (MQ=255)
                                                                                                                                      ccTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGAATg              >  1:4940687/1‑150 (MQ=255)
                                                                                                                                         gcgcATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGAATGcga           <  1:1666123/150‑1 (MQ=255)
                                                                                                                                             aTGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGAATGCgacgac       >  2:3150565/1‑150 (MQ=255)
                                                                                                                                              tGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCt                                   >  2:3880243/1‑121 (MQ=255)
                                                                                                                                              tGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCt                                   <  1:3880243/121‑1 (MQ=255)
                                                                                                                                                ccTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGAATGCGACGACCCg    <  1:2948027/150‑1 (MQ=255)
                                                                                                                                                 cTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGAATGCGACGACCCgg   <  1:2498929/150‑1 (MQ=255)
                                                                                                                                                  ttGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCtc                       >  2:4956117/1‑129 (MQ=255)
                                                                                                                                                  ttGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGAATGCGACGACCCGGt  >  2:362657/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCTGTTGATTGCCTCGATCGAATGCGACGACCCGGT  >  NC_002947/4995500‑4995795

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: