Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,897,615 T→C 7.6% L37P (CTG→CCG) ‡ cbrB → PRK09823 family inner membrane protein, creBC regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,897,6150TC 7.6% 328.8 / 13.5 107L37P (CTG→CCG) ‡cbrBPRK09823 family inner membrane protein, creBC regulon
Reads supporting (aligned to +/- strand):  ref base T (50/47);  new base C (4/4);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

Only 100 of 107 total aligned reads displayed.

ATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCAT  >  NC_000913/3897472‑3897762
                                                                                                                                                |                                                                                                                                                   
   ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCttt                                                                                                                                             <  2:328984/149‑1 (MQ=255)
   ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCttt                                                                                                                                             <  3:82278/149‑1 (MQ=255)
   ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑AGCGGTTCttt                                                                                                                                             <  2:559318/149‑1 (MQ=255)
    gcgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGC‑GGTTCtttt                                                                                                                                            <  2:354586/149‑2 (MQ=255)
    gcgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTg                                                                                                                                            <  3:429678/149‑1 (MQ=255)
       cTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGt                                                                                                                                         <  2:300292/149‑1 (MQ=255)
       cTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGt                                                                                                                                         <  4:103201/149‑1 (MQ=255)
        tCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTg                                                                                                                                        >  2:507350/1‑149 (MQ=255)
        tCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTg                                                                                                                                        >  4:176034/1‑149 (MQ=255)
          aaaGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGat                                                                                                                                      <  4:444220/149‑1 (MQ=255)
              gtaagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGataata                                                                                                                                  >  2:547676/1‑149 (MQ=255)
               taagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATAc                                                                                                                                 >  4:543622/1‑149 (MQ=255)
               taagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATAc                                                                                                                                 <  1:18750/149‑1 (MQ=255)
                aagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACa                                                                                                                                <  3:409123/149‑1 (MQ=255)
                aagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACa                                                                                                                                >  3:312924/1‑149 (MQ=255)
                aagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACa                                                                                                                                <  4:17854/149‑1 (MQ=255)
                  gtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑AGCGGTTCtttt                                                                                                                                            <  4:156396/135‑2 (MQ=255)
                  gtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑AGCGGTTCtttt                                                                                                                                            >  3:156396/1‑134 (MQ=255)
                          aGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAAttttt                                                                                                                      <  1:216718/149‑1 (MQ=255)
                           gTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTtc                                                                                                                     <  1:507350/149‑1 (MQ=255)
                            tttGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTtct                                                                                                                    >  4:303659/1‑149 (MQ=255)
                                atatGTCTGTTTCACGTCTGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                                                <  2:260761/149‑1 (MQ=255)
                                    gTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGt                                                                                                            <  4:312924/149‑1 (MQ=255)
                                     tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                           >  1:266485/1‑149 (MQ=255)
                                     tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                           >  3:53080/1‑149 (MQ=255)
                                     tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCGTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTgcc                                                                                                           >  4:532895/1‑147 (MQ=255)
                                       tGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGa                                                                                                         >  3:98778/1‑149 (MQ=255)
                                         tttCACGTCGGGTATTACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATa                                                                                                       >  3:43438/1‑149 (MQ=255)
                                           tCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAAc                                                                                                     >  3:219305/1‑149 (MQ=255)
                                           tCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAAc                                                                                                     >  2:456570/1‑149 (MQ=255)
                                             aCGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACa                                                                                                                                <  1:110589/120‑1 (MQ=255)
                                             aCGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACa                                                                                                                                >  2:110589/1‑120 (MQ=255)
                                              cGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTAc                                                                                                  <  2:469481/149‑1 (MQ=255)
                                               gTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACg                                                                                                 >  1:519413/1‑149 (MQ=255)
                                                  gggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATATCTACGGGa                                                                                              <  4:53080/149‑1 (MQ=255)
                                                   ggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                                             >  3:495651/1‑149 (MQ=255)
                                                   ggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                                             >  1:268843/1‑149 (MQ=255)
                                                         aCATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTc                                                                                       <  4:598287/149‑1 (MQ=255)
                                                            tCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTg                                                                                    >  1:279647/1‑149 (MQ=255)
                                                              aCGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCg                                                                                  >  3:264979/1‑149 (MQ=255)
                                                                ggACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                                                <  4:219305/149‑1 (MQ=255)
                                                                ggACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                                                >  4:231117/1‑149 (MQ=255)
                                                                ggACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                                                >  1:193247/1‑149 (MQ=255)
                                                                    tttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTa                                                                            >  1:73483/1‑149 (MQ=255)
                                                                       aTTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAcc                                                                         <  1:547676/149‑1 (MQ=255)
                                                                       aTTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAcc                                                                         <  2:462589/149‑1 (MQ=255)
                                                                        ttttGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCg                                                                        >  1:255615/1‑149 (MQ=255)
                                                                          ttGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgt                                                                      <  1:567756/149‑1 (MQ=255)
                                                                              aGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTaa                                                                  <  2:262313/149‑1 (MQ=255)
                                                                                  ttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGt                                                              <  2:73483/149‑1 (MQ=255)
                                                                                   tAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACTTCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGtt                                                             <  4:43438/149‑1 (MQ=255)
                                                                                    aGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACgg                                                                                                >  1:142854/1‑113 (MQ=255)
                                                                                    aGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACgg                                                                                                <  2:142854/113‑1 (MQ=255)
                                                                                     ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGc                                                           <  2:375510/149‑1 (MQ=255)
                                                                                      gACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                                          >  1:429669/1‑149 (MQ=255)
                                                                                      gACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                                          >  1:284276/1‑149 (MQ=255)
                                                                                         cGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTtctc                                                                                                                   <  4:578972/89‑1 (MQ=255)
                                                                                         cGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTtctc                                                                                                                   >  3:578972/1‑89 (MQ=255)
                                                                                         cGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGt                                                       <  3:543622/149‑1 (MQ=255)
                                                                                          gTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTc                                                      >  2:297669/1‑149 (MQ=255)
                                                                                           ttAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCt                                                     <  4:495651/149‑1 (MQ=255)
                                                                                           ttAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCt                                                     >  1:291888/1‑149 (MQ=255)
                                                                                               ttGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcc                                                  <  3:343689/148‑1 (MQ=255)
                                                                                               ttGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcc                                                  >  4:343689/1‑148 (MQ=255)
                                                                                                  gtgtTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaa                                              <  4:264979/149‑1 (MQ=255)
                                                                                                    gtTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaaa                                            <  2:350760/149‑1 (MQ=255)
                                                                                                     ttCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAg                                           >  2:334406/1‑149 (MQ=255)
                                                                                                       cTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGAt                                         >  1:407114/1‑149 (MQ=255)
                                                                                                       cTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTATAATGGTTGCCTGTCTGCCGGAAAAGAt                                         <  3:532895/149‑1 (MQ=255)
                                                                                                        tGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATc                                        >  3:371450/1‑149 (MQ=255)
                                                                                                          tttCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                      <  4:252472/149‑1 (MQ=255)
                                                                                                          tttCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                      >  1:607669/1‑149 (MQ=255)
                                                                                                           ttCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCggg                                     <  3:303659/149‑1 (MQ=255)
                                                                                                               tGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTcac                                 <  3:231117/149‑1 (MQ=255)
                                                                                                                  ccTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGa                              >  2:317743/1‑149 (MQ=255)
                                                                                                                   ctctACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTAc                                                                                                  <  4:106441/80‑1 (MQ=255)
                                                                                                                   ctctACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTAc                                                                                                  >  3:106441/1‑80 (MQ=255)
                                                                                                                    tctACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATCATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgt                                                                                 <  3:414416/96‑1 (MQ=255)
                                                                                                                    tctACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgt                                                                                 >  4:414416/1‑96 (MQ=255)
                                                                                                                    tctACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaaa                            >  2:258129/1‑149 (MQ=255)
                                                                                                                     ctACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAAc                           <  2:279647/149‑1 (MQ=255)
                                                                                                                      tACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACGTCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATAGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACg                          >  2:375890/1‑149 (MQ=255)
                                                                                                                        cATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaaa                            <  1:81436/145‑1 (MQ=255)
                                                                                                                        cATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaaa                            >  2:81436/1‑145 (MQ=255)
                                                                                                                        cATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGtt                        >  4:584583/1‑149 (MQ=255)
                                                                                                                         atatTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa                       >  3:291451/1‑149 (MQ=255)
                                                                                                                          tatTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTAt                      >  4:486050/1‑149 (MQ=255)
                                                                                                                           atTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATc                     >  4:525824/1‑149 (MQ=255)
                                                                                                                           atTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATc                     <  1:297669/149‑1 (MQ=255)
                                                                                                                             tcttcGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGt                   >  1:576181/1‑149 (MQ=255)
                                                                                                                                 cGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTc               >  1:19684/1‑149 (MQ=255)
                                                                                                                                  gCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTct              >  1:343650/1‑149 (MQ=255)
                                                                                                                                  gCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTct              <  1:334406/149‑1 (MQ=255)
                                                                                                                                   ccaaaGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg             >  3:377569/3‑149 (MQ=255)
                                                                                                                                   ccaaaGAACC‑AGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgtt                                                                             <  1:166657/83‑1 (MQ=255)
                                                                                                                                   ccaaaGAACC‑AGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                                                <  1:6072/80‑1 (MQ=255)
                                                                                                                                   ccaaaGAACC‑AGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                                                >  2:6072/3‑82 (MQ=255)
                                                                                                                                    aaaaGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctgg            >  2:178133/1‑149 (MQ=255)
                                                                                                                                          aCC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCtggtg      <  3:176034/149‑1 (MQ=255)
                                                                                                                                               cTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCAt  >  3:68393/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCAT  >  NC_000913/3897472‑3897762

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: ‑