Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,353,340 C→G 8.0% P141R (CCG→CGG)  yrpB → 2‑nitropropane dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,353,3400CG 8.0% 351.2 / 19.1 114P141R (CCG→CGG) yrpB2‑nitropropane dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (51/53);  new base G (4/5);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 118 total aligned reads displayed.

GAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATGT  >  NC_002947/4353191‑4353489
                                                                                                                                                       |                                                                                                                                                       
gAAGCAGGCACTTGCGCTTTATTACAGGGAAGTGGGCGCCGTTTTCACGGCCCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGcc                                                                                                                                                         <  2:2102652/150‑1 (MQ=255)
gAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGcc                                                                                                                                                         >  2:2266187/1‑150 (MQ=255)
    cAGGCACTGGGCCTTTATTACGGCAAGGTGGGGGCCCATTTCACGGCCCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAg                                                                                                                                                     <  1:2935181/150‑1 (MQ=255)
      ggCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCGGCAGGc                                                                                                                                                   >  1:2143441/1‑150 (MQ=255)
          cTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGaa                                                                                                                                               <  2:2414122/150‑1 (MQ=255)
          cTTCAGCCTTATTACAGCGAAGTGGGCGCCGATTTAACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGaa                                                                                                                                               <  2:734955/150‑1 (MQ=255)
           ttGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCGGCAGGCCGaaa                                                                                                                                              <  2:486042/150‑2 (MQ=255)
           ttGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAAc                                                                                                                                              >  2:222636/1‑150 (MQ=255)
            tGACCCTTCTCACAGCGCAGGCGGCGGCAGTTCCACGGCGCAAACACGGGTATCCAACCTTGCGCCTGTCGACGAGCAGAGCTGCCTGCTGGTCGACGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCGGCAGGCCGaaat                                                                                                                                             <  1:2636237/150‑3 (MQ=255)
             gAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAActg                                                                                                                                            <  2:1449383/150‑1 (MQ=255)
              aGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGCAGACGTTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCCGCCAGCCGCAGGCCGAActgc                                                                                                                                           >  1:1725457/1‑150 (MQ=255)
               gCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAActgct                                                                                                                                          <  1:984059/150‑1 (MQ=255)
                ccTTATTACAGCGAGGGGGGCGCCAATTCCCCGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCGGCAGGCCGaaatggaa                                                                                                                                         <  2:2143441/150‑7 (MQ=255)
                 cTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACtgctgc                                                                                                                                        >  2:2397984/1‑150 (MQ=255)
                     ttaCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCg                                                                                                                                    >  2:1494823/1‑150 (MQ=255)
                     ttaCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCg                                                                                                                                    >  2:1121828/1‑150 (MQ=255)
                      taCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                   <  2:339960/150‑1 (MQ=255)
                      taCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                   <  1:2266187/150‑1 (MQ=255)
                        cAGCGAAGTGGGCGCCGATTTCACGGGTCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGt                                                                                                                                 >  1:415227/1‑150 (MQ=255)
                         aGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGGGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTa                                                                                                                                >  2:2830017/1‑150 (MQ=255)
                            gAAGTGGGCGCTGACTGCACGGCGCCAACGCCTGTGTCCAATCGCGCCCCTTTCGGCGAGCGGAGCTGCCTGCTGGTCGAAGCGTGGCGCCCGGAGGTGGTAAGCTTCCATTTCGGCCTGCCGCAGGACTCACTGCGGCAGCGGCTaaac                                                                                                                             <  2:613490/150‑2 (MQ=255)
                            gAAGGGGGCGACGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAg                                                                                                                             <  2:1105276/150‑1 (MQ=255)
                             gagggggCGCCCGTTTCCGCGCGCCCGACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGTTATAg                                                                                                                             <  2:3061575/146‑1 (MQ=255)
                                    cgcCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCgg                                                                                                                     >  1:653317/1‑150 (MQ=255)
                                     gcCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTATCACCAACAAGGGGTCGCTTGTCTTGCAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGt                                                                                                                    <  1:1135061/150‑1 (MQ=255)
                                          tttCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                   <  2:181786/130‑1 (MQ=255)
                                          tttCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                   >  1:181786/1‑130 (MQ=255)
                                           ttCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAg                                                                                                              <  1:1121828/150‑1 (MQ=255)
                                              tcGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                          <  2:1813087/149‑1 (MQ=255)
                                              aCGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                          <  2:1539733/150‑1 (MQ=255)
                                               cGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTAc                                                                                                         <  1:494337/150‑1 (MQ=255)
                                                    cccACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTcc                                                                                                    <  2:1783828/150‑1 (MQ=255)
                                                        aCGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAgc                                                                                                 >  1:737196/1‑150 (MQ=255)
                                                        aCGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAgc                                                                                                 <  2:2917709/150‑1 (MQ=255)
                                                                    aaCCGTGCGCCTTTCGACGAGCAGAGCTGTCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                    >  1:2800148/1‑150 (MQ=255)
                                                                    aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                    >  2:1075937/1‑150 (MQ=255)
                                                                    aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                    >  2:2083174/1‑150 (MQ=255)
                                                                    aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                    >  2:2716047/1‑150 (MQ=255)
                                                                     aCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTgg                                                                                   <  1:1075937/150‑1 (MQ=255)
                                                                       cGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaa                                                                                 <  2:1183474/150‑1 (MQ=255)
                                                                         tGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaaga                                                                               >  2:1685828/1‑150 (MQ=255)
                                                                          gcgcCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCATGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaagaa                                                                              >  1:2422134/1‑150 (MQ=255)
                                                                            gcCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGc                                                                            <  1:222636/150‑1 (MQ=255)
                                                                             ccTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCg                                                                           <  1:1494823/150‑1 (MQ=255)
                                                                              cTTTCGACTAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCgg                                                                          >  1:2432886/1‑150 (MQ=255)
                                                                               tttCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGt                                                                                     <  2:2004251/138‑1 (MQ=255)
                                                                               tttCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGc                                                                         <  1:768162/150‑1 (MQ=255)
                                                                                     cgaGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctggct                                                                   <  2:415227/150‑1 (MQ=255)
                                                                                      gaGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctggctg                                                                  <  1:2716047/150‑1 (MQ=255)
                                                                                        gCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCaa                                                                                                               >  1:1190588/1‑105 (MQ=255)
                                                                                        gCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCaa                                                                                                               <  2:1190588/105‑1 (MQ=255)
                                                                                        gCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGa                                                                <  1:2324654/150‑1 (MQ=255)
                                                                                             agCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAgcgc                                                            >  2:1787458/1‑150 (MQ=255)
                                                                                                tgcctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAAt                                                         <  1:370208/150‑1 (MQ=255)
                                                                                                   ctgctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTTGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCCCGACCGTGGAAGCAGCGGCCTGGCTCGAGAGCactcgc                                                      >  1:2033763/1‑145 (MQ=255)
                                                                                                    tgctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                          >  2:15324/1‑98 (MQ=255)
                                                                                                    tgctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                          <  1:15324/98‑1 (MQ=255)
                                                                                                     gctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtg                                                    >  1:1336605/1‑150 (MQ=255)
                                                                                                         gTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgat                                                <  1:2397984/150‑1 (MQ=255)
                                                                                                          tCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatg                                               <  2:2782125/150‑1 (MQ=255)
                                                                                                           cGAAGCGTGGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCATTGGCtgcgatgg                                              <  2:1725457/150‑2 (MQ=255)
                                                                                                            gAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCCGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatgcg                                             >  1:916804/1‑150 (MQ=255)
                                                                                                               gCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGatc                                          >  1:2588872/1‑150 (MQ=255)
                                                                                                                        ccGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGa                                                                <  1:1840776/119‑1 (MQ=255)
                                                                                                                        ccGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGa                                                                >  2:1840776/1‑119 (MQ=255)
                                                                                                                         cGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGTCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATgg                                >  2:582591/1‑150 (MQ=255)
                                                                                                                          ggAGGTGGTGAGCTTCCATTTCGGCCTGCTTTAGTCCTAACCGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCCGCCTGGCTGGAGCGCCATGGGTGCGACGCGATCATCGCCATggg                               <  2:805622/150‑1 (MQ=255)
                                                                                                                          ggAGGTGGTGAGCTTCCATTTCGGCCTGCCGCATGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGCTGCGATCTTCGCCTTggg                               <  2:2422134/150‑1 (MQ=255)
                                                                                                                            aGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGt                             <  1:1787458/150‑1 (MQ=255)
                                                                                                                              gtggtgAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTAt                           >  1:2241920/1‑150 (MQ=255)
                                                                                                                                 gtgAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGc                                                                         >  2:2778232/1‑101 (MQ=255)
                                                                                                                                 gtgAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGc                                                                         <  1:2778232/101‑1 (MQ=255)
                                                                                                                                 gtgAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGaa                        <  2:653317/150‑1 (MQ=255)
                                                                                                                                   gAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctggc                                                                    >  1:1317136/1‑104 (MQ=255)
                                                                                                                                   gAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctggc                                                                    <  2:1317136/104‑1 (MQ=255)
                                                                                                                                   gAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGc                      >  1:670959/1‑150 (MQ=255)
                                                                                                                                     gCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCg                    <  1:2083174/150‑1 (MQ=255)
                                                                                                                                      cTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTCAAGGCCAGCGGGGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCgg                   <  2:2756170/150‑1 (MQ=255)
                                                                                                                                      cTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaa                                                                                 >  1:1127050/1‑88 (MQ=255)
                                                                                                                                      cTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCgg                   >  1:1421378/1‑150 (MQ=255)
                                                                                                                                       ttCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaaga                                                                               <  2:226328/89‑1 (MQ=255)
                                                                                                                                       ttCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGt                  >  2:422887/1‑150 (MQ=255)
                                                                                                                                       ttCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGt                  <  2:1243660/150‑1 (MQ=255)
                                                                                                                                       ttCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGt                  >  2:388564/1‑150 (MQ=255)
                                                                                                                                        tCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                 <  2:2588872/150‑1 (MQ=255)
                                                                                                                                        tCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                 >  2:1796232/1‑150 (MQ=255)
                                                                                                                                        agcagttCGGCCTGCGGCAGGCCGAAATGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                 >  1:2957339/6‑150 (MQ=255)
                                                                                                                                         ccATTTCGGCCTGGCGCAGGCCGAACTGCTGCAGCGGGTAAAAGCCCGCGGTGCCCAGGTTATTGTCAGCGCCCACGCCCGGGAGAGAGCAGCCTGGCTTGAGGGCCATGGGCTCGTTGCGATCAACGCCCTGGGGGATGAAGACCGTgg                >  2:566666/1‑150 (MQ=255)
                                                                                                                                          cATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTTCGATGCTATCATCGCCATGGGGGATGAAGCCGCTGGc               >  1:1120513/1‑150 (MQ=255)
                                                                                                                                            tttCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGc                      >  2:2054550/1‑141 (MQ=255)
                                                                                                                                            tttCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGc                      <  1:2054550/141‑1 (MQ=255)
                                                                                                                                            tttCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCa             >  2:3071609/1‑150 (MQ=255)
                                                                                                                                            tttCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCa             >  2:1091742/1‑150 (MQ=255)
                                                                                                                                            tttCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCa             >  1:918974/1‑150 (MQ=255)
                                                                                                                                            tttCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCa             >  1:1289542/1‑150 (MQ=255)
                                                                                                                                             ttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGa                                                                >  1:1786469/1‑98 (MQ=255)
                                                                                                                                             ttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgc                                                   <  2:2748861/111‑1 (MQ=255)
                                                                                                                                             ttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgc                                                   >  1:2748861/1‑111 (MQ=255)
                                                                                                                                                     gCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCAt    <  2:1421378/150‑1 (MQ=255)
                                                                                                                                                     gCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCAt    <  2:918974/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
GAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATGT  >  NC_002947/4353191‑4353489

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: