Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,888,032 T→A 5.8% intergenic (‑41/+29) 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,0320TA 5.8% 333.4 / 5.4 111intergenic (‑41/+29)rnd/fadDribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)
Reads supporting (aligned to +/- strand):  ref base T (51/53);  new base A (5/2);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

AGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCC  >  NC_000913/1887884‑1888181
                                                                                                                                                     |                                                                                                                                                     
tgTGCGCGTACGAAGAAATTCAGTATCCCGGGCTATGGCCGGAAAGGCACGGACGGCTTCACACATCGTAGCCAGCGCATCGTCCGTGGTAATCTTTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTa                                                                                                                                                      <  1:1665714/149‑1 (MQ=255)
aGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTa                                                                                                                                                      >  1:398359/1‑150 (MQ=255)
aGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTa                                                                                                                                                      >  1:580334/1‑150 (MQ=255)
     gcgTACGATCAAATTCAGTATCCAGGGCTATCGCCGGAAATGCAAGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATAATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCgg                                                                                                                                                 <  1:98163/150‑1 (MQ=255)
     gcgTACGAACAAATTCAGTAGACAGGGCTATAGCCGGAAAGGCACTGACGGCTTCAAACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAACCGCCGGATTAACCgg                                                                                                                                                 <  1:102594/150‑1 (MQ=255)
       gTACGCACAAATTCAGTATCCAGGGCTAGCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCg                                                                                                                                               <  2:45778/150‑1 (MQ=255)
        tACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTATCCGGCGt                                                                                                                                              >  1:160953/1‑150 (MQ=255)
        tACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGACCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGt                                                                                                                                              >  2:1608460/1‑150 (MQ=255)
         aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGTACGGCTTCACACAAAGGTGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTGTTGTGGGCGTCAAAAAAAAAGCCGGATTAAACGGCGTc                                                                                                                                             >  2:786199/1‑150 (MQ=255)
          cGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTAAAAAAAAACGCCGGATTAACCGGCGTCt                                                                                                                                            >  1:1092492/1‑150 (MQ=255)
          cGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACAGACGGCTTCACACTAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCt                                                                                                                                            >  1:111080/1‑150 (MQ=255)
           gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                           <  2:587877/150‑1 (MQ=255)
           gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                           <  2:707991/150‑1 (MQ=255)
           gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                           <  2:1156267/150‑1 (MQ=255)
           gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACAAAAAGAAGCCAGCGCATCGTCGGTGGTAAACATTTTGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                           <  1:1146507/150‑1 (MQ=255)
           gAACAAATTCAGTATCCAGGGCTAACGCCGGAAAGGCACGGACGGCGTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                           <  2:932696/150‑1 (MQ=255)
            aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTgc                                                                                                                                          >  2:1194507/1‑149 (MQ=255)
              cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTgacg                                                                                                                                        <  1:1211596/150‑1 (MQ=255)
              cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTCGTAATTCAAATTGTTTTCTCTTTAGTGGGCGGCAAAAAAAACGCCGGTTAATCCGGCGttttttt                                                                                                                                        <  1:103138/150‑7 (MQ=255)
              cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTGCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCCTTTGGTAATTCAAATTGTGTTCGCTTTAGGGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTgccg                                                                                                                                        >  2:932931/1‑150 (MQ=255)
              cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGCACGGCTTCACACAAAGAAGCCAGCGCATCGGCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTCGTGGGCGTCAAAAAAAACGCCGGAATAACCGGCGTCTgacg                                                                                                                                        <  2:1390179/150‑1 (MQ=255)
               aaaTTCACTATCCAGGGCTAGCGCCGGAACGGCACGGACGGCGTCACCCAACGAAGCCCGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAATCCGGCGtttttttt                                                                                                                                       <  2:1488271/150‑8 (MQ=255)
                aaTTCAGTTTCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAACAAGCCAGCGCAGCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGGCAAAAAAACGCCGGATTAACCGGCGTCTGACgact                                                                                                                                     <  2:1709872/150‑1 (MQ=255)
                aaTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCGGTGGGAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTgacgac                                                                                                                                      >  1:1561709/1‑150 (MQ=255)
                 aTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGACGCCAGCGCATCGTCCGTGGTTATCAGTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAACACGCCGGATTAGCCGGCGTCTGACgcct                                                                                                                                     >  1:537775/1‑150 (MQ=255)
                  ttCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactg                                                                                                                                    >  1:1463055/1‑150 (MQ=255)
                  ttCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactg                                                                                                                                    <  1:37031/150‑1 (MQ=255)
                    cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAGCGCCGGATTAACCGGCGTCTGACgactgac                                                                                                                                  >  1:1421728/1‑150 (MQ=255)
                        aTCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTGGTTAGTGGGCTTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTaa                                                                                                                              >  2:367895/1‑150 (MQ=255)
                        aTCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTaa                                                                                                                              >  2:556787/1‑150 (MQ=255)
                                      ccGGAAAGGCACGGACTGCTTCACACCAAGAATCCAGCGCATCGTCCGTTGTAATCATTTGGTATTTAAAATTTTTTTCTCTTTAGTGGGCCTCAAAAATAACGCCGGTTTAACCGGCTTCTGACGACGGACTTAACGCTCAGGCTTTAt                                                                                                                <  1:742780/150‑1 (MQ=255)
                                      ccGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCAGCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTAt                                                                                                                <  1:751949/150‑1 (MQ=255)
                                        ggAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTg                                                                                                              <  1:1774725/150‑1 (MQ=255)
                                         gAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGt                                                                                                             <  1:1198756/150‑1 (MQ=255)
                                         gAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGt                                                                                                             <  2:1378136/150‑1 (MQ=255)
                                           aaGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTcc                                                                                                           <  2:1061529/150‑1 (MQ=255)
                                           aaGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTcc                                                                                                           <  1:574943/150‑1 (MQ=255)
                                            aGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCa                                                                                                          >  2:597462/1‑150 (MQ=255)
                                               cACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACtt                                                                                                       >  2:780637/1‑150 (MQ=255)
                                                acggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACttt                                                                                                      <  2:260047/150‑1 (MQ=255)
                                                acggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCGTTATTGTCCACttt                                                                                                      >  2:1417132/1‑150 (MQ=255)
                                                  ggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGc                                                                                                    <  2:1081883/150‑1 (MQ=255)
                                                  ggacggCTTCACAAAAAGAATCCAGCGCTTGTTCCGTGTTAATCATTTTGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGc                                                                                                    <  2:1719855/150‑1 (MQ=255)
                                                   gacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGc                                                                                                    >  2:1509169/1‑150 (MQ=255)
                                                    acggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcg                                                                                                  <  1:325570/150‑1 (MQ=255)
                                                    acggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcg                                                                                                  <  1:1727548/150‑1 (MQ=255)
                                                     cggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgc                                                                                                 >  2:1661535/1‑150 (MQ=255)
                                                        cTTCACACAAAGAAGCCAGCGCATCGTCCGTGTTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                              <  1:746134/150‑1 (MQ=255)
                                                        cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGCCTTAACGCTAAGGCTTTATTGTCCACTTTGCcccgcg                                                                                              >  2:1668566/1‑150 (MQ=255)
                                                         ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCgcgcgc                                                                                             <  2:398359/150‑1 (MQ=255)
                                                            acacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCtta                                                                                          >  1:1460549/1‑149 (MQ=255)
                                                             cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGAGGGCGTCAAAAAAAACGCCGGATTACCCGGCGTCTGCCGACTGACTTATCGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCGTCg                                                                                         >  2:26248/1‑150 (MQ=255)
                                                              acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGt                                                                                        >  1:1535926/1‑150 (MQ=255)
                                                               caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTATAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTc                                                                                       <  2:600720/150‑1 (MQ=255)
                                                                aaaGACGCCAGCGCATCGTCCGTGGTAATCCTTTGTTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAACAACGCCGGCTTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCa                                                                                      <  1:1330113/150‑1 (MQ=255)
                                                                      gCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAAAGCTAAGGCTTTATTGTCCACTTTTCCGCCCGCTTCGTCACGTaaa                                                                                >  1:1366468/1‑149 (MQ=255)
                                                                         aGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACTCTCAGGCTTTATTGTCCACTTTGCCGCGCTCTTCTTCACGTAAtaca                                                                             >  1:1103711/1‑147 (MQ=255)
                                                                           ggcATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcg                                                                           <  1:1401958/149‑1 (MQ=255)
                                                                              aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGAAGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgtcg                                                                        >  2:1408321/1‑150 (MQ=255)
                                                                                    cGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCATTAGTGGGCGTCAAAAAAAACGCAGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGGTATATTGTCCACTTTGCCGCGCGCTGCGTCACGTTATTCTCGTCGCAAAAt                                                                  >  2:1228100/1‑150 (MQ=255)
                                                                                         tAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTc                                                             <  2:239881/150‑1 (MQ=255)
                                                                                            tCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGa                                                          >  2:1061916/1‑150 (MQ=255)
                                                                                            tCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGa                                                          <  1:1619347/150‑1 (MQ=255)
                                                                                              aTTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACg                                                        <  2:207372/150‑1 (MQ=255)
                                                                                                  ggTAATTCCAATTGTTTTCTCAATAGTGGGCGTCAAAAAAAACGCCGGATAAACCGGCGTCTGACGACTGACTTAACGCTCAGTCTTTATTGTCCACTTTGACCCGCGCTTCGTCACGTAATTCTCGTCGCAAAATATTTCCGAAGTTAg                                                    >  2:1219107/1‑150 (MQ=255)
                                                                                                     aaTTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGCCGTTAGAtt                                                 >  2:826518/1‑150 (MQ=255)
                                                                                                         cAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                             <  1:1083902/150‑1 (MQ=255)
                                                                                                         cAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                             >  1:580148/1‑150 (MQ=255)
                                                                                                         aaaaTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                             <  1:1588690/149‑1 (MQ=255)
                                                                                                          aaaTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGCCTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGGGCTTCGTCACGTAAGTCTCGGCGCAAAATTTTTGCGACGTTAGATTTCGGt                                            >  1:351783/1‑150 (MQ=255)
                                                                                                          aaaTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                            >  1:569033/1‑150 (MQ=255)
                                                                                                             ttGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAAc                                         <  1:1474435/150‑1 (MQ=255)
                                                                                                              tGTTTTCTCTTTAGTGGGCGTCAAAGAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                        >  1:1432386/1‑150 (MQ=255)
                                                                                                               gTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtc                                       <  1:1194507/150‑1 (MQ=255)
                                                                                                                  ttCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcatc                                    >  1:1713292/1‑150 (MQ=255)
                                                                                                                  ttCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcatc                                    >  2:480073/1‑150 (MQ=255)
                                                                                                                    ctctTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                                  >  2:512772/1‑150 (MQ=255)
                                                                                                                    ctctTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                                  <  1:599131/150‑1 (MQ=255)
                                                                                                                     tctTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCAACGTTAGATTTCGGTAACTCATCACg                                 <  2:111080/150‑1 (MQ=255)
                                                                                                                       tttAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGAAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaa                               >  1:1262786/1‑150 (MQ=255)
                                                                                                                          aGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAAc                             <  1:1122245/150‑1 (MQ=255)
                                                                                                                            tGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTcc                          <  2:210369/150‑1 (MQ=255)
                                                                                                                             gggCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGAGTTAACGCGCAGGCTTTTTTGTCCACTTTGACGCGCGCTTCGTCACGAAATTCTCGTCGCAAAATTTTTCCGACGGTAGATTTCGGTAAGTCATCACGAACCTccc                         >  2:77437/1‑149 (MQ=255)
                                                                                                                              ggCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccac                        >  1:440294/1‑150 (MQ=255)
                                                                                                                               gCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                       <  1:451260/150‑1 (MQ=255)
                                                                                                                                cGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTAGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacca                      <  1:1205429/150‑1 (MQ=255)
                                                                                                                                 gTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAg                     >  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)
                                                                                                                                          aaaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa             <  2:1380055/150‑1 (MQ=255)
                                                                                                                                           gaCGCCGGTTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTaa            >  1:1075722/2‑149 (MQ=255)
                                                                                                                                               ccGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTg        <  2:580334/150‑1 (MQ=255)
                                                                                                                                               ccGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTg        >  1:865602/1‑150 (MQ=255)
                                                                                                                                                cGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGt       >  2:1294087/1‑150 (MQ=255)
                                                                                                                                                 ggATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTa      <  2:268007/150‑1 (MQ=255)
                                                                                                                                                  gATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCCCGTAATTCTCGTCTCAAAATTTTTCAAACGTTAGATTTCGGTAACTCATCACGAAACTCCAACAGATTCGGTACTTTGTAt     >  2:1420966/1‑150 (MQ=255)
                                                                                                                                                   aTTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATc    <  1:684127/150‑1 (MQ=255)
                                                                                                                                                    ttAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATcc   <  1:556787/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCC  >  NC_000913/1887884‑1888181

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: