Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,858 G→C 6.8% intergenic (‑100/+15) yecJ ← / ← azuC DUF2766 family protein/acid‑inducible small membrane‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,8580GC 6.8% 423.0 / 12.4 125intergenic (‑100/+15)yecJ/azuCDUF2766 family protein/acid‑inducible small membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (65/51);  new base C (3/6);  total (68/57)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 125 total aligned reads displayed.

CCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGC  >  NC_000913/1987713‑1988005
                                                                                                                                                  |                                                                                                                                                    
 cacaTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGc                                                                                                                                                  >  1:314478/1‑149 (MQ=255)
  acaTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCa                                                                                                                                                 <  4:414796/149‑1 (MQ=255)
   caTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCATCATCAGGACACAAAAAAACCTGC‑CGGAGCAg                                                                                                                                                <  2:634236/149‑1 (MQ=255)
    aTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCAGCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAgg                                                                                                                                               >  4:684812/1‑149 (MQ=255)
    aTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAgg                                                                                                                                               >  4:338827/1‑149 (MQ=255)
     tCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑t                                                                                                                                             >  4:262108/1‑149 (MQ=255)
      cAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tt                                                                                                                                            >  2:266117/1‑149 (MQ=255)
        gAAACCAATTCCTGATCGGCATTCAGCGTCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCTACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttt                                                                                                                                          >  2:149857/1‑149 (MQ=255)
        gAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttt                                                                                                                                          >  2:941049/1‑149 (MQ=255)
          aaCCAATTCCTGGTCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATTGCATTTTGAAAGCCGGTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttttt                                                                                                                                        <  1:97316/149‑1 (MQ=255)
          aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttttt                                                                                                                                        <  3:444321/149‑1 (MQ=255)
          aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttttt                                                                                                                                        <  4:19520/149‑1 (MQ=255)
          aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttttt                                                                                                                                        <  2:846415/149‑1 (MQ=255)
          aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttttt                                                                                                                                        <  3:62241/149‑1 (MQ=255)
          aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttttt                                                                                                                                        <  3:386689/149‑1 (MQ=255)
          aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGAATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGG‑tttttt                                                                                                                                        <  4:196511/149‑1 (MQ=255)
           aCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGG‑ttttttt                                                                                                                                       <  1:253174/149‑2 (MQ=255)
            ccAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGt                                                                                                                                      >  4:531129/1‑149 (MQ=255)
            ccAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACAC‑AAAAAACCTGCTCCG‑GCAGG‑tttttttgt                                                                                                                                     <  2:106415/149‑4 (MQ=255)
             cAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTCATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGG‑tttttttgt                                                                                                                                     <  3:896988/149‑4 (MQ=255)
              aaTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAGCTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTa                                                                                                                                    <  3:684812/149‑1 (MQ=255)
              aaTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTa                                                                                                                                    >  1:493188/1‑149 (MQ=255)
              aaTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTa                                                                                                                                    >  4:145378/1‑149 (MQ=255)
               aTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTAt                                                                                                                                   >  3:158938/1‑149 (MQ=255)
                     gATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAAc                                                                                                                             >  1:85913/1‑149 (MQ=255)
                       tCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACat                                                                                                                           <  2:85913/149‑1 (MQ=255)
                          gCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATAtt                                                                                                                        >  1:218981/1‑149 (MQ=255)
                                   ggCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCTTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCgg                                                                                                               >  1:825529/1‑149 (MQ=255)
                                   ggCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCgg                                                                                                               <  4:878564/149‑1 (MQ=255)
                                           cATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCtt                                                                                                       >  1:35422/1‑149 (MQ=255)
                                             tCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTg                                                                                                     >  2:1002466/1‑149 (MQ=255)
                                               gTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGaa                                                                                                   >  3:749907/1‑149 (MQ=255)
                                                     tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCtt                                                                                             >  2:1031130/1‑149 (MQ=255)
                                                     tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCtt                                                                                             <  4:750994/149‑1 (MQ=255)
                                                          aTAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAAt                                                                                        >  2:816721/1‑149 (MQ=255)
                                                          aTAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAAt                                                                                        >  3:399402/1‑149 (MQ=255)
                                                           tAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                                       <  2:120738/149‑1 (MQ=255)
                                                           tAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                                       >  1:120738/1‑149 (MQ=255)
                                                             gCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGt                                                                                     >  4:631865/1‑149 (MQ=255)
                                                              cATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGaa                                                                                                   <  3:1017106/134‑1 (MQ=255)
                                                              cATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGaa                                                                                                   >  4:1017106/1‑134 (MQ=255)
                                                              cATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGtt                                                                                    <  1:941049/149‑1 (MQ=255)
                                                               aTTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTa                                                                                   <  4:749907/149‑1 (MQ=255)
                                                               aTTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTa                                                                                   >  2:6695/1‑149 (MQ=255)
                                                                ttttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTAt                                                                                  <  4:158938/149‑1 (MQ=255)
                                                                ttttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTAt                                                                                  >  2:141399/1‑149 (MQ=255)
                                                                ttttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTAt                                                                                  >  2:184648/1‑149 (MQ=255)
                                                                 tttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTAtt                                                                                 <  3:682096/149‑1 (MQ=255)
                                                                   tGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGa                                                                               >  4:190785/1‑149 (MQ=255)
                                                                    gAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGaa                                                                              <  3:531129/149‑1 (MQ=255)
                                                                    gAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGaa                                                                              <  2:352444/149‑1 (MQ=255)
                                                                      aaGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACa                                                                            >  2:317734/1‑149 (MQ=255)
                                                                        gCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCtt                                                                                                       >  2:926560/1‑120 (MQ=255)
                                                                        gCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCtt                                                                                                       <  1:926560/120‑1 (MQ=255)
                                                                           gTTATAGTAGCGACTTCACATCTTCAGCGATAGTGACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACtt                                                                       <  1:873511/149‑1 (MQ=255)
                                                                           gTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCG‑GCAGG‑tttttt                                                                                                                                        <  4:427315/84‑1 (MQ=37)
                                                                           gTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCG‑GCAGG‑tttttt                                                                                                                                        >  3:427315/1‑84 (MQ=37)
                                                                            ttATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACttt                                                                      >  1:786614/1‑149 (MQ=255)
                                                                            ttATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACttt                                                                      <  2:527029/149‑1 (MQ=255)
                                                                               tagtagCGACTTCACATCTTCCGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCa                                                                   <  2:810614/149‑1 (MQ=255)
                                                                               tagtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCa                                                                   <  1:184648/149‑1 (MQ=255)
                                                                               tagtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCa                                                                   >  2:189646/1‑149 (MQ=255)
                                                                                 gtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAgg                                                                 <  1:927671/149‑1 (MQ=255)
                                                                                  tagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGa                                                                >  2:862403/1‑149 (MQ=255)
                                                                                    gCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtt                                                              >  1:11349/1‑149 (MQ=255)
                                                                                    gCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtt                                                              <  2:62024/149‑1 (MQ=255)
                                                                                     cGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAttt                                                             <  4:127583/149‑1 (MQ=255)
                                                                                       aCTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTg                                                           >  1:247547/1‑149 (MQ=255)
                                                                                         ttCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGA‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTgcg                                                         >  4:719362/1‑149 (MQ=255)
                                                                                             caTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGt                                                     <  2:517094/149‑1 (MQ=255)
                                                                                              aTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTa                                                                                   >  2:909169/1‑118 (MQ=255)
                                                                                               tCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGttt                                                   >  4:157077/1‑149 (MQ=255)
                                                                                                     gCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATcgc                                             >  1:147976/1‑149 (MQ=255)
                                                                                                     gCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATcgc                                             <  4:399402/149‑1 (MQ=255)
                                                                                                           gTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACt                                       >  2:160166/1‑149 (MQ=255)
                                                                                                            tCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACata                                                                                                                          <  3:534223/65‑1 (MQ=255)
                                                                                                            tCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACata                                                                                                                          >  4:534223/1‑65 (MQ=255)
                                                                                                              acaTCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCg                                    >  2:5471/1‑149 (MQ=255)
                                                                                                              acaTCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCg                                    >  3:165938/1‑149 (MQ=255)
                                                                                                                aTCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                                  <  3:631865/149‑1 (MQ=255)
                                                                                                                 tCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAcc                                 >  3:166940/1‑149 (MQ=255)
                                                                                                                 tCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAcc                                 >  4:693101/1‑149 (MQ=255)
                                                                                                                    aCCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAtt                              <  3:1032712/149‑1 (MQ=255)
                                                                                                                      cGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTg                            <  1:317734/149‑1 (MQ=255)
                                                                                                                            caGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATAc                      >  4:105028/1‑149 (MQ=255)
                                                                                                                              ggACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAg                    <  1:5471/149‑1 (MQ=255)
                                                                                                                                acacAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGTGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgt                  <  2:314478/149‑1 (MQ=255)
                                                                                                                                acacAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgt                  <  4:854709/149‑1 (MQ=255)
                                                                                                                                acacAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgt                  <  4:1028051/149‑1 (MQ=255)
                                                                                                                                acacAAAAAAACCTGC‑CGGAGCAGG‑TTTTGTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgt                  <  2:218981/149‑1 (MQ=255)
                                                                                                                                 cacaAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGTGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtg                 >  1:1060783/1‑149 (MQ=255)
                                                                                                                                  acaAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtgt                >  3:699043/1‑149 (MQ=255)
                                                                                                                                  acaAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtgt                >  4:630013/1‑149 (MQ=255)
                                                                                                                                  acaAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtgt                >  4:709327/1‑149 (MQ=255)
                                                                                                                                    aaaaaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTCGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa              >  1:308235/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCCTGC‑TCCGGCAGGTTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa              >  1:489502/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAt             >  2:214619/1‑149 (MQ=255)
                                                                                                                                       aaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTg           >  2:137222/1‑149 (MQ=255)
                                                                                                                                        aaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGt          >  1:583375/1‑149 (MQ=255)
                                                                                                                                             tGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTg     >  2:960193/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGC  >  NC_000913/1987713‑1988005

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