Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,207 T→G 100% V53V (GTT→GTG)  stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,2070TG97.1% 4964.5 / 62.1 1734V53V (GTT→GTG) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (24/31);  new base G (832/845);  total (857/877)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 1736 total aligned reads displayed.

ATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATG‑CCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAG  >  NC_000913/1429061‑1429351
                                                                                                                                                  |                                                                                                                                                 
 tttCAGGATTCCTTAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTAATCTTAAAGCCAGACGTAACAGCACCACGGTGGTTGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAg                                                                                                                                                <  7:87281/149‑1 (MQ=255)
    cAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTa                                                                                                                                              <  4:365250/148‑1 (MQ=255)
    cAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTa                                                                                                                                              <  6:382387/148‑1 (MQ=255)
     aGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTAc                                                                                                                                             <  4:121264/148‑1 (MQ=255)
     aGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTAAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACg                                                                                                                                            <  4:92917/149‑1 (MQ=255)
       gagtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGt                                                                                                                                          <  4:154942/144‑1 (MQ=255)
       gagtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGt                                                                                                                                          >  4:272046/6‑149 (MQ=255)
         gtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCa                                                                                                                                        <  8:190056/146‑1 (MQ=255)
            cTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTACCAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGGCAGTa                                                                                                                                     <  1:311827/149‑1 (MQ=255)
            cTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGt                                                                                                                                      <  5:310248/148‑1 (MQ=255)
            cTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTa                                                                                                                                     >  8:73033/1‑149 (MQ=255)
             tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCCCGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTa                                                                                                                                     >  1:157225/1‑148 (MQ=255)
             tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGa                                                                                                                                                 >  5:203668/1‑136 (MQ=255)
             tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTa                                                                                                                                     >  2:144000/1‑148 (MQ=255)
             tGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTa                                                                                                                                     >  5:52761/1‑148 (MQ=255)
              gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTa                                                                                                                                     >  8:300678/1‑147 (MQ=255)
               aaaGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAg                                                                                                                                  <  6:192270/149‑1 (MQ=255)
                 aGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTTAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGAGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtg                                                                                                                                <  2:279336/149‑3 (MQ=255)
                   aCGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACGCGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGTGTACGGTCAGTACAgtgt                                                                                                                               >  6:264048/1‑148 (MQ=255)
                   aCGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTTGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtc                                                                                                                              <  6:158767/149‑2 (MQ=255)
                   aCGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTTGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtc                                                                                                                              <  1:275642/149‑2 (MQ=255)
                     ggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTTGGGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCAt                                                                                                                            <  6:114885/149‑1 (MQ=255)
                      gCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCAtc                                                                                                                           <  5:36050/149‑2 (MQ=255)
                          aGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTgc                                                                                                                       <  8:269142/149‑2 (MQ=255)
                          aGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTg                                                                                                                        >  7:332710/1‑148 (MQ=255)
                          aGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGGGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTgc                                                                                                                       <  8:368811/149‑2 (MQ=255)
                             aaaaCCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTgcagg                                                                                                                    <  7:268335/149‑5 (MQ=255)
                              aaaCCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTgcaggt                                                                                                                   <  5:127314/149‑6 (MQ=255)
                              aaaCCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTgcaggt                                                                                                                   <  1:196751/149‑6 (MQ=255)
                                                       aGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCCGGTTGACGGTTTTCCACCATCGCACG‑cc                                                                                         >  2:221257/1‑149 (MQ=255)
                                                       aGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑cc                                                                                         >  3:182795/1‑149 (MQ=255)
                                                          tGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCgg                                                                                       >  4:311021/1‑148 (MQ=255)
                                                            aaaGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGa                                                                                     <  7:73033/148‑1 (MQ=255)
                                                              agccagaCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGAcc                                                                                   >  4:58805/7‑148 (MQ=255)
                                                                ccagaCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCAtc                                                                                <  7:306443/145‑2 (MQ=255)
                                                                 cagaCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCAtca                                                                               >  7:244541/4‑147 (MQ=255)
                                                                   gaCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCAcc                                                                             >  2:163489/2‑149 (MQ=255)
                                                                    aCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCg                                                                            >  7:362634/1‑149 (MQ=255)
                                                                    aCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCg                                                                            >  3:9345/1‑149 (MQ=255)
                                                                     cGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCTGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCgt                                                                           <  6:143225/149‑1 (MQ=255)
                                                                     cGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCgt                                                                           <  8:149778/149‑1 (MQ=255)
                                                                      gTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAg                                                                                                                                                >  1:132501/1‑80 (MQ=255)
                                                                      gTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCCCCgtg                                                                          >  5:257835/1‑149 (MQ=255)
                                                                      gTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCgtg                                                                          >  5:171238/1‑149 (MQ=255)
                                                                       tAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCgtgt                                                                         <  8:21812/149‑1 (MQ=255)
                                                                       tAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCgtgt                                                                         <  5:333969/149‑1 (MQ=255)
                                                                       tAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCgtg                                                                          >  3:288421/1‑148 (MQ=255)
                                                                         aCAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACa                                                                                                                                   <  6:220816/90‑1 (MQ=255)
                                                                         aCAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCgtgt                                                                         <  1:2385/147‑1 (MQ=255)
                                                                         aCAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTAt                                                                       >  1:101534/1‑149 (MQ=255)
                                                                         aCAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGGTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTAt                                                                       <  3:124104/149‑1 (MQ=255)
                                                                           aGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATg                                                                      >  8:329663/1‑148 (MQ=255)
                                                                           aGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGa                                                                     <  1:395459/149‑1 (MQ=255)
                                                                            gCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCATGCAGGTTGACGGTTTTCCACCATAGCACG‑ACGGGACCATCACCGTGTATg                                                                      >  2:239557/1‑147 (MQ=255)
                                                                               ccacGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGAt                                                                 <  3:353692/149‑1 (MQ=255)
                                                                               ccacGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGAt                                                                 >  5:146644/1‑149 (MQ=255)
                                                                                cacGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtg                                                                                                                                >  7:218372/1‑84 (MQ=255)
                                                                                cacGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGAtt                                                                >  3:234351/1‑149 (MQ=255)
                                                                                cacGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGAtt                                                                <  2:385158/149‑1 (MQ=255)
                                                                                cacGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGAt                                                                 >  7:144685/1‑148 (MQ=255)
                                                                                 acGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGGTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATtc                                                               <  1:331529/149‑1 (MQ=255)
                                                                                   ggtggtggtGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATtcac                                                             <  6:336447/149‑3 (MQ=255)
                                                                                          gtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAAccgg                                                       <  5:80961/148‑3 (MQ=255)
                                                                                          gtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAAccgg                                                       <  1:220099/148‑3 (MQ=255)
                                                                                          gtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGCTCAGTAAAGTGTCATCCTGCAGGTTGACGGTTTTCCAACTTCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAAcgggg                                                      >  8:347149/1‑145 (MQ=255)
                                                                                          gtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGAGGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAAccgg                                                       <  6:304956/148‑3 (MQ=255)
                                                                                          gtgAACACGGGGGGCGCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGGCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAAccgg                                                       <  5:287034/148‑3 (MQ=255)
                                                                                           tgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTTTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCgggg                                                     <  4:180803/149‑2 (MQ=255)
                                                                                           tgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCgggg                                                     <  5:381548/149‑2 (MQ=255)
                                                                                             aaCACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCgggga                                                    <  3:220663/148‑3 (MQ=255)
                                                                                              acacGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACg                                                  >  2:366496/1‑149 (MQ=255)
                                                                                              acacGGTGGGCTCAGAGAATCCGCATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGTTTTTCCACCATCGCAAG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACg                                                  >  2:105732/1‑149 (MQ=255)
                                                                                                acGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGCCCATCACCGTGTATGAAGATTCACAACCGGGGACGCt                                                >  6:61853/1‑149 (MQ=255)
                                                                                                acGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCt                                                <  3:163842/149‑1 (MQ=255)
                                                                                                   gtGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGaa                                             >  3:71818/2‑149 (MQ=255)
                                                                                                    tGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAAt                                            <  2:198042/149‑1 (MQ=255)
                                                                                                         tCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACGACCGGGGACGCTGAATGAttt                                       >  2:345177/1‑149 (MQ=255)
                                                                                                         tCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAttt                                       >  7:48544/1‑149 (MQ=255)
                                                                                                         tCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAttt                                       >  8:10830/1‑149 (MQ=255)
                                                                                                           agagAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAtttt                                      >  2:320649/1‑149 (MQ=255)
                                                                                                           agagAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTtc                                     >  7:117350/1‑149 (MQ=255)
                                                                                                              gaaTCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTtctct                                  <  3:104164/148‑2 (MQ=255)
                                                                                                               aaTCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTtctctg                                 >  1:365983/1‑147 (MQ=255)
                                                                                                                 tCCGGATGAAGCCGGGCGTTACAGCTTGGATGTTGAGTACGTTCAGTTCAGTGTCATACTGCAGGTTGACGTTTTTCCACCATGCCACT‑CCGGGACCATGCCCGTGTATGAAGATTCACATCCGGGGACGCTGAATGATTTTCTCTGt                                <  6:89090/148‑1 (MQ=255)
                                                                                                                 tCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCtgtg                               >  8:251676/1‑149 (MQ=255)
                                                                                                                  ccGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCCCCATCGCACG‑CCGGGACCATCACCGTGTATGATGATTCACCACCGGGGACGCTGAATGATTTTCTATGTGc                              >  4:225342/1‑149 (MQ=255)
                                                                                                                   cGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGc                              <  8:397767/148‑1 (MQ=255)
                                                                                                                    ggATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACAGGGGACGCTGAATGATTTTCTCTGTGCCa                            >  3:45091/1‑149 (MQ=255)
                                                                                                                     gATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGc                              >  1:62361/1‑146 (MQ=255)
                                                                                                                     gATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGc                              <  2:62361/146‑1 (MQ=255)
                                                                                                                     gATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCAt                           <  7:178806/149‑1 (MQ=255)
                                                                                                                        gAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGAc                        >  6:35892/1‑149 (MQ=255)
                                                                                                                         aaGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACg                       >  3:319545/1‑149 (MQ=255)
                                                                                                                           gCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACg                       >  5:226409/1‑147 (MQ=255)
                                                                                                                            ccGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACggag                    >  4:248373/1‑149 (MQ=255)
                                                                                                                                 cGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGgatga               <  1:84408/149‑1 (MQ=255)
                                                                                                                                      cAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGccc          <  3:258683/149‑1 (MQ=255)
                                                                                                                                          aTGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGcc      >  7:383435/1‑149 (MQ=255)
                                                                                                                                           tGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCg     <  4:335317/149‑1 (MQ=255)
                                                                                                                                            ggATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACG‑CCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCccggccgg    >  1:352391/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATG‑CCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAG  >  NC_000913/1429061‑1429351

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