Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,888,034 A→T 7.1% 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 7.1% 588.4 / 26.3 168intergenic (‑43/+27)rnd/fadDribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)
Reads supporting (aligned to +/- strand):  ref base A (80/76);  new base T (7/5);  total (87/81)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 175 total aligned reads displayed.

TGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCG  >  NC_000913/1887886‑1888182
                                                                                                                                                    |                                                                                                                                                     
   gcgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCCAATTGTTTTCTCTTTAGTGGGCGTCAAAAAACACGCCGGATTAA‑CCg                                                                                                                                                   <  1:1158158/149‑1 (MQ=255)
       aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                               <  1:1320995/149‑1 (MQ=255)
         gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCt                                                                                                                                             >  2:1106419/1‑149 (MQ=255)
          aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAAT‑CCGGCGttt                                                                                                                                             <  1:662541/148‑3 (MQ=255)
          aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTg                                                                                                                                            >  3:1311611/1‑149 (MQ=255)
          aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTg                                                                                                                                            >  2:567652/1‑149 (MQ=255)
           aCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTaa                                                                                                                                                       <  2:713059/138‑1 (MQ=255)
           aCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTga                                                                                                                                           <  3:989581/149‑1 (MQ=255)
            cAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAAT‑CCGGCGtttttt                                                                                                                                          <  3:328846/149‑6 (MQ=255)
             aaaTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTgacg                                                                                                                                         >  2:119685/1‑149 (MQ=255)
              aaTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGG‑TTAATCCGGCGtttttttt                                                                                                                                        <  2:1629486/149‑8 (MQ=255)
                ttCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgact                                                                                                                                      >  3:833515/1‑149 (MQ=255)
                  cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgactga                                                                                                                                    <  3:1467422/149‑1 (MQ=255)
                  cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgactga                                                                                                                                    >  3:268778/1‑149 (MQ=255)
                       tCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTaa                                                                                                                               <  2:1295640/149‑1 (MQ=255)
                         cAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCg                                                                                                                                                   <  4:496305/127‑1 (MQ=255)
                         cAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCg                                                                                                                                                   >  3:496305/1‑127 (MQ=255)
                          aGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGc                                                                                                                            >  3:1393435/1‑149 (MQ=255)
                             gCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCa                                                                                                                         <  2:429727/149‑1 (MQ=255)
                                   gCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCttt                                                                                                                   >  3:1025733/1‑149 (MQ=255)
                                   gCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTACGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCttt                                                                                                                   <  1:502834/149‑1 (MQ=255)
                                    ccGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTa                                                                                                                  <  4:183470/149‑1 (MQ=255)
                                      ggAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTAtt                                                                                                                >  2:1648050/1‑149 (MQ=255)
                                              acggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACtt                                                                                                        >  4:722907/1‑149 (MQ=255)
                                                ggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTg                                                                                                      >  4:1520234/1‑149 (MQ=255)
                                                  acggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGAAATTCAAATTGTATTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGcc                                                                                                    >  4:1606983/1‑149 (MQ=255)
                                                  acggCTTCACAAAAAGAAGCCAGCGCATCGTCCGTAGTAATCATTTGGTAATTCAAATTGTTTTCGCTTTAATGGCCGACAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACAACTGACGTAACGCTCAGGCTTTATTGTCCACTTTGcc                                                                                                    <  4:3233/149‑1 (MQ=255)
                                                    ggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgc                                                                                                  >  3:854210/1‑149 (MQ=255)
                                                      cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑c                                                                                                                                                     >  1:971068/1‑96 (MQ=255)
                                                      cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑c                                                                                                                                                     <  2:971068/96‑1 (MQ=255)
                                                       ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                               >  4:1195972/1‑149 (MQ=255)
                                                       ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                               >  1:1002753/1‑149 (MQ=255)
                                                         cacacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCt                                                                                             <  3:1520234/149‑1 (MQ=255)
                                                         cacacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGAGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCt                                                                                                                     >  1:306089/1‑125 (MQ=255)
                                                           cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTc                                                                                           <  1:1635642/149‑1 (MQ=255)
                                                           cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTc                                                                                           <  3:1225403/149‑1 (MQ=255)
                                                            acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                          >  2:1209719/1‑149 (MQ=255)
                                                              aaaGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTc                                                                                        >  4:38537/1‑149 (MQ=255)
                                                                agaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                      >  4:418549/1‑149 (MQ=255)
                                                                 gaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACg                                                                                     >  2:1375171/1‑149 (MQ=255)
                                                                  aagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGt                                                                                    >  3:1259330/1‑149 (MQ=255)
                                                                    gCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTaa                                                                                  >  4:873593/1‑149 (MQ=255)
                                                                     ccAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑c                                                                                                                                                     <  3:1264663/81‑1 (MQ=255)
                                                                     ccAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAt                                                                                 >  3:990431/1‑149 (MQ=255)
                                                                      cAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAtt                                                                                >  2:792437/1‑149 (MQ=255)
                                                                      cAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAtt                                                                                <  1:1281003/149‑1 (MQ=255)
                                                                           cATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgt                                                                           <  4:1613301/149‑1 (MQ=255)
                                                                            aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                               >  1:73661/1‑80 (MQ=255)
                                                                            aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                               <  2:73661/80‑1 (MQ=255)
                                                                                   gTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAAT‑CCGGCGtt                                                                                                                                              <  2:1573054/74‑2 (MQ=39)
                                                                                   gTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAt                                                                   <  3:418549/149‑1 (MQ=255)
                                                                                      gTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgact                                                                                                                                      <  3:1625483/79‑1 (MQ=255)
                                                                                      gTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgact                                                                                                                                      >  4:1625483/1‑79 (MQ=255)
                                                                                       tAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                               <  4:119831/149‑1 (MQ=255)
                                                                                         aTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTcc                                                             <  4:854210/149‑1 (MQ=255)
                                                                                          tCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCg                                                            >  1:4230/1‑149 (MQ=255)
                                                                                             tttGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACg                                                         >  4:689058/1‑149 (MQ=255)
                                                                                                ggTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTa                                                      <  4:1382927/149‑1 (MQ=255)
                                                                                                  tAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGa                                                    >  4:165101/1‑149 (MQ=255)
                                                                                                    aTTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAtt                                                  <  4:1311611/149‑1 (MQ=255)
                                                                                                     ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAttt                                                 >  3:936001/1‑149 (MQ=255)
                                                                                                      tCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAACAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTc                                                <  4:526291/149‑1 (MQ=255)
                                                                                                          aTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGc                                                                                                                                                 >  1:671788/1‑48 (MQ=255)
                                                                                                          aTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTaa                                                                                  >  4:586177/1‑111 (MQ=255)
                                                                                                          aTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTaa                                                                                  <  3:586177/111‑1 (MQ=255)
                                                                                                          aTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTa                                            >  1:1016004/1‑149 (MQ=255)
                                                                                                             gTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                         <  3:1606983/149‑1 (MQ=255)
                                                                                                                ttCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcat                                      >  1:812204/1‑149 (MQ=255)
                                                                                                                ttCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcat                                      <  4:268778/149‑1 (MQ=255)
                                                                                                                 tctcTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcatc                                     >  1:293685/1‑149 (MQ=255)
                                                                                                                   tctTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                                   >  4:96072/1‑149 (MQ=255)
                                                                                                                     tttAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                      >  4:1292004/1‑96 (MQ=255)
                                                                                                                     tttAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCg                                               >  1:1146169/1‑135 (MQ=255)
                                                                                                                     tttAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGa                                 >  4:424517/1‑149 (MQ=255)
                                                                                                                     tttAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGa                                 >  2:230095/1‑149 (MQ=255)
                                                                                                                              cGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                        <  3:722907/149‑1 (MQ=255)
                                                                                                                                tcgaaaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCATAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAg                      >  4:555243/4‑149 (MQ=255)
                                                                                                                                tCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAg                      >  4:1057415/1‑149 (MQ=255)
                                                                                                                                  aaaaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCt                    >  2:666086/1‑149 (MQ=255)
                                                                                                                                   aaaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccac                         >  1:925692/1‑143 (MQ=255)
                                                                                                                                   aaaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccac                         <  2:925692/143‑1 (MQ=255)
                                                                                                                                    aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                         >  1:1134923/1‑126 (MQ=255)
                                                                                                                                    aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                         <  2:1134923/126‑1 (MQ=255)
                                                                                                                                    aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc                  <  2:410940/149‑1 (MQ=255)
                                                                                                                                    aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc                  >  4:1118385/1‑149 (MQ=255)
                                                                                                                                    aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCACAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc                  >  1:1129374/1‑149 (MQ=255)
                                                                                                                                       agaCGCCGGTTAAT‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGt               >  3:921732/3‑149 (MQ=255)
                                                                                                                                        gaCGCCGGTTAAT‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa              >  3:1671936/2‑149 (MQ=255)
                                                                                                                                        gaCGCCGGTTAAT‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa              >  2:607385/2‑149 (MQ=255)
                                                                                                                                         aCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTAc             <  2:1705873/149‑1 (MQ=255)
                                                                                                                                           gCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACtt           <  1:1462467/149‑1 (MQ=255)
                                                                                                                                            ccGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACttt          <  2:812204/149‑1 (MQ=255)
                                                                                                                                            ccGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACACATCACGAAACTCCACCAGCTTCGGTACttt          <  1:1502710/149‑1 (MQ=255)
                                                                                                                                             cGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                               <  3:494478/95‑1 (MQ=255)
                                                                                                                                              ggATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGt        >  2:763097/1‑149 (MQ=255)
                                                                                                                                               gcttAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTa       >  2:1720622/3‑149 (MQ=255)
                                                                                                                                                   aa‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATccc   >  1:1434357/1‑149 (MQ=255)
                                                                                                                                                   aa‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATccc   <  4:1440072/149‑1 (MQ=255)
                                                                                                                                                   aa‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATccc   <  4:930106/149‑1 (MQ=255)
                                                                                                                                                    a‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCg  <  4:1241353/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCG  >  NC_000913/1887886‑1888182

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑