Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,897,613 G→A 5.9% P36P (CCG→CCA)  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,6130GA 5.9% 547.7 / 14.7 168P36P (CCG→CCA) cbrBPRK09823 family inner membrane protein, creBC regulon
Reads supporting (aligned to +/- strand):  ref base G (76/82);  new base A (8/2);  total (84/84)
Fisher's exact test for biased strand distribution p-value = 9.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

CGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCA  >  NC_000913/3897465‑3897761
                                                                                                                                                     |                                                                                                                                                    
 gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑Gc                                                                                                                                                     >  3:381633/1‑149 (MQ=255)
   ttAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTg                                                                                                                                                   >  4:779979/1‑149 (MQ=255)
    tAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGGCTGTTTCACGTCGGGTAATACATCGCGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTgg                                                                                                                                                  >  4:623905/1‑149 (MQ=255)
          ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGC‑GGTTCttt                                                                                                                                            <  4:77178/149‑1 (MQ=255)
          ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCttt                                                                                                                                            <  3:779979/149‑1 (MQ=255)
          ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCttt                                                                                                                                            <  4:198030/149‑1 (MQ=255)
          ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑AGCGGTTCttt                                                                                                                                            <  4:531768/149‑1 (MQ=255)
            cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTgg                                                                                                                                          <  1:165545/149‑1 (MQ=255)
            cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGCTCTTTgg                                                                                                                                          <  2:599764/149‑1 (MQ=255)
             gcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTggg                                                                                                                                         >  3:380924/1‑149 (MQ=255)
               tCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTg                                                                                                                                       >  1:912195/1‑149 (MQ=255)
                cAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGa                                                                                                                                      <  3:553317/149‑1 (MQ=255)
                  aaGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGCGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGata                                                                                                                                    <  3:623905/149‑1 (MQ=255)
                   atgtaagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGataa                                                                                                                                   <  4:928518/147‑1 (MQ=255)
                     gtaagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGataata                                                                                                                                 >  3:355622/1‑149 (MQ=255)
                        agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACAt                                                                                                                              <  1:404353/149‑1 (MQ=255)
                        agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACAt                                                                                                                              <  1:660030/149‑1 (MQ=255)
                            aaGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACa                                                                                                                               <  2:169156/144‑1 (MQ=255)
                            aaGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCaa                                                                                                                          <  4:795686/149‑1 (MQ=255)
                            aaGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCaa                                                                                                                          >  2:129638/1‑149 (MQ=255)
                                gAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAAtttt                                                                                                                      <  4:105842/149‑1 (MQ=255)
                                   tttGATATATCTGTGTGACCTTGGGTAAAACATCACGGACTTTTTTTTGCACATTTAGCACCGTTAATTGGTGTTCTGTTTCTTGTTTTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATACAATTTTtct                                                                                                                   >  2:863200/1‑149 (MQ=255)
                                       atatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                                               <  3:834195/149‑1 (MQ=255)
                                         atGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTgg                                                                                                                                          >  1:758940/1‑120 (MQ=255)
                                         atGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTgg                                                                                                                                          <  2:758940/120‑1 (MQ=255)
                                         atGTCTGTTTAACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAtt                                                                                                             <  3:448527/149‑1 (MQ=255)
                                           gTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGt                                                                                                           <  4:401990/149‑1 (MQ=255)
                                            tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGt                                                                                                           <  2:68907/148‑1 (MQ=255)
                                            tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                          <  3:131718/149‑1 (MQ=255)
                                            tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                          <  3:138143/149‑1 (MQ=255)
                                            tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑AGCGGTTCtttt                                                                                                                                           <  1:486492/116‑2 (MQ=255)
                                            tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATCGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑AGCGGTTCtttt                                                                                                                                           >  2:486492/1‑115 (MQ=255)
                                              tGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGa                                                                                                        <  1:868110/149‑1 (MQ=255)
                                                    aCGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATaa                                                                                                     <  3:533566/146‑1 (MQ=255)
                                                       tCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGGCCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACgg                                                                                               >  2:849599/1‑149 (MQ=255)
                                                       tCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACgg                                                                                               <  3:1076516/149‑1 (MQ=255)
                                                         gggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGa                                                                                             <  3:1009834/149‑1 (MQ=255)
                                                           gTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATACTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGc                                                                                           <  2:485159/149‑1 (MQ=255)
                                                            tAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCt                                                                                          >  1:686288/1‑149 (MQ=255)
                                                                aCATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTc                                                                                      <  4:353228/149‑1 (MQ=255)
                                                                   tCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTg                                                                                   >  2:414546/1‑149 (MQ=255)
                                                                    cACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGc                                                                                  <  1:1035013/149‑1 (MQ=255)
                                                                    cACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGc                                                                                  <  2:1036542/149‑1 (MQ=255)
                                                                     aCGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCg                                                                                 >  3:970333/1‑149 (MQ=255)
                                                                       ggACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGTGCTATTCCTGCgat                                                                               >  3:647542/1‑147 (MQ=255)
                                                                       ggACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCAATAACTACGGGAGCTATTCCTGCgtt                                                                               <  2:379087/149‑1 (MQ=255)
                                                                        gACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttg                                                                              >  1:184800/1‑149 (MQ=255)
                                                                         aCTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgt                                                                             <  1:309514/149‑1 (MQ=255)
                                                                           tttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTa                                                                           >  2:38033/1‑149 (MQ=255)
                                                                           tttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTa                                                                           >  4:761287/1‑149 (MQ=255)
                                                                           attATTGTGTAGTTTTAGTATCGTTAATTGGTGTTCTGTTCCTTGTTCTCTACATATTCTTCGCAAAAGAACT‑GTTGGTTGTTTGGATGATAATACATCTAATTTTTCTGATATTGTTGGTAACTACGGGAGCTATGCCTGAgtgtttc                                                                           >  4:263934/2‑144 (MQ=255)
                                                                             tATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAc                                                                         <  4:528919/149‑1 (MQ=255)
                                                                                      gTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAAt                                                                <  4:889071/149‑1 (MQ=255)
                                                                                           aGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTg                                                           <  2:832386/149‑1 (MQ=255)
                                                                                            ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGc                                                          >  4:82430/1‑149 (MQ=255)
                                                                                            ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGc                                                          <  3:720249/149‑1 (MQ=255)
                                                                                             gACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                                         >  2:633968/1‑149 (MQ=255)
                                                                                             gACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                                         >  4:178162/1‑149 (MQ=255)
                                                                                             gACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                                         <  3:231235/149‑1 (MQ=255)
                                                                                             gACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                                         >  4:399728/1‑149 (MQ=255)
                                                                                                cGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATg                                                               >  4:567473/1‑140 (MQ=255)
                                                                                                cGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGt                                                      >  2:868538/1‑149 (MQ=255)
                                                                                                  ttAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCt                                                    <  1:745171/149‑1 (MQ=255)
                                                                                                          tgtTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTCTTCTCTTTTTGTCGATAATAAAGGGCGCTATTCCTACGTTTTTAACCGGTGTAATGGTTGCATGTCTGCCGGaaa                                            >  1:837094/1‑149 (MQ=255)
                                                                                                               tGTTTCTTGTCCTCTACATATTCTTCTCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTCGCCTGTCTGCCGGAAAAGATc                                       >  4:688135/1‑149 (MQ=255)
                                                                                                               tGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTaa                                                                 <  3:1078131/123‑1 (MQ=255)
                                                                                                               tGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTaa                                                                 >  4:1078131/1‑123 (MQ=255)
                                                                                                                gTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCAGTCGGAAAAGATCg                                      <  3:82430/149‑1 (MQ=255)
                                                                                                                 tttCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                                               <  2:351766/75‑1 (MQ=255)
                                                                                                                  ttCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCggg                                    <  1:849599/149‑1 (MQ=255)
                                                                                                                   tCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGt                                   <  4:980607/149‑1 (MQ=255)
                                                                                                                   tCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGt                                   <  3:76224/149‑1 (MQ=255)
                                                                                                                   tCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGt                                   <  2:27341/149‑1 (MQ=255)
                                                                                                                   tCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGt                                   <  3:178162/149‑1 (MQ=255)
                                                                                                                        tCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGAt                                        <  4:375693/139‑1 (MQ=255)
                                                                                                                         ccTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAAGGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGa                             >  4:705334/1‑149 (MQ=255)
                                                                                                                          ctctACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCTTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaa                            >  3:752551/1‑149 (MQ=255)
                                                                                                                           tctACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaaa                           <  2:166748/149‑1 (MQ=255)
                                                                                                                                atatTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTTCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa                      >  1:263029/1‑149 (MQ=255)
                                                                                                                                atatTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa                      <  1:38033/149‑1 (MQ=255)
                                                                                                                                atatTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa                      >  1:63963/1‑149 (MQ=255)
                                                                                                                                 tatTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTAt                     >  1:314346/1‑149 (MQ=255)
                                                                                                                                   ttcttcGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCg                   >  1:409594/1‑149 (MQ=255)
                                                                                                                                   ttcttcGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCg                   >  2:233036/1‑149 (MQ=255)
                                                                                                                                         gCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTct             <  2:263029/149‑1 (MQ=255)
                                                                                                                                         gCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTct             <  1:1037490/149‑1 (MQ=255)
                                                                                                                                         gCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTct             <  2:184800/149‑1 (MQ=255)
                                                                                                                                          ccaaaGAACC‑AGCGGTTCTTTGGGTGATAATATATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg            >  2:942082/3‑149 (MQ=255)
                                                                                                                                          ccaaaGAACC‑AGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGttagctg            >  4:1012485/3‑144 (MQ=255)
                                                                                                                                          ccaaaGAACC‑AGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg            >  1:942057/3‑149 (MQ=255)
                                                                                                                                          ccaaaGAACC‑AGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg            >  3:21836/3‑149 (MQ=255)
                                                                                                                                          cAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                                               <  4:466803/50‑1 (MQ=255)
                                                                                                                                          cAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                                               >  3:466803/1‑50 (MQ=255)
                                                                                                                                          cAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATc                                       >  1:182199/1‑122 (MQ=255)
                                                                                                                                           caaaGAACC‑AGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctgg           >  4:239602/2‑149 (MQ=255)
                                                                                                                                           aaaaGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctgg           >  4:224482/1‑149 (MQ=255)
                                                                                                                                            aaaGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAGACGTTATCGTTGTctggc          <  3:670539/149‑1 (MQ=255)
                                                                                                                                                aaCC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCtggt      >  1:826373/1‑149 (MQ=255)
                                                                                                                                                   c‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGc   >  4:352871/1‑149 (MQ=255)
                                                                                                                                                     gCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGGTGGTGGCa  >  1:343367/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACC‑GCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCA  >  NC_000913/3897465‑3897761

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