Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,995,613 C→T 6.9% G109G (GGC→GGT 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,6130CT 6.9% 479.4 / 17.8 158G109G (GGC→GGTbkdABbranched‑chain alpha‑keto acid dehydrogenase complex subunit beta
Reads supporting (aligned to +/- strand):  ref base C (70/77);  new base T (5/6);  total (75/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 168 total aligned reads displayed.

ATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCT  >  NC_002947/4995464‑4995762
                                                                                                                                                     |                                                                                                                                                     
aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                                                       <  1:2759034/150‑1 (MQ=255)
aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                                                       <  1:3242831/150‑1 (MQ=255)
aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                                                       >  1:1742710/1‑150 (MQ=255)
aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                                                       >  1:973114/1‑150 (MQ=255)
aTCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTcggccggc                                                                                                                                                       <  1:200868/150‑1 (MQ=255)
   tcCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAg                                                                                                                                                    >  2:458384/1‑150 (MQ=255)
   tcCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAg                                                                                                                                                    >  2:2191798/1‑150 (MQ=255)
   tcCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAg                                                                                                                                                    >  1:2040095/1‑150 (MQ=255)
     cGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGtt                                                                                                                                                  <  1:1409203/150‑1 (MQ=255)
         agCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCAtt                                                                                                                                              <  1:2126525/150‑1 (MQ=255)
           cGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGc                                                                                                                                            <  1:4948026/150‑1 (MQ=255)
           cGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGc                                                                                                                                            <  1:966108/150‑1 (MQ=255)
            ggCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGcc                                                                                                                                           <  2:392689/150‑1 (MQ=255)
             gCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGccc                                                                                                                                          <  1:1706394/150‑1 (MQ=255)
              cATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGcccc                                                                                                                                         >  2:1237090/1‑150 (MQ=255)
               aTCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCg                                                                                                                                        >  1:2958879/1‑150 (MQ=255)
                  gtcgGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTg                                                                                                                                     >  2:2498929/1‑150 (MQ=255)
                   tcgGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGa                                                                                                                                    <  2:1135061/150‑1 (MQ=255)
                        aCCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTg                                                                                                                               >  2:3054941/1‑150 (MQ=255)
                         ccgccgTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTgc                                                                                                                              <  2:3110205/150‑1 (MQ=255)
                         ccgccgTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTgc                                                                                                                              <  1:4970579/150‑1 (MQ=255)
                                    aTGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTgc                                                                                                                   >  2:661287/1‑150 (MQ=255)
                                    aTGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTgc                                                                                                                   <  2:2443795/150‑1 (MQ=255)
                                            cTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCAGACTACTTCTACCCGGCATCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGcggcggcgg                                                                                                           >  1:3795084/1‑150 (MQ=255)
                                             tATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCggcggcggc                                                                                                          >  2:1939577/1‑150 (MQ=255)
                                             tATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCggcggcggc                                                                                                          >  2:1939523/1‑150 (MQ=255)
                                                 gCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCt                                                                                                      >  1:253211/1‑150 (MQ=255)
                                                  ccTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTa                                                                                                     <  1:2641026/150‑1 (MQ=255)
                                                   cTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTAt                                                                                                    >  1:4548092/1‑150 (MQ=255)
                                                   cTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTAt                                                                                                    <  2:2017130/150‑1 (MQ=255)
                                                   cTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTAt                                                                                                    <  1:925484/150‑1 (MQ=255)
                                                    tGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATg                                                                                                   >  1:4664032/1‑150 (MQ=255)
                                                       gcCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcg                                                                                                >  2:2092308/1‑150 (MQ=255)
                                                       gcCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcg                                                                                                <  2:4106834/150‑1 (MQ=255)
                                                       gcCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATGGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcg                                                                                                <  2:4106657/150‑1 (MQ=255)
                                                        cccGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcgg                                                                                               >  1:3134409/1‑150 (MQ=255)
                                                        cccGGTGGTGGAGATCCAGTTCGCCGACTAATTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcgg                                                                                               <  1:4764443/150‑1 (MQ=255)
                                                         ccGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcggc                                                                                              >  1:3752820/1‑150 (MQ=255)
                                                         ccGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATggcggc                                                                                              >  1:3008840/1‑150 (MQ=255)
                                                           ggtggtGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCa                                                                                            <  2:738320/150‑1 (MQ=255)
                                                              ggtggAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGAc                                                                                         <  2:2784765/150‑1 (MQ=255)
                                                                     gtCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGc                                                                                  <  1:843708/149‑1 (MQ=255)
                                                                      tgcAGTTTGCGGATTATGTCTATCCGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGcc                                                                                 >  1:3957254/3‑150 (MQ=255)
                                                                      tgcAGTTTGCGGATTATGTCTATCCGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGcc                                                                                 >  2:1495558/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:3550703/1‑150 (MQ=255)
                                                                          gTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCagag                                                                             >  1:2462877/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)
                                                                                  actactTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAg                                                                     <  1:723961/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)
                                                                                     actTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGa                                                                  <  1:4629170/150‑1 (MQ=255)
                                                                                     actTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGa                                                                  <  1:1766764/150‑1 (MQ=255)
                                                                                      ctTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGAt                                                                 <  2:3248886/150‑1 (MQ=255)
                                                                                      ctTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGAt                                                                 <  2:3248705/150‑1 (MQ=255)
                                                                                      ctTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGAt                                                                 >  1:2045979/1‑150 (MQ=255)
                                                                                       ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATg                                                                >  2:1715854/1‑150 (MQ=255)
                                                                                       ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATg                                                                <  1:4262922/150‑1 (MQ=255)
                                                                                       ttCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATg                                                                >  2:3422007/1‑150 (MQ=255)
                                                                                        tCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGt                                                               <  2:1353789/150‑1 (MQ=255)
                                                                                        tCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGt                                                               <  2:253211/150‑1 (MQ=255)
                                                                                         cTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGtt                                                              <  2:1648699/150‑1 (MQ=255)
                                                                                         cTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGtt                                                              <  2:2464567/150‑1 (MQ=255)
                                                                                           atccGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCa                                                            <  1:1016386/148‑1 (MQ=255)
                                                                                             ccGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCAcc                                                          >  2:720908/1‑150 (MQ=255)
                                                                                              cGGCCTATGACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCAccc                                                         >  1:2017743/1‑150 (MQ=255)
                                                                                                gcctatgACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCACCCat                                                       <  2:2700475/144‑2 (MQ=255)
                                                                                                gcctatgACCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCACCCat                                                       <  2:3590014/144‑2 (MQ=255)
                                                                                                gCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAg                                                       <  1:4293393/150‑1 (MQ=255)
                                                                                                  cTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGgt                                                     <  2:3567172/150‑1 (MQ=255)
                                                                                                   tCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGgtg                                                    >  2:4911614/1‑150 (MQ=255)
                                                                                                   tCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGAGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGgtg                                                    >  2:3072271/1‑150 (MQ=255)
                                                                                                     cGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACACAGGtgtg                                                  >  2:1352952/1‑150 (MQ=255)
                                                                                                       aCCAGATCGTCTCGGAAGCGGCGCGCTTGCGCTATCGTTCGGCGGGTGAATTTACATGCCCTCTGGTGATCCGTATGCCGTCAGGCGGCGGCATCTATGGCGGGCAGACACACAGCCAGAGCCCGGAGGCGCTTTTCACCCATgtttcgg                                                <  2:1172034/150‑6 (MQ=255)
                                                                                                          aGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGcc                                             <  1:2377351/150‑1 (MQ=255)
                                                                                                           gATCGTCTCCGAGCTGGCCCGCCTGCGTGACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCt                                            <  2:3434783/150‑1 (MQ=255)
                                                                                                                  tcCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCAcc                                     >  1:2811567/1‑150 (MQ=255)
                                                                                                                  tcCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCAcc                                     >  2:1091820/1‑150 (MQ=255)
                                                                                                                  tcCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCAcc                                     >  1:4434922/1‑150 (MQ=255)
                                                                                                                    cGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGt                                   >  2:4867128/1‑150 (MQ=255)
                                                                                                                          ggCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGcc                             <  2:2040264/150‑1 (MQ=255)
                                                                                                                             ccGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTc                          >  2:1373562/1‑150 (MQ=255)
                                                                                                                                  tGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAAcc                     <  2:4249204/150‑1 (MQ=255)
                                                                                                                                     gTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCtt                  >  1:728447/1‑150 (MQ=255)
                                                                                                                                     gTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCtt                  <  1:3123017/150‑1 (MQ=255)
                                                                                                                                      ttACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCTGTCCAACCCTTa                 >  2:3830602/1‑150 (MQ=255)
                                                                                                                                      ttACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTa                 >  1:1946922/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)
                                                                                                                                          cGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGAc             >  2:2527665/1‑150 (MQ=255)
                                                                                                                                                gccggcGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCTTCTATGGCTGCCAGACTCACAGCCAGAGCCCGGAAGCGTTGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCaaa       <  2:3795084/150‑1 (MQ=255)
                                                                                                                                                   ggcGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGc    >  1:2694364/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCTCCGAGAGCGGCATCGTCGGTACCGCCGTGGGCATGGGTGCCTATGGCCTGCGCCCGGTGGTGGAGATCCAGTTCGCCGACTACTTCTACCCGGCCTCCGACCAGATCGTCTCCGAGCTGGCCCGCCTGCGTTACCGTTCGGCCGGCGAGTTCATTGCCCCGCTGACCCTGCGCATGCCTTGCGGCGGCGGCATCTATGGCGGCCAGACTCACAGCCAGAGCCCGGAAGCGATGTTCACCCAGGTGTGCGGCCTGCGCACCGTGATGCCGTCCAACCCTTATGACGCCAAAGGCCT  >  NC_002947/4995464‑4995762

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: