Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,581,360 A→G 5.0% Y389H (TAT→CAT)  ydeN ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,581,3600AG 5.0% 248.1 / 2.7 79Y389H (TAT→CAT) ydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (33/42);  new base G (2/2);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

TCAGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAGGTAAAGAGAA  >  NC_000913/1581213‑1581506
                                                                                                                                                   |                                                                                                                                                  
tCAGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATaa                                                                                                                                                   <  7:47014/149‑1 (MQ=255)
    aTTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTAc                                                                                                                                               <  2:241611/149‑1 (MQ=255)
     tttttATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTAcc                                                                                                                                              <  1:220912/149‑1 (MQ=255)
        ttATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCgg                                                                                                                                            <  5:49638/148‑1 (MQ=255)
           tgtgGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGtt                                                                                                                                         >  5:212166/1‑148 (MQ=255)
           tgtgGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGtt                                                                                                                                         <  6:198011/148‑1 (MQ=255)
             tgGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGaa                                                                                                                                      <  7:178965/149‑1 (MQ=255)
             tgGCTCGCCGTGTTTGTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGGCATAATTACCGGGTTGaa                                                                                                                                      <  2:59882/149‑1 (MQ=255)
             tgGCTAGACTTGTTTCTTATATTGCAACCATGGCAGCAAGGAAACGCCATACAGCTTAATGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTTTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGaa                                                                                                                                      <  1:78226/149‑1 (MQ=255)
              ggCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAg                                                                                                                                     >  7:350033/1‑149 (MQ=255)
               gCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAg                                                                                                                                     <  8:322593/148‑1 (MQ=255)
               gCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGt                                                                                                                                    >  1:378370/1‑149 (MQ=255)
               gCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGt                                                                                                                                    >  7:211045/1‑149 (MQ=255)
                 tCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGttt                                                                                                                                  >  2:355431/1‑149 (MQ=255)
                   gCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTtttc                                                                                                                                <  6:348918/149‑1 (MQ=255)
                     cTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTtttcct                                                                                                                              >  5:197341/1‑149 (MQ=255)
                           tCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTccac                                                                                                                        >  8:307688/1‑149 (MQ=255)
                           tCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTccac                                                                                                                        >  5:198807/1‑149 (MQ=255)
                           tCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTcca                                                                                                                         <  2:378370/148‑1 (MQ=255)
                           tCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTcca                                                                                                                         <  3:39617/148‑1 (MQ=255)
                            cTTATCTTGCAACCAGGGCAGCAAGGAAACCCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCggg                                                                                                                                           <  4:16984/129‑1 (MQ=255)
                            cTTATCTTGCAACCAGGGCAGCAAGGAAACCCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCggg                                                                                                                                           >  3:16984/1‑129 (MQ=255)
                               aTCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCcaccac                                                                                                                     <  6:198807/148‑1 (MQ=255)
                                      aaCCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACAt                                                                                                             >  4:90038/1‑149 (MQ=255)
                                       aCCAGGGCTGCAAGGAAACGCCATCCAGCTTAAGGTATTTTGGATTGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTGATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACAtt                                                                                                            <  1:210846/149‑1 (MQ=255)
                                         cAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTgg                                                                                                          <  1:355431/149‑1 (MQ=255)
                                         cAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTgg                                                                                                          <  5:138233/149‑1 (MQ=255)
                                          aGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTggg                                                                                                         >  5:363600/1‑149 (MQ=255)
                                          aGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTggg                                                                                                         <  3:90038/149‑1 (MQ=255)
                                          aGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATAAATAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTggg                                                                                                         <  2:7695/149‑1 (MQ=255)
                                          aGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCGGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTggg                                                                                                         <  7:72869/149‑1 (MQ=255)
                                              cagcaAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGgtgt                                                                                                     <  2:286812/149‑1 (MQ=255)
                                                gcaAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGa                                                                                                   <  4:132318/149‑1 (MQ=255)
                                                    ggAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACGACATAAAGATTGTGGGGTGAGTAc                                                                                               >  5:125039/1‑149 (MQ=255)
                                                         cGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCa                                                                                          <  8:211045/149‑1 (MQ=255)
                                                           ccATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACa                                                                                                              <  7:350717/127‑1 (MQ=255)
                                                           ccATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACa                                                                                                              >  8:350717/1‑127 (MQ=255)
                                                             aTCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATGATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGa                                                                                       >  5:111570/1‑148 (MQ=255)
                                                             aTCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATGATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGa                                                                                       >  7:244643/1‑148 (MQ=255)
                                                                  gCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATGAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTc                                                                                 <  6:197341/149‑1 (MQ=255)
                                                                  gCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAg                                                                                                  >  4:271964/1‑132 (MQ=255)
                                                                  gCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAg                                                                                                  <  3:271964/132‑1 (MQ=255)
                                                                    ttAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCt                                                                                >  1:242867/1‑148 (MQ=255)
                                                                     tAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGa                                                                              >  8:92996/1‑149 (MQ=255)
                                                                      aaGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGAc                                                                             >  4:263987/1‑149 (MQ=255)
                                                                        ggTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGa                                                                                                   >  6:283340/1‑125 (MQ=255)
                                                                        ggTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGa                                                                                                   <  5:283340/125‑1 (MQ=255)
                                                                            ttttGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTta                                                                        >  4:166239/1‑148 (MQ=255)
                                                                            ttttGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCAGTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTta                                                                        >  5:287664/1‑148 (MQ=255)
                                                                              ttGGAATGCTGATCTCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTtata                                                                      <  2:95123/148‑1 (MQ=255)
                                                                              ttGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATa                                                                                                                  <  6:139474/104‑1 (MQ=255)
                                                                              ttGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTtata                                                                      <  1:27526/148‑1 (MQ=255)
                                                                              ttGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTtata                                                                      <  8:350033/148‑1 (MQ=255)
                                                                                  aaTGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCtt                                                                 <  6:363600/149‑1 (MQ=255)
                                                                                       tGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTgc                                                            >  4:296303/1‑149 (MQ=255)
                                                                                           atCGGCTGCATCAAGAGCTGTCGGGTAGAAATAAATTGCGGAAAACAGATTGTCATAATAACCGGGTTGAAGATTTCATTTCCACCACATACACATTGGGGTGTGCGTACCGCCAGGATACGACAGACTCTTAAAGCCTTTTTGCGacc                                                        >  2:144510/1‑149 (MQ=255)
                                                                                                   cATCAAGAGCTGTCGGGTAGAAATCCACTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGc                                                <  6:287664/149‑1 (MQ=255)
                                                                                                        agagCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCgg                                              >  1:286213/1‑146 (MQ=255)
                                                                                                        agagCTGTCGGGTAGAAAAACATTGCGTAAATAAGATTCTCATAATTACCGGGATGAAGTTTAAAATAACAAAACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTGATAGACGTGTTGCGCCCCGTGCAGCGGCag                                           >  6:230313/1‑149 (MQ=255)
                                                                                                         gagCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCaga                                          <  6:212166/149‑1 (MQ=255)
                                                                                                            cTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGa                                       <  2:242867/149‑1 (MQ=255)
                                                                                                                    tAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGAt                               <  3:263987/149‑1 (MQ=255)
                                                                                                                    tAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGAt                               <  2:286213/149‑1 (MQ=255)
                                                                                                                      gAAATCCATTTCGGTAATGAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATaa                             <  6:125039/149‑1 (MQ=255)
                                                                                                                         aTCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGt                                                                                  >  7:281776/1‑93 (MQ=255)
                                                                                                                         aTCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGt                                                                                  <  8:281776/93‑1 (MQ=255)
                                                                                                                            cATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCaga                                          >  1:154862/1‑130 (MQ=255)
                                                                                                                            cATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCaga                                          <  2:154862/130‑1 (MQ=255)
                                                                                                                            cATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCAc                       >  8:317013/1‑149 (MQ=255)
                                                                                                                                  ggAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTAt                 >  5:346327/1‑149 (MQ=255)
                                                                                                                                  ggAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTAt                 >  3:243770/1‑149 (MQ=255)
                                                                                                                                     aaTCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCg               >  2:3694/1‑148 (MQ=255)
                                                                                                                                         aGCTTGTCATGATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAgg           <  8:244643/148‑1 (MQ=255)
                                                                                                                                         aGCTTGTCATGATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAgg           <  6:111570/148‑1 (MQ=255)
                                                                                                                                            ttGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAg            >  1:228178/1‑144 (MQ=255)
                                                                                                                                            ttGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAg            <  2:228178/144‑1 (MQ=255)
                                                                                                                                              gTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAGGTaaa       <  7:92996/147‑1 (MQ=255)
                                                                                                                                                cATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCa                        >  8:217891/1‑128 (MQ=255)
                                                                                                                                                cATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCa                        <  7:217891/128‑1 (MQ=255)
                                                                                                                                                 aTAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAGGTAAAGAGaa  >  7:310496/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCAGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAGGTAAAGAGAA  >  NC_000913/1581213‑1581506

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