Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 495,026 A→T 5.6% intergenic (+16/‑94) recR → / → htpG gap repair protein/protein refolding molecular co‑chaperone Hsp90, Hsp70‑dependent; heat‑shock protein; ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913495,0260AT 5.6% 241.6 / 3.6 77intergenic (+16/‑94)recR/htpGgap repair protein/protein refolding molecular co‑chaperone Hsp90, Hsp70‑dependent; heat‑shock protein; ATPase
Reads supporting (aligned to +/- strand):  ref base A (32/40);  new base T (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

CGAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCAC  >  NC_000913/494878‑495173
                                                                                                                                                    |                                                                                                                                                   
cGAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGa                                                                                                                                                     >  2:259213/1‑149 (MQ=255)
cGAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGa                                                                                                                                                     <  1:370854/149‑1 (MQ=255)
                tatGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGTGATCCTGCTCTCGttt                                                                                                                                     <  6:383018/149‑1 (MQ=255)
                   gACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGa                                                                                                                                   >  2:392836/1‑148 (MQ=255)
                    aCGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaaa                                                                                                                                 <  8:270818/149‑1 (MQ=255)
                    aCGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaaa                                                                                                                                 <  2:101259/149‑1 (MQ=255)
                        ggAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAattat                                                                                                                             >  7:369472/1‑149 (MQ=255)
                        ggAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAattat                                                                                                                             >  2:50793/1‑149 (MQ=255)
                             ccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTcc                                                                                                                        <  1:392836/149‑1 (MQ=255)
                              cagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTccc                                                                                                                       <  7:1330/149‑1 (MQ=255)
                                gccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCtt                                                                                                                     <  5:392043/149‑1 (MQ=255)
                                 ccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTg                                                                                                                    <  1:449725/149‑1 (MQ=255)
                                   gAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTc                                                                                                                  <  5:421396/149‑1 (MQ=255)
                                         cTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtc                                                                                                            >  4:312373/1‑149 (MQ=255)
                                             tGGCGTTCCGGTTGGAGGCGAGCTGGAAATGGTCGAAGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtctctc                                                                                                        >  6:341160/1‑149 (MQ=255)
                                                 gTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGGGCCCCTGCTGTCGCGTGAAATTATTCGCCCTTGTCCCCCTCTCTCCcac                                                                                                    >  2:388954/1‑149 (MQ=255)
                                                        ttGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCTTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGt                                                                                             <  8:169550/149‑1 (MQ=255)
                                                          ggcggcGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGtt                                                                                            <  8:194308/148‑1 (MQ=255)
                                                              gcgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTaa                                                                                       <  1:50793/149‑1 (MQ=255)
                                                                  gCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCt                                                                                    >  2:188730/1‑148 (MQ=255)
                                                                      gAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGAt                                                                                                                                                    >  8:256168/1‑80 (MQ=255)
                                                                      gAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTaaaa                                                                               <  8:369472/149‑1 (MQ=255)
                                                                        aaTGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATg                                                                             >  2:50253/1‑149 (MQ=255)
                                                                          tGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCATGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGc                                                                           <  6:194257/149‑1 (MQ=255)
                                                                               gACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCattat                                                                      <  1:388954/149‑1 (MQ=255)
                                                                                 cGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGttttt                                                                                         <  5:262957/128‑1 (MQ=255)
                                                                                 cGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGttttt                                                                                         >  6:262957/1‑128 (MQ=255)
                                                                                  ggCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCattatt                                                                     >  2:209700/1‑147 (MQ=255)
                                                                                  ggCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGa                                                                   <  7:271443/149‑1 (MQ=255)
                                                                                   gCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAg                                                                  <  4:124103/149‑1 (MQ=255)
                                                                                    caccacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAgg                                                                 >  3:48148/1‑149 (MQ=255)
                                                                                            tGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCtt                                                                                                                     >  8:191105/1‑89 (MQ=255)
                                                                                            tGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCtt                                                                                                                     <  7:191105/89‑1 (MQ=255)
                                                                                               cacaCTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTccc                                                                                                      <  1:299740/101‑1 (MQ=255)
                                                                                               cacaCTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTccc                                                                                                      >  2:299740/1‑101 (MQ=255)
                                                                                                acacTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTccc                                                                                                                       >  6:173765/1‑83 (MQ=255)
                                                                                                acacTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTccc                                                                                                                       <  5:173765/83‑1 (MQ=255)
                                                                                                 cacTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATg                                                     >  8:307616/1‑148 (MQ=255)
                                                                                                  acTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGa                                                    >  8:161808/1‑148 (MQ=255)
                                                                                                   cTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGaaa                                                  <  5:341160/149‑1 (MQ=255)
                                                                                                   cTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGaaa                                                  <  7:307616/149‑1 (MQ=255)
                                                                                                      ccTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGa                                               >  5:188822/1‑149 (MQ=255)
                                                                                                      ccTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGa                                               >  6:79683/1‑149 (MQ=255)
                                                                                                       cTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGAc                                              <  3:312373/149‑1 (MQ=255)
                                                                                                        ttGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGAc                                              >  3:388397/1‑148 (MQ=255)
                                                                                                             gggCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGaaa                                        <  1:50253/149‑1 (MQ=255)
                                                                                                                cGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTc                                     <  3:161638/149‑1 (MQ=255)
                                                                                                                   cATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGaaa                                                  >  8:272657/1‑133 (MQ=255)
                                                                                                                   cATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGaaa                                                  <  7:272657/133‑1 (MQ=255)
                                                                                                                   cATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTg                                  <  8:228926/149‑1 (MQ=255)
                                                                                                                    aTAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTAc                                                        <  2:439858/126‑1 (MQ=255)
                                                                                                                    aTAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCACTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTAc                                                        >  1:439858/1‑126 (MQ=255)
                                                                                                                     tAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGt                                >  6:369372/1‑149 (MQ=255)
                                                                                                                        taTTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGtttt                             <  1:129535/148‑1 (MQ=255)
                                                                                                                          ttCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGtttt                             >  2:129535/1‑147 (MQ=255)
                                                                                                                          ttCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtca                           <  1:259213/149‑1 (MQ=255)
                                                                                                                              tttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAgg                                                                 >  3:170855/1‑107 (MQ=255)
                                                                                                                              tttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAgg                                                                 <  4:170855/107‑1 (MQ=255)
                                                                                                                               ttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGttt                              >  8:323457/1‑141 (MQ=255)
                                                                                                                               ttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGttt                              <  7:323457/141‑1 (MQ=255)
                                                                                                                                  tAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGaa                    <  4:48148/148‑1 (MQ=255)
                                                                                                                                   aaGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGt                  >  6:418662/1‑149 (MQ=255)
                                                                                                                                     gCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGa                <  7:161808/149‑1 (MQ=255)
                                                                                                                                      cAAGCGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGa                >  5:29097/1‑148 (MQ=255)
                                                                                                                                      cAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGa                >  7:2319/1‑148 (MQ=255)
                                                                                                                                       aagcGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtca                           <  8:46006/133‑1 (MQ=255)
                                                                                                                                       aagcGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtca                           >  7:46006/4‑136 (MQ=255)
                                                                                                                                       aaaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATgg                                                                            >  7:434874/1‑87 (MQ=255)
                                                                                                                                       aaaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATgg                                                                            <  8:434874/87‑1 (MQ=255)
                                                                                                                                       aaaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGaaa              <  2:379374/149‑1 (MQ=255)
                                                                                                                                        aaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAAc             <  6:188822/149‑1 (MQ=255)
                                                                                                                                           gagaGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGc          >  7:466309/1‑149 (MQ=255)
                                                                                                                                             gagCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGATGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                          >  4:258384/1‑131 (MQ=255)
                                                                                                                                             gagCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                          <  3:258384/131‑1 (MQ=255)
                                                                                                                                             gagCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCtt        <  4:388397/149‑1 (MQ=255)
                                                                                                                                             gagCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCtt        <  1:209700/149‑1 (MQ=255)
                                                                                                                                                  ggATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTg                                                                    <  2:421011/84‑1 (MQ=255)
                                                                                                                                                  ggATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTg                                                                    >  1:421011/1‑84 (MQ=255)
                                                                                                                                                   gATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCAc  >  7:388558/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCAC  >  NC_000913/494878‑495173

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: