Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,897,615 T→C 6.5% 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 6.5% 575.0 / 18.2 182L37P (CTG→CCG) ‡cbrBPRK09823 family inner membrane protein, creBC regulon
Reads supporting (aligned to +/- strand):  ref base T (78/92);  new base C (7/5);  total (85/97)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 183 total aligned reads displayed.

TTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATA  >  NC_000913/3897468‑3897763
                                                                                                                                                   |                                                                                                                                                    
ttAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTg                                                                                                                                                     <  2:792476/149‑1 (MQ=255)
ttAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTg                                                                                                                                                     >  1:586320/1‑149 (MQ=255)
 tAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTgg                                                                                                                                                    <  1:661018/149‑1 (MQ=255)
   aaTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGtt                                                                                                                                                  >  3:151271/1‑149 (MQ=255)
      tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCtt                                                                                                                                               <  3:687751/149‑1 (MQ=255)
        gcgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTg                                                                                                                                             <  1:38642/149‑1 (MQ=255)
        gcgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGCGGTTCtttt                                                                                                                                             <  3:345087/149‑2 (MQ=255)
        gcgcTAAAGGGTAAGTAGGGGGAGTTAGATATGTCTGTTGCACGTCGGGTAACACATCACGGACTATATTTTCCAGATTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTg                                                                                                                                             >  2:181912/1‑149 (MQ=255)
         cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTgg                                                                                                                                            <  3:850597/149‑1 (MQ=255)
         cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTgg                                                                                                                                            <  4:20080/149‑1 (MQ=255)
          gcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTggg                                                                                                                                           <  2:837751/149‑1 (MQ=255)
            tCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTgg                                                                                                                                            <  3:508854/146‑1 (MQ=255)
            tCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTg                                                                                                                                         >  2:537091/1‑149 (MQ=255)
            tCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTg                                                                                                                                         >  3:115754/1‑149 (MQ=255)
            tCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTg                                                                                                                                         >  4:129964/1‑149 (MQ=255)
            tCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTg                                                                                                                                         >  2:50379/1‑149 (MQ=255)
                aGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGataa                                                                                                                                     <  1:79047/149‑1 (MQ=255)
                 ggTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGataat                                                                                                                                    <  4:561824/149‑1 (MQ=255)
                       taaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTACAGTTTTAGGACAGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATcc                                                                                                                              >  4:1019820/1‑149 (MQ=255)
                         gaGGGAGTTTGATGTGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTACTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACCTCCaa                                                                                                                            >  4:682054/2‑149 (MQ=255)
                          aGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAt                                                                                                                           <  1:149896/149‑1 (MQ=255)
                                tttGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTtct                                                                                                                     >  1:161464/1‑149 (MQ=255)
                                tttGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTtct                                                                                                                     <  2:196321/149‑1 (MQ=255)
                                    atatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                                                 <  1:549359/149‑1 (MQ=255)
                                         tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGTCCGTTACTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                            >  1:982198/1‑149 (MQ=255)
                                         tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                            >  1:1012654/1‑149 (MQ=255)
                                         tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                            >  4:927114/1‑149 (MQ=255)
                                            gTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGAt                                                                                                         <  3:546797/149‑1 (MQ=255)
                                               tCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                            >  1:543839/1‑143 (MQ=255)
                                               tCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                            <  2:543839/143‑1 (MQ=255)
                                                   gTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACg                                                                                                  <  3:637228/149‑1 (MQ=255)
                                                   gTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACg                                                                                                  >  3:79292/1‑149 (MQ=255)
                                                    tCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACgg                                                                                                 <  1:199187/149‑1 (MQ=255)
                                                      gggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGa                                                                                               >  1:312241/1‑149 (MQ=255)
                                                      gggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGa                                                                                               <  3:940570/149‑1 (MQ=255)
                                                       ggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                                              >  1:642132/1‑149 (MQ=255)
                                                         tAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTGTTGTCCTCTACCTATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTAAGGGAGCt                                                                                            <  2:11500/149‑1 (MQ=255)
                                                             acgtCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATCTTTCTCTTATTGTCGATAACTACGGGAGCTACTc                                                                                        <  3:927114/146‑1 (MQ=255)
                                                                 cACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGc                                                                                    <  4:151271/149‑1 (MQ=255)
                                                                  aCGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGataa                                                                                                                                     <  2:694872/99‑1 (MQ=255)
                                                                  aCGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAttttt                                                                                                                       <  1:920355/113‑1 (MQ=255)
                                                                    ggACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                                                 <  2:1012406/149‑1 (MQ=255)
                                                                     gACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttg                                                                                >  1:947900/1‑149 (MQ=255)
                                                                      aCTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                                                 >  1:218485/1‑147 (MQ=255)
                                                                      aCTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                                                 <  2:218485/147‑1 (MQ=255)
                                                                       cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgtt                                                                              <  1:725409/149‑1 (MQ=255)
                                                                        tttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTa                                                                             >  2:945253/1‑149 (MQ=255)
                                                                         ttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTaa                                                                            <  2:6729/149‑1 (MQ=255)
                                                                          tATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAc                                                                           <  3:488243/149‑1 (MQ=255)
                                                                               tGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTCATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgtg                                                                      <  2:143165/149‑1 (MQ=255)
                                                                                   gTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAAt                                                                  >  4:398864/1‑149 (MQ=255)
                                                                                    ttttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgt                                                                       >  3:574082/1‑143 (MQ=255)
                                                                                    ttttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATg                                                                 >  4:19046/1‑149 (MQ=255)
                                                                                    ttttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATg                                                                 >  1:543321/1‑149 (MQ=255)
                                                                                    ttttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATg                                                                 <  2:586320/149‑1 (MQ=255)
                                                                                         gtaCCGTTAATTGGTCTTCTGGCTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTTGTTCTTTGGGTGTTAATCTATACACTTTTTCTCTTAGTGTCCATAACTACGGGAGCGATACCCGCGTTGTTGCCCGGTGTAATGGTTgg                                                            <  2:982198/147‑2 (MQ=255)
                                                                                         ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGc                                                            <  3:19046/149‑1 (MQ=255)
                                                                                           aCCGTTAATTGGTGTTCTGTTTCTCGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCt                                                          <  2:766937/149‑1 (MQ=255)
                                                                                            ccGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTg                                                         >  1:366530/1‑149 (MQ=255)
                                                                                              gTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTc                                                       >  2:828285/1‑149 (MQ=255)
                                                                                                  aTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcc                                                   >  4:638219/1‑149 (MQ=255)
                                                                                                  aTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTCTAATGGTTGCCTGTCTGcc                                                   <  4:917773/149‑1 (MQ=255)
                                                                                                   ttGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCg                                                  >  2:105867/1‑149 (MQ=255)
                                                                                                         ttCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAg                                            >  1:720293/1‑149 (MQ=255)
                                                                                                         ttCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAg                                            >  2:991946/1‑149 (MQ=255)
                                                                                                            tGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATc                                         <  3:943229/149‑1 (MQ=255)
                                                                                                             gTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCg                                        <  4:79292/149‑1 (MQ=255)
                                                                                                              tttCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                       >  3:303874/1‑149 (MQ=255)
                                                                                                              tttCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                       >  3:373276/1‑149 (MQ=255)
                                                                                                              tttCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                       <  4:327140/149‑1 (MQ=255)
                                                                                                              tttCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                       <  2:161464/149‑1 (MQ=255)
                                                                                                              tttCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTGCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                       <  4:832712/149‑1 (MQ=255)
                                                                                                                   tGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTcac                                  <  4:507316/149‑1 (MQ=255)
                                                                                                                     tCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgt                                                                                  <  4:624917/99‑1 (MQ=255)
                                                                                                                     tCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgt                                                                                  >  3:624917/1‑99 (MQ=255)
                                                                                                                     tCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAAt                                                                  >  3:97083/1‑115 (MQ=255)
                                                                                                                     tCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAg                                >  4:51042/1‑149 (MQ=255)
                                                                                                                       ctctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                       <  2:287726/140‑1 (MQ=255)
                                                                                                                       ctctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaa                              >  4:960123/1‑149 (MQ=255)
                                                                                                                         ctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTc                                                                                        >  2:317736/1‑89 (MQ=255)
                                                                                                                         ctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAAc                            >  1:619934/1‑149 (MQ=255)
                                                                                                                          tACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACg                           >  3:87507/1‑149 (MQ=255)
                                                                                                                             atatTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGc                                                    <  2:721059/121‑1 (MQ=255)
                                                                                                                             atatTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa                        <  1:1061059/149‑1 (MQ=255)
                                                                                                                               atTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATc                      <  2:642132/149‑1 (MQ=255)
                                                                                                                               atTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATc                      <  4:687455/149‑1 (MQ=255)
                                                                                                                                ttcttcGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCg                     >  1:110559/1‑149 (MQ=255)
                                                                                                                                  cttcGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGtt                   <  2:619934/149‑1 (MQ=255)
                                                                                                                                  atatcccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGtt                   <  3:287087/142‑1 (MQ=255)
                                                                                                                                    tcGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGt                 <  4:263273/149‑1 (MQ=255)
                                                                                                                                      gCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAATGATCGGGTCACAGAAACGTTATCGTTGTct               >  1:691207/1‑149 (MQ=255)
                                                                                                                                       ccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg              >  2:402270/3‑149 (MQ=255)
                                                                                                                                       ccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg              >  2:573449/3‑149 (MQ=255)
                                                                                                                                        caaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTGCGGGAGCTATTCCTGCGTTGTTAACCGGTGTAAt                                                                  >  1:815107/2‑96 (MQ=255)
                                                                                                                                        caaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAAt                                                                  <  2:815107/95‑1 (MQ=255)
                                                                                                                                        caaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctgg             >  2:530583/2‑149 (MQ=255)
                                                                                                                                         aaaGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctggc            <  2:366530/149‑1 (MQ=255)
                                                                                                                                         aaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctggc            >  1:10502/1‑149 (MQ=255)
                                                                                                                                             aaCCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGCGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCtggt        <  3:1040636/149‑1 (MQ=255)
                                                                                                                                                   tGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATa  >  4:741537/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATA  >  NC_000913/3897468‑3897763

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