Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,888,034 A→T 6.6% intergenic (‑43/+27) rnd ← / ← fadD ribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,888,0340AT 6.6% 301.3 / 7.5 110intergenic (‑43/+27)rnd/fadDribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)
Reads supporting (aligned to +/- strand):  ref base A (47/54);  new base T (5/2);  total (54/56)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 111 total aligned reads displayed.

GCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCG  >  NC_000913/1887889‑1888182
                                                                                                                                                  |                                                                                                                                                    
gcgTACGATCAAATTCAGTATCCAGGGCTATCGCCGGAAATGCAAGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATAATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCgg                                                                                                                                                  <  1:98163/150‑1 (MQ=255)
gcgTACGAACAAATTCAGTAGACAGGGCTATAGCCGGAAAGGCACTGACGGCTTCAAACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAACCGCCGGATTAACCgg                                                                                                                                                  <  1:102594/150‑1 (MQ=255)
  gTACGCACAAATTCAGTATCCAGGGCTAGCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCg                                                                                                                                                <  2:45778/150‑1 (MQ=255)
   tACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGACCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGt                                                                                                                                               >  2:1608460/1‑150 (MQ=255)
    aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGTACGGCTTCACACAAAGGTGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTGTTGTGGGCGTC‑AAAAAAAAAGCCGGATTAAACGGCGTc                                                                                                                                              >  2:786199/1‑150 (MQ=255)
     cGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTA‑AAAAAAAACGCCGGATTAACCGGCGTCt                                                                                                                                             >  1:1092492/1‑150 (MQ=255)
     cGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACAGACGGCTTCACACTAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCt                                                                                                                                             >  1:111080/1‑150 (MQ=255)
      gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                            <  2:707991/150‑1 (MQ=255)
      gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                            <  2:587877/150‑1 (MQ=255)
      gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                            <  2:1156267/150‑1 (MQ=255)
      gAACAAATTCAGTATCCAGGGCTAACGCCGGAAAGGCACGGACGGCGTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                            <  2:932696/150‑1 (MQ=255)
       aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTgc                                                                                                                                           >  2:1194507/1‑149 (MQ=255)
         cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTgacg                                                                                                                                         <  1:1211596/150‑1 (MQ=255)
         cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTCGTAATTCAAATTGTTTTCTCTTTAGTGGGCGGC‑AAAAAAAACGCCGGTTAATCCGGCGttttttt                                                                                                                                         <  1:103138/150‑7 (MQ=255)
         cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTGCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCCTTTGGTAATTCAAATTGTGTTCGCTTTAGGGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTgccg                                                                                                                                         >  2:932931/1‑150 (MQ=255)
         cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGCACGGCTTCACACAAAGAAGCCAGCGCATCGGCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTCGTGGGCGTC‑AAAAAAAACGCCGGAATAACCGGCGTCTgacg                                                                                                                                         <  2:1390179/150‑1 (MQ=255)
          aaaTTCACTATCCAGGGCTAGCGCCGGAACGGCACGGACGGCGTCACCCAACGAAGCCCGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAATCCGGCGtttttttt                                                                                                                                        <  2:1488271/150‑8 (MQ=255)
           aaTTCAGTTTCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAACAAGCCAGCGCAGCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGGC‑‑AAAAAAACGCCGGATTAACCGGCGTCTGACgact                                                                                                                                      <  2:1709872/150‑1 (MQ=255)
           aaTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCGGTGGGAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTgacgac                                                                                                                                       >  1:1561709/1‑150 (MQ=255)
            aTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGACGCCAGCGCATCGTCCGTGGTTATCAGTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAACACGCCGGATTAGCCGGCGTCTGACgcct                                                                                                                                      >  1:537775/1‑150 (MQ=255)
             ttCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACgactg                                                                                                                                     <  1:37031/150‑1 (MQ=255)
             ttCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACgactg                                                                                                                                     >  1:412209/1‑150 (MQ=255)
             ttCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACgactg                                                                                                                                     >  1:1463055/1‑150 (MQ=255)
               cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAGCGCCGGATTAACCGGCGTCTGACgactgac                                                                                                                                   >  1:1421728/1‑150 (MQ=255)
                aGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACgactgact                                                                                                                                  <  1:1283430/150‑1 (MQ=255)
                   aTCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTGGTTAGTGGGCTTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTaa                                                                                                                               >  2:367895/1‑150 (MQ=255)
                   aTCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTaa                                                                                                                               >  2:556787/1‑150 (MQ=255)
                             aTCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACTCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCt                                                                                                                     >  2:684127/1‑150 (MQ=255)
                                 ccGGAAAGGCACGGACTGCTTCACACCAAGAATCCAGCGCATCGTCCGTTGTAATCATTTGGTATTTAAAATTTTTTTCTCTTTAGTGGGCCTC‑AAAAATAACGCCGGTTTAACCGGCTTCTGACGACGGACTTAACGCTCAGGCTTTAt                                                                                                                 <  1:742780/150‑1 (MQ=255)
                                 ccGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCAGCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTAt                                                                                                                 <  1:751949/150‑1 (MQ=255)
                                   ggAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTg                                                                                                               <  1:1774725/150‑1 (MQ=255)
                                    gAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGt                                                                                                              <  1:1198756/150‑1 (MQ=255)
                                    gAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGt                                                                                                              <  2:1378136/150‑1 (MQ=255)
                                      aaGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTcc                                                                                                            <  2:1061529/150‑1 (MQ=255)
                                      aaGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTcc                                                                                                            <  1:574943/150‑1 (MQ=255)
                                          cACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACtt                                                                                                        >  2:780637/1‑150 (MQ=255)
                                           acggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACttt                                                                                                       <  2:260047/150‑1 (MQ=255)
                                           acggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCGTTATTGTCCACttt                                                                                                       >  2:1417132/1‑150 (MQ=255)
                                             ggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGc                                                                                                     <  2:1081883/150‑1 (MQ=255)
                                             ggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGc                                                                                                     <  2:1763439/150‑1 (MQ=255)
                                             ggacggCTTCACAAAAAGAATCCAGCGCTTGTTCCGTGTTAATCATTTTGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGc                                                                                                     <  2:1719855/150‑1 (MQ=255)
                                              gacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGc                                                                                                     >  2:1509169/1‑150 (MQ=255)
                                               acggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcg                                                                                                   <  1:1727548/150‑1 (MQ=255)
                                               acggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcg                                                                                                   <  1:325570/150‑1 (MQ=255)
                                                   cTTCACACAAAGAAGCCAGCGCATCGTCCGTGTTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                               <  1:746134/150‑1 (MQ=255)
                                                   cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGCCTTAACGCTAAGGCTTTATTGTCCACTTTGCcccgcg                                                                                               >  2:1668566/1‑150 (MQ=255)
                                                    ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCgcgcgc                                                                                              <  2:398359/150‑1 (MQ=255)
                                                       acacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCtta                                                                                           >  1:1460549/1‑149 (MQ=255)
                                                        cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGAGGGCGTC‑AAAAAAAACGCCGGATTACCCGGCGTCTGCCGACTGACTTATCGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCGTCg                                                                                          >  2:26248/1‑150 (MQ=255)
                                                         acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGt                                                                                         >  2:1273044/1‑150 (MQ=255)
                                                         acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGt                                                                                         >  1:1535926/1‑150 (MQ=255)
                                                          caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTATAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTc                                                                                        <  2:600720/150‑1 (MQ=255)
                                                           aaaGACGCCAGCGCATCGTCCGTGGTAATCCTTTGTTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAACAACGCCGGCTTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCa                                                                                       <  1:1330113/150‑1 (MQ=255)
                                                                 gCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAt                                                                                 >  1:1116924/1‑150 (MQ=255)
                                                                 gCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAAAGCTAAGGCTTTATTGTCCACTTTTCCGCCCGCTTCGTCACGTaaa                                                                                 >  1:1366468/1‑149 (MQ=255)
                                                                    aGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACTCTCAGGCTTTATTGTCCACTTTGCCGCGCTCTTCTTCACGTAAtaca                                                                              >  1:1103711/1‑147 (MQ=255)
                                                                         aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGAAGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgtcg                                                                         >  2:1408321/1‑150 (MQ=255)
                                                                               cGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCATTAGTGGGCGTC‑AAAAAAAACGCAGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGGTATATTGTCCACTTTGCCGCGCGCTGCGTCACGTTATTCTCGTCGCAAAAt                                                                   >  2:1228100/1‑150 (MQ=255)
                                                                                    tAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTc                                                              <  2:239881/150‑1 (MQ=255)
                                                                                     aaTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGCCTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTcc                                                             >  2:1013415/1‑150 (MQ=255)
                                                                                      aTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCg                                                            <  2:1463055/150‑1 (MQ=255)
                                                                                       tCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGa                                                           <  1:1619347/150‑1 (MQ=255)
                                                                                         aTTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACg                                                         <  2:207372/150‑1 (MQ=255)
                                                                                             ggTAATTCCAATTGTTTTCTCAATAGTGGGCGTC‑AAAAAAAACGCCGGATAAACCGGCGTCTGACGACTGACTTAACGCTCAGTCTTTATTGTCCACTTTGACCCGCGCTTCGTCACGTAATTCTCGTCGCAAAATATTTCCGAAGTTAg                                                     >  2:1219107/1‑150 (MQ=255)
                                                                                                aaTTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGCCGTTAGAtt                                                  >  2:826518/1‑150 (MQ=255)
                                                                                                    cAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                              >  1:580148/1‑150 (MQ=255)
                                                                                                    cAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                              <  1:1083902/150‑1 (MQ=255)
                                                                                                    aaaaTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                              <  1:1588690/149‑1 (MQ=255)
                                                                                                     aaaTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                             >  1:569033/1‑150 (MQ=255)
                                                                                                        ttGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAAc                                          <  1:1474435/150‑1 (MQ=255)
                                                                                                         tGTTTTCTCTTTAGTGGGCGTC‑AAAGAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                         >  1:1432386/1‑150 (MQ=255)
                                                                                                          gTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtc                                        <  1:1194507/150‑1 (MQ=255)
                                                                                                             ttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcatc                                     >  1:1713292/1‑150 (MQ=255)
                                                                                                             ttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcatc                                     >  2:480073/1‑150 (MQ=255)
                                                                                                               ctctTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                                   <  1:599131/150‑1 (MQ=255)
                                                                                                               ctctTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                                   >  2:512772/1‑150 (MQ=255)
                                                                                                                tctTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCAACGTTAGATTTCGGTAACTCATCACg                                  <  2:111080/150‑1 (MQ=255)
                                                                                                                  tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGAAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaa                                >  1:1262786/1‑150 (MQ=255)
                                                                                                                     aGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAAc                              <  1:1122245/150‑1 (MQ=255)
                                                                                                                       tGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTcc                           <  2:210369/150‑1 (MQ=255)
                                                                                                                        gggCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGAGTTAACGCGCAGGCTTTTTTGTCCACTTTGACGCGCGCTTCGTCACGAAATTCTCGTCGCAAAATTTTTCCGACGGTAGATTTCGGTAAGTCATCACGAACCTccc                          >  2:77437/1‑149 (MQ=255)
                                                                                                                         ggCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccac                         >  1:440294/1‑150 (MQ=255)
                                                                                                                          gCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                        <  1:451260/150‑1 (MQ=255)
                                                                                                                           cGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTAGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacca                       <  1:1205429/150‑1 (MQ=255)
                                                                                                                            gTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAg                      >  1:1844802/1‑150 (MQ=255)
                                                                                                                                 cgtcagaCGCCGGTTAATCCGGCGTCTGCCGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTGGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc                  >  1:122331/7‑150 (MQ=255)
                                                                                                                                  aaaaaaCGCCGGATTTACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCg                 <  1:1729113/150‑1 (MQ=255)
                                                                                                                                   aaaaaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCgg                <  2:1418780/150‑1 (MQ=255)
                                                                                                                                     agaCGCCGGTTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTGTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTGAGATTTCGGTACCTCATCACGCAACTCCACCAGCTTCGGTa              >  1:1039637/3‑150 (MQ=255)
                                                                                                                                      gaCGCCGGTTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTaa             >  1:1075722/2‑149 (MQ=255)
                                                                                                                                        cGCCGGTTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACtt           >  2:43130/1‑150 (MQ=255)
                                                                                                                                          ccGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTg         <  2:580334/150‑1 (MQ=255)
                                                                                                                                          ccGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTg         >  1:865602/1‑150 (MQ=255)
                                                                                                                                           cGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGt        >  2:1294087/1‑150 (MQ=255)
                                                                                                                                             gATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCCCGTAATTCTCGTCTCAAAATTTTTCAAACGTTAGATTTCGGTAACTCATCACGAAACTCCAACAGATTCGGTACTTTGTAt      >  2:1420966/1‑150 (MQ=255)
                                                                                                                                              aTTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATc     <  1:684127/150‑1 (MQ=255)
                                                                                                                                              aTTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATc     <  1:1509169/150‑1 (MQ=255)
                                                                                                                                               ttAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATcc    <  1:556787/150‑1 (MQ=255)
                                                                                                                                                tAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATccc   <  1:171913/150‑1 (MQ=255)
                                                                                                                                                 aaCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACGCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCg  <  1:1871705/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCG  >  NC_000913/1887889‑1888182

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑