Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,851,450 G→C 8.4% L59L (CTG→CTC ilvD → dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,851,4500GC 8.4% 775.1 / 44.7 238L59L (CTG→CTCilvDdihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base G (95/123);  new base C (18/2);  total (113/125)
Fisher's exact test for biased strand distribution p-value = 7.20e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 253 total aligned reads displayed.

ACCCAAGGCCGCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCA  >  NC_002947/5851301‑5851598
                                                                                                                                                     |                                                                                                                                                    
aCCCAAGGCCGCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCCACTCGTTCACCCAGTTCGTACCGGGCCCTGTGCACCTGAAGGACCTGGGCCAGCtg                                                                                                                                                      >  1:2087264/1‑150 (MQ=255)
  ccACGGCCGCAATATGGCTGGCGCGCGCGGCCTATGGCGCGCCACCGGCATGAAGGATGAGGATTTTGGCAAGCCAATTATTGCTGTGGTGAACTCGTTCACCCAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCtcgt                                                                                                                                                    >  2:884198/1‑150 (MQ=255)
  ccACGGCCGCAATATGGCTGGCGCGCGCGGCCTATGGCGCGCCACCGGCATGAAGGATGAGGATTTTGGCAAGCCAATTATTGCTGTGGTGAACTCGTTCACCCAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCtcgt                                                                                                                                                    >  2:775419/1‑150 (MQ=255)
  ccACGGCCGCAATATGGCTGGCGCGCGCGGCCTATGGCGCGCCACCGGCATGAAGGATGAGGATTTTGGCAAGCCAATTATTGCTGTGGTGAACTCGTTCACCCAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCtcgt                                                                                                                                                    >  1:2259300/1‑150 (MQ=255)
  ccAAGGCCGCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCtggt                                                                                                                                                    >  2:1021278/1‑150 (MQ=255)
     aGGCCGCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGc                                                                                                                                                 >  1:1374163/1‑150 (MQ=255)
     aGGCCGCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGc                                                                                                                                                 >  1:4511188/1‑150 (MQ=255)
       gCCGCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTc                                                                                                                                               >  1:5414235/1‑150 (MQ=255)
         cGCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTcgc                                                                                                                                             <  2:4596213/150‑1 (MQ=255)
          gCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTcgcg                                                                                                                                            >  1:3627836/1‑150 (MQ=255)
          gCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGGGCACCTGAAGGACCTGGGCCAGCTGGTGGCTcgcg                                                                                                                                            <  2:39207/150‑1 (MQ=255)
            aaCATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGaa                                                                                                                                          >  1:4487851/1‑150 (MQ=255)
             aCATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGaaa                                                                                                                                         <  2:3774303/150‑1 (MQ=255)
                 ggCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGaaa                                                                                                                                         >  2:2124061/1‑146 (MQ=255)
                  gCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGaa                                                                                                                                    <  1:1156370/150‑1 (MQ=255)
                  gCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGaa                                                                                                                                    >  1:1851697/1‑150 (MQ=255)
                          ccGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCgg                                                                                                                            <  1:2449412/150‑1 (MQ=255)
                            gTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGgtg                                                                                                                          >  2:3939924/1‑150 (MQ=255)
                                 cTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGGGGCTCGCGAAATCGAACGCGCCGgtggcgtg                                                                                                                     <  2:3647445/150‑1 (MQ=255)
                                      gcgcgaCACCGGGATGAAGGACGAAGAATTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCaa                                                                                                                <  1:4003592/150‑1 (MQ=255)
                                        gcgcCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAgg                                                                                                              <  1:1586938/150‑1 (MQ=255)
                                        gcgcCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAgg                                                                                                              >  2:2778654/1‑150 (MQ=255)
                                          gcCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCACGGaa                                                                                                            <  2:433549/150‑1 (MQ=255)
                                          gcCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGaa                                                                                                            >  2:1859869/1‑150 (MQ=255)
                                            cACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAAtt                                                                                                          >  2:3825233/1‑150 (MQ=255)
                                             aCCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTc                                                                                                         <  1:882314/150‑1 (MQ=255)
                                               cGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGACCGCGCCGGTGGCGTGGCCAAGGAATTCaa                                                                                                       <  2:5023156/150‑1 (MQ=255)
                                               cGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCaa                                                                                                       <  2:3018301/150‑1 (MQ=255)
                                               cGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCaa                                                                                                       <  1:3476321/150‑1 (MQ=255)
                                                gggATGAAGGATGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAac                                                                                                      <  1:4383812/150‑1 (MQ=255)
                                                     gAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGaa                                                                                                                                    <  1:4567169/115‑1 (MQ=255)
                                                     gAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCAt                                                                                                 >  1:5255872/1‑150 (MQ=255)
                                                     gAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCAAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGaa                                                                                                                                    >  2:4567169/1‑115 (MQ=255)
                                                        ggACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATcgc                                                                                              <  2:3456881/150‑1 (MQ=255)
                                                           cGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGaa                                                                                                            <  1:1418590/133‑1 (MQ=255)
                                                           cGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGaa                                                                                                            >  2:1418548/1‑133 (MQ=255)
                                                           cGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTCCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGt                                                                                           <  1:4830544/150‑1 (MQ=255)
                                                            gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                          <  1:4061872/150‑1 (MQ=255)
                                                            gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                          >  2:501734/1‑150 (MQ=255)
                                                            gAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTc                                                                                          <  1:2563202/150‑1 (MQ=255)
                                                                  ttCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGTGTGGCCAAGGAATTCAACACCATCGCGGTCGATGAc                                                                                    >  2:4895411/1‑150 (MQ=255)
                                                                  ttCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGAc                                                                                    >  2:2501513/1‑150 (MQ=255)
                                                                  ttCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGAc                                                                                    >  2:1269567/1‑150 (MQ=255)
                                                                    cAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACgg                                                                                  <  2:3375758/150‑1 (MQ=255)
                                                                        aaaCCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATc                                                                              <  1:1936842/150‑1 (MQ=255)
                                                                           ccGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGcc                                                                           >  2:2371807/1‑150 (MQ=255)
                                                                            cGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCa                                                                          >  2:4968924/1‑150 (MQ=255)
                                                                             gATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCAt                                                                         >  2:4448345/1‑150 (MQ=255)
                                                                             gATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCAt                                                                         >  1:4775321/1‑150 (MQ=255)
                                                                              atcatcGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATg                                                                        <  1:3954847/150‑1 (MQ=255)
                                                                               tcatcGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATgg                                                                       <  1:617716/150‑1 (MQ=255)
                                                                               tcatcGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATgg                                                                       <  1:3344016/150‑1 (MQ=255)
                                                                                 atcgccatcgccAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGc                                                                     <  2:4648330/150‑1 (MQ=255)
                                                                                 atcgccatcgccAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGc                                                                     <  1:3076841/150‑1 (MQ=255)
                                                                                 atcgccatcgccAAATCGTTCCCCCAGTTCGTACCGGGCCATGTGCACCTGCAGGACACGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAATGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGc                                                                     <  1:746/150‑1 (MQ=255)
                                                                                        tcgccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCacgacg                                                              <  2:3166548/150‑1 (MQ=255)
                                                                                         cgccaACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACgg                                                             <  2:2708055/150‑1 (MQ=255)
                                                                                             aaCTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATg                                                         <  1:5024956/150‑1 (MQ=255)
                                                                                               cTCGTTCACCCAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCt                                                       >  1:5001597/1‑150 (MQ=255)
                                                                                               cTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCt                                                       <  2:2909851/150‑1 (MQ=255)
                                                                                               cTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCt                                                       <  1:2846652/150‑1 (MQ=255)
                                                                                                  gTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGTGTGGCCAAGGAATTc                                                                                                         <  1:2659105/97‑1 (MQ=255)
                                                                                                  gTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTa                                                    >  1:247799/1‑150 (MQ=255)
                                                                                                        ccAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCt                                              >  1:3466051/1‑150 (MQ=255)
                                                                                                        ccAGTTCGTGCCGGGCCATGTGCATCTGAAGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCt                                              >  2:3599147/1‑150 (MQ=255)
                                                                                                        ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                              <  1:5263871/150‑1 (MQ=255)
                                                                                                        ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                              <  2:1734000/150‑1 (MQ=255)
                                                                                                        ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                              >  1:3911234/1‑150 (MQ=255)
                                                                                                        ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                              >  2:3535861/1‑150 (MQ=255)
                                                                                                        ccAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCt                                              <  1:3835354/150‑1 (MQ=255)
                                                                                                         cAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTg                                             >  1:4927375/1‑150 (MQ=255)
                                                                                                            ttCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCa                                          >  2:717293/1‑150 (MQ=255)
                                                                                                                tACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGcc                                      <  1:456485/150‑1 (MQ=255)
                                                                                                                tACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGcc                                      <  1:2271229/150‑1 (MQ=255)
                                                                                                                     ggCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGaa                                 >  2:2108601/1‑150 (MQ=255)
                                                                                                                       ccATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAAt                               >  2:2832672/1‑150 (MQ=255)
                                                                                                                       ccATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAAt                               <  2:1904078/150‑1 (MQ=255)
                                                                                                                          tgtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAt                            <  1:1269567/150‑1 (MQ=255)
                                                                                                                          tgtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAt                            <  1:4448345/150‑1 (MQ=255)
                                                                                                                          tgtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAt                            >  2:2166269/1‑150 (MQ=255)
                                                                                                                           gtgCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCAtt                           <  2:4699411/150‑1 (MQ=255)
                                                                                                                             gCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGc                         <  1:4231439/150‑1 (MQ=255)
                                                                                                                                ccTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGa                      <  2:3306019/150‑1 (MQ=255)
                                                                                                                                 cTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCCTGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGAc                     >  2:1423691/1‑150 (MQ=255)
                                                                                                                                  tGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACg                    >  1:3769772/1‑150 (MQ=255)
                                                                                                                                    aaGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGcc                  >  2:1547191/1‑150 (MQ=255)
                                                                                                                                    aaGGACCTCGGCCAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGActcg                  >  2:1276178/1‑147 (MQ=255)
                                                                                                                                      ggACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGt                <  2:4178206/150‑1 (MQ=255)
                                                                                                                                       gACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTg               <  2:4511188/150‑1 (MQ=255)
                                                                                                                                       gACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGACATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTg               <  1:3302289/150‑1 (MQ=255)
                                                                                                                                        aCCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGATGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTgg              >  2:4148064/1‑150 (MQ=255)
                                                                                                                                        aCCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGATGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTgg              >  2:3688929/1‑150 (MQ=255)
                                                                                                                                           tGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAAt           >  1:3540324/1‑150 (MQ=255)
                                                                                                                                              gCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa        <  1:1005161/150‑1 (MQ=255)
                                                                                                                                              gCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa        <  2:5214760/150‑1 (MQ=255)
                                                                                                                                              gCCAGCGGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACa        >  2:2290805/1‑150 (MQ=255)
                                                                                                                                               ccAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACAt       <  1:4598443/150‑1 (MQ=255)
                                                                                                                                               ccAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACAt       <  2:584823/150‑1 (MQ=255)
                                                                                                                                               ccAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAatatat       >  1:2813624/1‑150 (MQ=255)
                                                                                                                                               ccAGCTCGTCGCGCGCGAAATTGAAGCAGCTGGCGGCGTTGCCAAGGAATTCAACACCATTGCGGTCGATGACGGTATCGCGATGGGCCATGACGGCATGCTCTATTCGCTGCCATCGCGCGAGCTGATTTCCGACTCGGTTGAatatat       >  1:1759186/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCCAAGGCCGCAACATGGCCGGCGCCCGTGCCCTGTGGCGCGCCACCGGGATGAAGGACGAAGACTTCAAGAAACCGATCATCGCCATCGCCAACTCGTTCACCCAGTTCGTACCGGGCCATGTGCACCTGAAGGACCTGGGCCAGCTGGTGGCTCGCGAAATCGAACGCGCCGGTGGCGTGGCCAAGGAATTCAACACCATCGCGGTCGATGACGGCATCGCCATGGGCCACGACGGCATGCTGTACTCGCTGCCAAGCCGCGAAATCATTGCCGACGCCGTGGAATACATGGTCA  >  NC_002947/5851301‑5851598

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: