Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,851,471 C→A 8.3% R66R (CGC→CGA) ‡ ilvD → dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,851,4710CA 8.3% 789.2 / 43.2 242R66R (CGC→CGA) ‡ilvDdihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base C (95/127);  new base A (17/3);  total (112/130)
Fisher's exact test for biased strand distribution p-value = 2.91e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 246 total aligned reads displayed.

CCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCA  >  NC_002947/5851326‑5851619
                                                                                                                                                 |                                                                                                                                                     
 ccGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCgg                                                                                                                                                 <  1:2449412/150‑1 (MQ=255)
   gTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtg                                                                                                                                               >  2:3939924/1‑150 (MQ=255)
    tGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtgg                                                                                                                                              <  1:141070/150‑1 (MQ=255)
       ccTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtggcgt                                                                                                                                           >  2:2124941/1‑150 (MQ=255)
       ccTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtggcgt                                                                                                                                           >  2:4522177/1‑150 (MQ=255)
       ccTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtggcgt                                                                                                                                           <  2:479891/150‑1 (MQ=255)
        cTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtggcgtg                                                                                                                                          <  2:282693/150‑1 (MQ=255)
        cTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGGGGCTCGCGAAATCGAACGCGCCGgtggcgtg                                                                                                                                          <  2:3647445/150‑1 (MQ=255)
             gcgcgaCACCGGGATGAAGGACGAAGAATTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCaa                                                                                                                                     <  1:4003592/150‑1 (MQ=255)
                cgcCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGa                                                                                                                                  <  2:5052827/150‑1 (MQ=255)
                      cGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCaa                                                                                                                            <  2:3018301/150‑1 (MQ=255)
                       gggATGAAGGAAGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAac                                                                                                                           <  1:4383833/150‑1 (MQ=255)
                              aGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATcg                                                                                                                    >  2:3807758/1‑150 (MQ=255)
                                  cGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGaa                                                                                                                                 <  1:1418590/133‑1 (MQ=255)
                                  cGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGaa                                                                                                                                 >  2:1418548/1‑133 (MQ=255)
                                   gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                                               >  2:1431728/1‑150 (MQ=255)
                                   gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                                               >  2:501734/1‑150 (MQ=255)
                                   gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                                               >  1:5171544/1‑150 (MQ=255)
                                   gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                                               >  1:4549498/1‑150 (MQ=255)
                                   gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                                               <  1:4061872/150‑1 (MQ=255)
                                        cTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGa                                                                                                          <  2:1274754/150‑1 (MQ=255)
                                        cTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGa                                                                                                          <  2:4590629/150‑1 (MQ=255)
                                         ttCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGAc                                                                                                         >  2:2501513/1‑150 (MQ=255)
                                           cAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACgg                                                                                                       <  2:3375758/150‑1 (MQ=255)
                                            aaGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGc                                                                                                      >  2:431698/1‑150 (MQ=255)
                                               aaaCCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATc                                                                                                   <  1:1936842/150‑1 (MQ=255)
                                                    gATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCAt                                                                                              >  1:4775321/1‑150 (MQ=255)
                                                     atcatcGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATg                                                                                             <  1:3954847/150‑1 (MQ=255)
                                                      tcatcGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATgg                                                                                            <  1:617716/150‑1 (MQ=255)
                                                        atcgccatcgccAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGc                                                                                          <  1:3076841/150‑1 (MQ=255)
                                                        atcgccatcgccAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGc                                                                                          <  2:4648330/150‑1 (MQ=255)
                                                         tcgccatcgccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGcc                                                                                         <  2:389708/150‑1 (MQ=255)
                                                          cgccatcaccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCa                                                                                        >  2:4156069/1‑150 (MQ=255)
                                                               tcgccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCacgacg                                                                                   >  2:2250593/1‑150 (MQ=255)
                                                                cgccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACgg                                                                                  <  2:2708055/150‑1 (MQ=255)
                                                                    aaCTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATg                                                                              <  1:5024956/150‑1 (MQ=255)
                                                                    aaCTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGGGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATg                                                                              <  2:4925519/150‑1 (MQ=255)
                                                                      cTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCt                                                                            <  1:2846652/150‑1 (MQ=255)
                                                                         gTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTa                                                                         >  1:247858/1‑150 (MQ=255)
                                                                            cACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTc                                                                      <  2:2339511/150‑1 (MQ=255)
                                                                              cccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGGGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGc                                                                    <  1:5261630/150‑1 (MQ=255)
                                                                               ccAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCt                                                                   >  2:1380313/1‑150 (MQ=255)
                                                                               ccAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCt                                                                   >  1:3466051/1‑150 (MQ=255)
                                                                               ccAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCt                                                                   >  2:3599147/1‑150 (MQ=255)
                                                                               ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                                                   >  1:1865496/1‑150 (MQ=255)
                                                                               ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                                                   <  2:1734000/150‑1 (MQ=255)
                                                                               ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                                                   <  1:5263872/150‑1 (MQ=255)
                                                                               ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                                                   >  1:3911234/1‑150 (MQ=255)
                                                                               ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                                                   <  1:3835354/150‑1 (MQ=255)
                                                                               caaGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGGGGCCACGGACTTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                                                   >  1:2627968/3‑150 (MQ=255)
                                                                                cAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTg                                                                  >  1:4927375/1‑150 (MQ=255)
                                                                                   ttCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCa                                                               >  1:2820247/1‑150 (MQ=255)
                                                                                       tACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGcc                                                           <  1:2271229/150‑1 (MQ=255)
                                                                                       tACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGcc                                                           >  2:2481413/1‑150 (MQ=255)
                                                                                          cGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCcgcg                                                        >  1:4166005/1‑150 (MQ=255)
                                                                                             gCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGaaa                                                     <  2:4999441/150‑1 (MQ=255)
                                                                                              ccATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAAt                                                    <  2:1904078/150‑1 (MQ=255)
                                                                                               cATGTGCACCTGAAGGACCTGGGCCAGCTGGGGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATc                                                   <  2:2853470/150‑1 (MQ=255)
                                                                                                 tgtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAt                                                 >  1:2714424/1‑150 (MQ=255)
                                                                                                 tgtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAt                                                 <  1:4448345/150‑1 (MQ=255)
                                                                                                  gtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAtt                                                >  1:458948/1‑150 (MQ=255)
                                                                                                  gtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAtt                                                <  2:4699411/150‑1 (MQ=255)
                                                                                                   tgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTg                                               <  1:1721428/150‑1 (MQ=255)
                                                                                                   tgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTg                                               <  2:3505098/150‑1 (MQ=255)
                                                                                                     cACCTGAAGGACCAGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGcc                                             <  2:2754602/150‑1 (MQ=255)
                                                                                                       ccTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGa                                           <  2:3306019/150‑1 (MQ=255)
                                                                                                           aaGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGcc                                       >  2:1547191/1‑150 (MQ=255)
                                                                                                             ggACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAATGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGt                                     <  2:1111374/150‑1 (MQ=255)
                                                                                                             ggACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGt                                     >  1:3251488/1‑150 (MQ=255)
                                                                                                             ggACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGt                                     <  2:646879/150‑1 (MQ=255)
                                                                                                              gACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTg                                    <  1:496351/150‑1 (MQ=255)
                                                                                                               aCCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTgg                                   >  1:839135/1‑150 (MQ=255)
                                                                                                                  tGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAAt                                >  2:755326/1‑150 (MQ=255)
                                                                                                                  tGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAAt                                >  1:3540324/1‑150 (MQ=255)
                                                                                                                   gggCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATa                               >  1:3684996/1‑150 (MQ=255)
                                                                                                                    ggCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGTCGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATAc                              <  1:3358019/150‑1 (MQ=255)
                                                                                                                    ggCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATAc                              <  1:4509911/150‑1 (MQ=255)
                                                                                                                     gCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa                             <  1:1005161/150‑1 (MQ=255)
                                                                                                                     gCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa                             <  1:3492423/150‑1 (MQ=255)
                                                                                                                     gCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa                             <  2:5214760/150‑1 (MQ=255)
                                                                                                                      ccAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACAt                            <  1:4598443/150‑1 (MQ=255)
                                                                                                                      ccAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACAt                            <  2:1851697/150‑1 (MQ=255)
                                                                                                                      ccAGCTGGGGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACAt                            <  2:2087264/150‑1 (MQ=255)
                                                                                                                      ccAGCCCGTCGCGCGCGAAATTGACGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAatatat                            >  2:3755425/1‑150 (MQ=255)
                                                                                                                           tggtggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCa                       <  1:1233583/150‑1 (MQ=255)
                                                                                                                           tggtggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCa                       >  1:2495738/1‑150 (MQ=255)
                                                                                                                             gtggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCCACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAAt                     >  1:1818449/1‑150 (MQ=255)
                                                                                                                             gtggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAAt                     <  1:5363679/150‑1 (MQ=255)
                                                                                                                             gtggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAAt                     >  2:3488891/1‑150 (MQ=255)
                                                                                                                             ggggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAAt                     <  2:4396888/148‑1 (MQ=255)
                                                                                                                               ggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATgc                   <  1:3310718/150‑1 (MQ=255)
                                                                                                                                  tCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGt                                     >  1:1302451/1‑129 (MQ=255)
                                                                                                                                  tCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGt                                     <  2:1302451/129‑1 (MQ=255)
                                                                                                                                  tCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCa                >  1:89402/1‑150 (MQ=255)
                                                                                                                                      gAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCTGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTgc            >  2:2917708/1‑150 (MQ=255)
                                                                                                                                      gAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGaa                                                      <  2:4142079/108‑1 (MQ=255)
                                                                                                                                      gAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTgc            <  1:3825233/150‑1 (MQ=255)
                                                                                                                                             aaCGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACg     >  1:4291689/1‑150 (MQ=255)
                                                                                                                                              aCGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGc    <  2:3117623/150‑1 (MQ=255)
                                                                                                                                                gcgcCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCa  <  2:4487851/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCA  >  NC_002947/5851326‑5851619

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: