Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,888,032 T→A 7.0% 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 7.0% 590.5 / 24.1 170intergenic (‑41/+29)rnd/fadDribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)
Reads supporting (aligned to +/- strand):  ref base T (81/77);  new base A (7/5);  total (88/82)
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.93e-01

Only 100 of 173 total aligned reads displayed.

AGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATC  >  NC_000913/1887884‑1888179
                                                                                                                                                     |                                                                                                                                                    
 gTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTa                                                                                                                                                     >  1:994824/1‑149 (MQ=255)
     gcgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCCAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAACACGCCGGATTAA‑CCg                                                                                                                                                <  1:1158158/149‑1 (MQ=255)
      cgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCgg                                                                                                                                               <  2:793511/149‑1 (MQ=255)
         aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGt                                                                                                                                            <  2:781391/149‑1 (MQ=255)
         aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                            <  1:1320995/149‑1 (MQ=255)
         aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                            <  2:8518/149‑1 (MQ=255)
           gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCt                                                                                                                                          >  2:1106419/1‑149 (MQ=255)
           gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCt                                                                                                                                          >  3:1332082/1‑149 (MQ=255)
           gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCt                                                                                                                                          >  1:429727/1‑149 (MQ=255)
            aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGtttt                                                                                                                                         <  2:1239439/149‑4 (MQ=255)
            aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGttt                                                                                                                                          >  2:662541/1‑146 (MQ=255)
            aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGttt                                                                                                                                          <  1:662541/148‑3 (MQ=255)
            aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTg                                                                                                                                         <  2:1305018/149‑1 (MQ=255)
             aCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGGAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTga                                                                                                                                        <  3:349139/149‑1 (MQ=255)
             aCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTaa                                                                                                                                                    >  1:713059/1‑138 (MQ=255)
                aaTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGG‑TTAATCCGGCGtttttttt                                                                                                                                     <  2:1629486/149‑8 (MQ=255)
                    cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgactga                                                                                                                                 <  3:1467422/149‑1 (MQ=255)
                    cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgactga                                                                                                                                 <  3:965289/149‑1 (MQ=255)
                      gTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgactgact                                                                                                                               <  4:1232297/149‑1 (MQ=255)
                      gTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgactgact                                                                                                                               >  4:275560/1‑149 (MQ=255)
                            aGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGc                                                                                                                         >  3:1393435/1‑149 (MQ=255)
                               gCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCa                                                                                                                      <  4:701537/149‑1 (MQ=255)
                                     gCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTACGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCttt                                                                                                                <  1:502834/149‑1 (MQ=255)
                                        ggAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTAtt                                                                                                             <  4:398417/149‑1 (MQ=255)
                                           aaGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTc                                                                                                          <  3:288012/149‑1 (MQ=255)
                                                      ggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgc                                                                                               >  2:1713182/1‑149 (MQ=255)
                                                      ggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgc                                                                                               >  3:1613301/1‑149 (MQ=255)
                                                      ggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgc                                                                                               >  3:854210/1‑149 (MQ=255)
                                                         ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                            >  1:1002753/1‑149 (MQ=255)
                                                         ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                            >  4:1593668/1‑149 (MQ=255)
                                                           cacacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGAGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCt                                                                                                                  <  2:306089/125‑1 (MQ=255)
                                                             cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTc                                                                                        >  4:340480/1‑149 (MQ=255)
                                                             cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTc                                                                                        <  4:739746/149‑1 (MQ=255)
                                                              acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                       >  2:1631417/1‑149 (MQ=255)
                                                              acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                       >  2:1687341/1‑149 (MQ=255)
                                                              acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                       >  2:1209719/1‑149 (MQ=255)
                                                              acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                       >  3:119831/1‑149 (MQ=255)
                                                               caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGt                                                                                      <  4:1139655/149‑1 (MQ=255)
                                                               caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGt                                                                                      >  3:1241353/1‑149 (MQ=255)
                                                                aaaGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTAtt                                                                                                             <  3:889574/125‑1 (MQ=255)
                                                                aaaGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTAtt                                                                                                             >  4:889574/1‑125 (MQ=255)
                                                                aaaGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTc                                                                                     >  2:522336/1‑149 (MQ=255)
                                                                  agaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                   >  4:418549/1‑149 (MQ=255)
                                                                    aagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGt                                                                                 >  3:1259330/1‑149 (MQ=255)
                                                                       ccAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAt                                                                              >  3:990431/1‑149 (MQ=255)
                                                                        cAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAtt                                                                             <  1:1281003/149‑1 (MQ=255)
                                                                           cgcATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATtctc                                                                          <  1:792437/149‑1 (MQ=255)
                                                                             cATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgt                                                                        <  4:511530/149‑1 (MQ=255)
                                                                             cATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgt                                                                        <  4:1613301/149‑1 (MQ=255)
                                                                              aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                            <  2:73661/80‑1 (MQ=255)
                                                                               tCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgtcg                                                                      >  1:1377528/1‑149 (MQ=255)
                                                                                     gTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGtt                                                                                                                                           >  1:1573054/1‑73 (MQ=39)
                                                                                     gTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAt                                                                <  3:418549/149‑1 (MQ=255)
                                                                                        gTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgact                                                                                                                                   <  3:1625483/79‑1 (MQ=255)
                                                                                         tAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                            <  4:119831/149‑1 (MQ=255)
                                                                                          aaTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTc                                                           <  4:1641687/149‑1 (MQ=255)
                                                                                           aTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTcc                                                          <  4:854210/149‑1 (MQ=255)
                                                                                                 tGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGtt                                                    <  4:1620000/149‑1 (MQ=255)
                                                                                                 tGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGtt                                                    <  4:636835/149‑1 (MQ=255)
                                                                                                  ggTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTa                                                   <  4:1382927/149‑1 (MQ=255)
                                                                                                      aTTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAtt                                               <  4:1311611/149‑1 (MQ=255)
                                                                                                         cAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGATGTTAGATTTCg                                            <  4:550395/149‑1 (MQ=255)
                                                                                                           aaTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                          <  2:4230/149‑1 (MQ=255)
                                                                                                           aaTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                          >  2:1131365/1‑149 (MQ=255)
                                                                                                           aaTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                          >  3:686005/1‑149 (MQ=255)
                                                                                                            aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGc                                                                                                                                              <  2:671788/48‑1 (MQ=255)
                                                                                                            aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTaa                                                                               <  3:586177/111‑1 (MQ=255)
                                                                                                            aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTa                                         >  1:1016004/1‑149 (MQ=255)
                                                                                                             ttGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTaa                                        <  3:1593668/149‑1 (MQ=255)
                                                                                                               gTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                      >  3:572757/1‑149 (MQ=255)
                                                                                                                ttttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtc                                     >  4:172608/1‑149 (MQ=255)
                                                                                                                ttttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtc                                     >  1:416335/1‑149 (MQ=255)
                                                                                                                 tttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtca                                    <  4:1025733/149‑1 (MQ=255)
                                                                                                                  ttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcat                                   <  4:268778/149‑1 (MQ=255)
                                                                                                                  ttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcat                                   >  1:812204/1‑149 (MQ=255)
                                                                                                                  ttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcat                                   <  1:1141207/149‑1 (MQ=255)
                                                                                                                   tctcTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcatc                                  >  1:293685/1‑149 (MQ=255)
                                                                                                                     tctTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                                >  4:96072/1‑149 (MQ=255)
                                                                                                                       tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                   <  3:1292004/96‑1 (MQ=255)
                                                                                                                       tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                   >  4:1292004/1‑96 (MQ=255)
                                                                                                                       tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCg                                            <  2:1146169/135‑1 (MQ=255)
                                                                                                                       tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGa                              >  2:230095/1‑149 (MQ=255)
                                                                                                                         tAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaaa                            <  2:994824/149‑1 (MQ=255)
                                                                                                                         tAGTGAGCGTCAAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaa                             <  2:718734/149‑1 (MQ=255)
                                                                                                                                   tcgaaaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCATAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAg                   >  4:555243/4‑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‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc               >  3:930106/1‑149 (MQ=255)
                                                                                                                                       aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCACAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc               >  1:1129374/1‑149 (MQ=255)
                                                                                                                                        aaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                   >  1:1467181/1‑80 (MQ=255)
                                                                                                                                        aaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                   <  2:1467181/80‑1 (MQ=255)
                                                                                                                                           gaCGCCGGTTAAT‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa           >  2:607385/2‑149 (MQ=255)
                                                                                                                                           gaCGCCGGTTAAT‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa           >  3:1671936/2‑149 (MQ=255)
                                                                                                                                           gaCGCCGGTTAAT‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa           >  4:1068022/2‑149 (MQ=255)
                                                                                                                                           gaCGCCGG‑TTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa           >  3:821627/2‑149 (MQ=255)
                                                                                                                                               caggATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGGCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAGCTCCACCAGCTTCGGTACttt       <  4:1447263/147‑1 (MQ=255)
                                                                                                                                                cGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                            >  4:494478/1‑95 (MQ=255)
                                                                                                                                                cGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                            <  3:494478/95‑1 (MQ=255)
                                                                                                                                                    ttAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATc  <  1:230095/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATC  >  NC_000913/1887884‑1888179

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