Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,851,477 T→C 8.7% G68G (GGT→GGC ilvD → dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,851,4770TC 8.7% 721.0 / 39.7 242G68G (GGT→GGCilvDdihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base T (96/125);  new base C (18/3);  total (114/128)
Fisher's exact test for biased strand distribution p-value = 1.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 251 total aligned reads displayed.

GTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGT  >  NC_002947/5851329‑5851623
                                                                                                                                                     |                                                                                                                                                  
gTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtg                                                                                                                                                   >  2:3939924/1‑150 (MQ=255)
gTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCAGCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtg                                                                                                                                                   <  1:4611712/150‑1 (MQ=255)
    ccTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtggcgt                                                                                                                                               <  2:479891/150‑1 (MQ=255)
    ccTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtggcgt                                                                                                                                               >  2:4522177/1‑150 (MQ=255)
     cTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGGGGCTCGCGAAATCGAACGCGCCGgtggcgtg                                                                                                                                              <  2:3647445/150‑1 (MQ=255)
         ggCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCCTGTGCACCTGAAGGACCTGGGCCCGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCa                                                                                                                                          >  2:3310718/1‑150 (MQ=255)
          gcgcgaCACCGGGATGAAGGACGAAGAATTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCaa                                                                                                                                         <  1:4003592/150‑1 (MQ=255)
             cgcCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGa                                                                                                                                      <  2:5052827/150‑1 (MQ=255)
              gcCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGaa                                                                                                                                     >  2:1859869/1‑150 (MQ=255)
               ccACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAAt                                                                                                                                    <  1:5277146/150‑1 (MQ=255)
                 aCCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTc                                                                                                                                  <  1:882314/150‑1 (MQ=255)
                    gggATGAAGGAAGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAac                                                                                                                               <  1:4383833/150‑1 (MQ=255)
                         gAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCAt                                                                                                                          >  2:3108021/1‑150 (MQ=255)
                           aGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATcg                                                                                                                        >  2:3807758/1‑150 (MQ=255)
                               cGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGaa                                                                                                                                     <  1:1418590/133‑1 (MQ=255)
                               cGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGaa                                                                                                                                     >  2:1418548/1‑133 (MQ=255)
                               cGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGt                                                                                                                    <  2:1975859/150‑1 (MQ=255)
                               cGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTCCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGt                                                                                                                    <  1:4830544/150‑1 (MQ=255)
                                gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                                                   >  1:4366195/1‑150 (MQ=255)
                                gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                                                   <  1:2563202/150‑1 (MQ=255)
                                gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                                                   >  2:501804/1‑150 (MQ=255)
                                      ttCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGAc                                                                                                             >  2:1269567/1‑150 (MQ=255)
                                         aaGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGc                                                                                                          >  2:431698/1‑150 (MQ=255)
                                            aaaCCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATc                                                                                                       >  2:297996/1‑150 (MQ=255)
                                            aaaCCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATc                                                                                                       <  1:4002167/150‑1 (MQ=255)
                                               ccGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGcc                                                                                                    >  2:2371807/1‑150 (MQ=255)
                                                   tcatcGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATgg                                                                                                <  1:895970/150‑1 (MQ=255)
                                                     atcgccatcgccAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGc                                                                                              <  2:4648330/150‑1 (MQ=255)
                                                     atcgccatcgccAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGc                                                                                              <  1:3076841/150‑1 (MQ=255)
                                                      tcgccatcgccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGcc                                                                                             >  1:1633461/1‑150 (MQ=255)
                                                      tcgccatcgccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGcc                                                                                             <  2:389708/150‑1 (MQ=255)
                                                       cgccatcaccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCa                                                                                            >  2:4156069/1‑150 (MQ=255)
                                                            tggtgaaCTCGTTCACCCAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCatgacg                                                                                       >  1:1166327/6‑150 (MQ=255)
                                                            tcgccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCacgacg                                                                                       <  2:3166548/150‑1 (MQ=255)
                                                                gaaCTCGTTCACCCAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCAt                                                                                   >  2:1588682/2‑150 (MQ=255)
                                                                   cTCGTTCACCCAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCt                                                                                >  1:5001597/1‑150 (MQ=255)
                                                                   cTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCt                                                                                <  1:2846652/150‑1 (MQ=255)
                                                                   cTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCt                                                                                <  2:2909851/150‑1 (MQ=255)
                                                                      gTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGTGTGGCCAAGGAATTc                                                                                                                                  >  2:2659105/1‑97 (MQ=255)
                                                                            ccAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCt                                                                       >  2:3599147/1‑150 (MQ=255)
                                                                            ccAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCt                                                                       >  1:3466051/1‑150 (MQ=255)
                                                                            ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                                                       <  1:5263871/150‑1 (MQ=255)
                                                                            ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                                                       >  1:3911234/1‑150 (MQ=255)
                                                                              aGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGc                                                                     >  2:645670/1‑150 (MQ=255)
                                                                                ttCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCa                                                                   >  1:2820247/1‑150 (MQ=255)
                                                                                    tACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGcc                                                               <  1:456485/150‑1 (MQ=255)
                                                                                       cGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCcgcg                                                            >  1:4166005/1‑150 (MQ=255)
                                                                                         ggCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGaa                                                          >  2:2108601/1‑150 (MQ=255)
                                                                                          gCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGaaa                                                         <  2:4999441/150‑1 (MQ=255)
                                                                                           ccATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAAt                                                        >  2:2832672/1‑150 (MQ=255)
                                                                                              tgtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAt                                                     >  2:2166269/1‑150 (MQ=255)
                                                                                              tgtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAt                                                     >  1:2714424/1‑150 (MQ=255)
                                                                                               gtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAtt                                                    <  1:582697/150‑1 (MQ=255)
                                                                                               gtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAtt                                                    >  1:458948/1‑150 (MQ=255)
                                                                                                tgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTg                                                   <  2:3505098/150‑1 (MQ=255)
                                                                                                  cACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGcc                                                 >  2:2822861/1‑150 (MQ=255)
                                                                                                   aCCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCg                                                >  1:2464588/1‑150 (MQ=255)
                                                                                                     cTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCCTGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGAc                                              >  2:1423691/1‑150 (MQ=255)
                                                                                                        aaGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGcc                                           >  2:1547191/1‑150 (MQ=255)
                                                                                                          ggACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGt                                         >  1:3251488/1‑150 (MQ=255)
                                                                                                           gACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTg                                        <  1:4092056/150‑1 (MQ=255)
                                                                                                           gACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGActcggtt                                        >  1:1921664/1‑144 (MQ=255)
                                                                                                              cTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGaa                                     >  1:1961143/1‑150 (MQ=255)
                                                                                                               tGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAAt                                    >  1:3540324/1‑150 (MQ=255)
                                                                                                               tGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAAt                                    >  1:5214760/1‑150 (MQ=255)
                                                                                                               tGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAAt                                    >  2:755326/1‑150 (MQ=255)
                                                                                                                gggCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATa                                   >  1:4451196/1‑150 (MQ=255)
                                                                                                                gggCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATa                                   >  1:3684996/1‑150 (MQ=255)
                                                                                                                cggCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAATa                                   >  1:70179/2‑150 (MQ=255)
                                                                                                                 ggCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGTCGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATAc                                  <  1:3358019/150‑1 (MQ=255)
                                                                                                                 ggCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATAc                                  <  1:4509911/150‑1 (MQ=255)
                                                                                                                  gCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa                                 <  1:3492423/150‑1 (MQ=255)
                                                                                                                  gCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa                                 >  1:4928783/1‑150 (MQ=255)
                                                                                                                  gCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa                                 <  1:1005161/150‑1 (MQ=255)
                                                                                                                  gCCAGCGGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa                                 >  2:2290805/1‑150 (MQ=255)
                                                                                                                  gCCAGCGGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa                                 <  2:4597029/150‑1 (MQ=255)
                                                                                                                   ccAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACAt                                <  2:1851697/150‑1 (MQ=255)
                                                                                                                   ccAGCCCGTCGCGCGCGAAATTGACGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAatatat                                >  2:3755425/1‑150 (MQ=255)
                                                                                                                          gtggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCCACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAAt                         >  1:1818449/1‑150 (MQ=255)
                                                                                                                          gtggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAAt                         >  2:3488891/1‑150 (MQ=255)
                                                                                                                          ggggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAAt                         <  2:4396888/148‑1 (MQ=255)
                                                                                                                            ggCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATgc                       <  1:3310718/150‑1 (MQ=255)
                                                                                                                             gCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATgcg                      <  2:1374163/150‑1 (MQ=255)
                                                                                                                              cTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATgcgc                     <  2:4928783/150‑1 (MQ=255)
                                                                                                                               tCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCacga                                                                                         <  2:2494954/81‑1 (MQ=255)
                                                                                                                               tCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCacga                                                                                         >  1:2494954/1‑81 (MQ=255)
                                                                                                                               tCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGt                                         >  1:1302451/1‑129 (MQ=255)
                                                                                                                               tCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGt                                         <  2:1302451/129‑1 (MQ=255)
                                                                                                                               tCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa                                 >  1:3842304/1‑137 (MQ=255)
                                                                                                                                   gAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGaa                                                          <  2:4142079/108‑1 (MQ=255)
                                                                                                                                   gAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTgc                <  1:3825233/150‑1 (MQ=255)
                                                                                                                                   gAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTgc                <  1:5068917/150‑1 (MQ=255)
                                                                                                                                       tCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGACTACATGGTCAATGCGCACTGCGCCg            >  2:4727274/1‑150 (MQ=255)
                                                                                                                                         gAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGAc          <  2:2712821/150‑1 (MQ=255)
                                                                                                                                             gcgcCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCa                                                                                    >  1:1360478/1‑72 (MQ=255)
                                                                                                                                             gcgcCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCa      <  2:4487851/150‑1 (MQ=255)
                                                                                                                                                gcCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAtt                                                    >  1:1856093/1‑102 (MQ=255)
                                                                                                                                                 ccGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCg   >  2:3428248/1‑150 (MQ=255)
                                                                                                                                                  cGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGt  <  2:5414235/150‑1 (MQ=255)
                                                                                                                                                  cGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGAAGCCATCGt  <  1:2290805/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGT  >  NC_002947/5851329‑5851623

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: