Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,851,672 C→A 8.3% A133A (GCC→GCA ilvD → dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,851,6720CA 8.3% 748.6 / 40.1 229A133A (GCC→GCAilvDdihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base C (105/105);  new base A (10/9);  total (115/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

GCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGA  >  NC_002947/5851523‑5851821
                                                                                                                                                     |                                                                                                                                                     
      ggCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCATCTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTgcgc                                                                                                                                                 >  2:457998/1‑150 (MQ=255)
      ggCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTgcgc                                                                                                                                                 <  2:5255872/150‑1 (MQ=255)
         cACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTg                                                                                                                                              >  1:1145783/1‑150 (MQ=255)
            gacgGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACAGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAAc                                                                                                                                           >  1:1502290/1‑150 (MQ=255)
            gacgGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCaat                                                                                                                                           >  1:2265588/1‑149 (MQ=255)
                  aTGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCg                                                                                                                                     <  1:2108601/150‑1 (MQ=255)
                    gCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGt                                                                                                                                   >  1:3595100/1‑150 (MQ=255)
                      tGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGa                                                                                                                                 >  1:462901/1‑150 (MQ=255)
                      tGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGa                                                                                                                                 <  2:4366195/150‑1 (MQ=255)
                        tACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATc                                                                                                                               <  1:1431728/150‑1 (MQ=255)
                          cTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCtt                                                                                                                             >  2:3946881/1‑150 (MQ=255)
                           tCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTc                                                                                                                            >  2:187950/1‑150 (MQ=255)
                           tCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTc                                                                                                                            <  2:247858/150‑1 (MQ=255)
                           tCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTc                                                                                                                            <  2:3911272/150‑1 (MQ=255)
                            cGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCg                                                                                                                           <  1:645670/150‑1 (MQ=255)
                            cGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCg                                                                                                                           <  1:2481413/150‑1 (MQ=255)
                             gCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGt                                                                                                                          >  2:4393671/1‑150 (MQ=255)
                                  cAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCcg                                                                                                                     >  1:718571/1‑150 (MQ=255)
                                   aaGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCcgg                                                                                                                    >  1:3250791/1‑150 (MQ=255)
                                         cgagctgaTTTCCGACTCGGTTGAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGTTTCGGGCGGgcc                                                                                                              <  2:5001597/143‑1 (MQ=255)
                                           aaaTCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGa                                                                                                            <  1:4156069/150‑1 (MQ=255)
                                           aaaTCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGa                                                                                                            <  1:3428647/150‑1 (MQ=255)
                                            aaTCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGAt                                                                                                           <  1:717293/150‑1 (MQ=255)
                                            aaTCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGAt                                                                                                           <  1:494306/150‑1 (MQ=255)
                                              tCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATgg                                                                                                         <  1:3488891/150‑1 (MQ=255)
                                                aTTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGaa                                                                                                       <  2:89402/150‑1 (MQ=255)
                                                aTTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGaa                                                                                                       <  2:3216347/150‑1 (MQ=255)
                                                 ttGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAATATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAg                                                                                                      >  2:3283598/1‑150 (MQ=255)
                                                  tGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGc                                                                                                     <  1:2917708/150‑1 (MQ=255)
                                                    ccGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCg                                                                                                   >  1:3186152/1‑150 (MQ=255)
                                                    ccGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCg                                                                                                   <  2:2371923/150‑1 (MQ=255)
                                                    ccGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCg                                                                                                   <  1:2832672/150‑1 (MQ=255)
                                                      gACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGc                                                                                                 >  1:3685158/1‑150 (MQ=255)
                                                         tcggttgAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGTTTCGGGCGGGCCGATGGAAGCAGGCAAg                                                                                              <  2:4895703/144‑1 (MQ=255)
                                                         tcggttgAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGTTTCGGGCGGGCCGATGGAAGCAGGCAAg                                                                                              >  2:4151113/7‑150 (MQ=255)
                                                         tcggttgAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGTTTCGGGCGGGCCGATGGAAGCAGGCAAg                                                                                              <  2:3748003/144‑1 (MQ=255)
                                                          ccGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGa                                                                                             <  1:1547191/150‑1 (MQ=255)
                                                          ccGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGa                                                                                             <  1:5280636/150‑1 (MQ=255)
                                                           cGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGAc                                                                                            >  1:621157/1‑150 (MQ=255)
                                                            gTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGAcc                                                                                           >  1:1640253/1‑150 (MQ=255)
                                                              ggAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCaa                                                                                         >  1:1984610/1‑150 (MQ=255)
                                                               gAATATATGGTCAATGCCCATTGCGCCGACGCGATGGTCTGCATTTCCAACTGCGACAAGATCACTCCCGGCATGCTGATGGCAGCACTCCGCCTCAATATCCCGGTCGTGTTCGTTTCgggcgg                                                                                                                 >  1:1279463/1‑125 (MQ=255)
                                                                 aTACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCt                                                                                      >  2:2170572/1‑150 (MQ=255)
                                                                 aTACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCt                                                                                      <  2:289741/150‑1 (MQ=255)
                                                                  tACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTg                                                                                     <  2:3684996/150‑1 (MQ=255)
                                                                  tACATGGTCAAGGCGCACGGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTg                                                                                     <  2:3842920/150‑1 (MQ=255)
                                                                   aCATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTgg                                                                                    <  1:2170572/150‑1 (MQ=255)
                                                                     aTGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgcc                                                                                  <  1:1715489/150‑1 (MQ=255)
                                                                     aTGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgcc                                                                                  >  2:5152968/1‑150 (MQ=255)
                                                                     aTGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgcc                                                                                  >  1:2358662/1‑150 (MQ=255)
                                                                      tGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgcca                                                                                 >  2:3618142/1‑150 (MQ=255)
                                                                       ggTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgccag                                                                                >  2:2052697/1‑150 (MQ=255)
                                                                        gTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgccagc                                                                               <  2:2714424/150‑1 (MQ=255)
                                                                         tCAATGCGCACTGCGCCGACGCCATCGTCTATATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGgccagcc                                                                              >  2:4222880/1‑150 (MQ=255)
                                                                            aTGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACg                                                                           >  2:4496075/1‑150 (MQ=255)
                                                                             tGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACgg                                                                          >  1:3763169/1‑150 (MQ=255)
                                                                              gcgcACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGc                                                                         >  1:4562676/1‑150 (MQ=255)
                                                                              gcgcACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGc                                                                         <  2:4166005/150‑1 (MQ=255)
                                                                                gcACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCt                                                                       >  1:4875378/1‑150 (MQ=255)
                                                                                gcACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCt                                                                       >  1:4697685/1‑150 (MQ=255)
                                                                                    tGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGAc                                                                   <  2:3540324/150‑1 (MQ=255)
                                                                                       gcCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTg                                                                <  2:2495738/150‑1 (MQ=255)
                                                                                       gcCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTg                                                                >  1:1727893/1‑150 (MQ=255)
                                                                                        ccGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATcc                                                                                                                                       <  1:823434/78‑1 (MQ=255)
                                                                                        ccGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATcc                                                                                                                                       >  2:823434/1‑78 (MQ=255)
                                                                                           aCGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCg                                                            >  2:4600645/1‑150 (MQ=255)
                                                                                             gCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGAt                                                          <  2:1860766/150‑1 (MQ=255)
                                                                                               cATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGc                                                        >  2:4961022/1‑150 (MQ=255)
                                                                                               cATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGc                                                        >  2:1049081/1‑150 (MQ=255)
                                                                                               cATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCGGGTCGATGc                                                        >  2:5391603/1‑150 (MQ=255)
                                                                                                aTCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGcc                                                       >  1:1425199/1‑150 (MQ=255)
                                                                                                 tCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCa                                                      <  2:57377/150‑1 (MQ=255)
                                                                                                   gTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATg                                                    <  1:3920295/150‑1 (MQ=255)
                                                                                                          tctcCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATcg                                             >  1:2009509/1‑150 (MQ=255)
                                                                                                              cAACGGCGACAAGATCACCCCCGGCATGCTGCTGGCCGCCCTGCGCCTGAACCGCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGc                                         <  2:4291689/150‑1 (MQ=255)
                                                                                                               aaCTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGcc                                        >  1:1042236/1‑150 (MQ=255)
                                                                                                               aaCTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGcc                                        >  1:4104817/1‑150 (MQ=255)
                                                                                                               aaCTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGcc                                        <  1:4131661/150‑1 (MQ=255)
                                                                                                                  tGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGAc                                     <  2:1818449/150‑1 (MQ=255)
                                                                                                                     gACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTcc                                  >  1:5129495/1‑150 (MQ=255)
                                                                                                                     gACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTcc                                  >  1:1207623/1‑150 (MQ=255)
                                                                                                                     gACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTcc                                  >  1:1399791/1‑150 (MQ=255)
                                                                                                                      aCAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGaa                                                                                                       >  1:436754/1‑80 (MQ=255)
                                                                                                                         aGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACgg                              >  2:5198243/1‑150 (MQ=255)
                                                                                                                                ccccGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCcga                       >  1:4775850/1‑150 (MQ=255)
                                                                                                                                ccccGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCcga                       >  2:3767242/1‑150 (MQ=255)
                                                                                                                                  ccGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCcgacg                     >  2:4593288/1‑150 (MQ=255)
                                                                                                                                  ccGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCcgacg                     >  2:2358025/1‑150 (MQ=255)
                                                                                                                                   cGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCcgacga                    <  1:3115610/150‑1 (MQ=255)
                                                                                                                                       aTGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGaaa                <  1:4536933/150‑1 (MQ=255)
                                                                                                                                       aTGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGaaa                <  2:1555901/150‑1 (MQ=255)
                                                                                                                                         gCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCcggcgg                                                                                                                 <  2:2440140/51‑1 (MQ=255)
                                                                                                                                          cTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTc             <  2:1502290/150‑1 (MQ=255)
                                                                                                                                          cTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTc             <  2:64272/150‑1 (MQ=255)
                                                                                                                                          cTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTc             <  2:3989257/150‑1 (MQ=255)
                                                                                                                                            gATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGc           >  1:2472818/1‑150 (MQ=255)
                                                                                                                                             aTGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGcc          <  2:4537044/150‑1 (MQ=255)
                                                                                                                                             aTGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGcc          <  1:3416449/150‑1 (MQ=255)
                                                                                                                                             aTGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGcc          <  2:1145783/150‑1 (MQ=255)
                                                                                                                                                 ccgccCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGaaa                <  2:1626387/140‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCAATGCGCACTGCGCCGACGCCATCGTCTGTATCTCCAACTGCGACAAGATCACCCCCGGCATGCTGATGGCCGCCCTGCGCCTGAACATCCCGGTGATCTTCGTTTCCGGCGGCCCGATGGAAGCCGGCAAGACCAAGCTGGCCAGCCACGGCCTGGACCTGGTCGATGCCATGGTCATCGCGGCCGACTCCACGGCCTCCGACGAGAAAGTCGCCGAATACGA  >  NC_002947/5851523‑5851821

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: