Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,327 A→G 5.4% intergenic (‑55/+396) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,3270AG 5.4% 346.1 / 5.7 110intergenic (‑55/+396)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base A (56/48);  new base G (3/3);  total (59/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 111 total aligned reads displayed.

CCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAT  >  NC_000913/4094181‑4094475
                                                                                                                                                  |                                                                                                                                                    
ccgccgTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAg                                                                                                                                                    <  4:21540/149‑1 (MQ=255)
 cgccgTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAgg                                                                                                                                                   >  7:5756/1‑149 (MQ=255)
      ttGCGCTCATCCCAGCCTTTCAGGCAGGCGTCGGTGGGGGCTTTGATCATTCCCTGGACAATCCAGGACTGAGTATTGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTc                                                                                                                                              <  2:295771/149‑1 (MQ=255)
          gcTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATGCGGCGCGAgcgc                                                                                                                                          >  3:315806/1‑149 (MQ=255)
          gcTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCgcgc                                                                                                                                          >  5:52326/1‑149 (MQ=255)
          gcTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGCCGCTCgcgc                                                                                                                                          >  1:276326/1‑149 (MQ=255)
          gcTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTACCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATTCACTACCCGTAGGCCGGATAATGCTATCgcgc                                                                                                                                          >  1:283023/1‑149 (MQ=255)
                  cAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCg                                                                                                                                  >  7:375205/1‑149 (MQ=255)
                    tccTTTCAGCCAGGCGGCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGTACGGAGTAATGTTTTTCATGTTATGTGTTCCTGTAACTTCGGTGTTGTCGGATTCACTCACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                                                                                <  6:5936/148‑1 (MQ=38)
                     ccTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                                                                                >  8:62953/1‑148 (MQ=255)
                     ccTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                                                                                >  6:46747/1‑148 (MQ=255)
                       tttCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcaat                                                                                                                             >  3:210275/1‑149 (MQ=255)
                       tttCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcaat                                                                                                                             >  8:405520/1‑149 (MQ=255)
                         tCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGccta                                                                                                                             <  4:315806/147‑4 (MQ=255)
                          cAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGcctacg                                                                                                                           <  6:183391/148‑6 (MQ=38)
                           aGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcaatcaat                                                                                                                         <  6:416607/149‑1 (MQ=255)
                            gCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATGCGGCGCGAGcgccttatc                                                                                                                                    <  5:271898/137‑6 (MQ=255)
                            gCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcaatcaat                                                                                                                         <  8:304911/148‑1 (MQ=255)
                            gCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGcaatcaat                                                                                                                         >  7:304911/1‑148 (MQ=255)
                            gCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATg                                                                                                                        >  5:46745/1‑149 (MQ=255)
                             ccAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGc                                                                                                                       <  6:52326/149‑1 (MQ=255)
                               aGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCt                                                                                                                     <  1:266540/149‑1 (MQ=255)
                               aGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCt                                                                                                                     >  2:191573/1‑149 (MQ=255)
                                ggCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGATCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTg                                                                                                                    <  1:155397/149‑1 (MQ=255)
                                 gCGTCGGTGGTGGCTTTGATCATTCCCTGGACAACCCAGGACTGAGTAATGTTTTGCAGGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGAGCCGCATCCGGCAATCAATGCCTg                                                                                                                    >  4:11715/1‑148 (MQ=255)
                                 gCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTc                                                                                                                                              >  7:276043/1‑122 (MQ=255)
                                 gCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTc                                                                                                                                              <  8:276043/122‑1 (MQ=255)
                                   gTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAAGGTTTTGCAGGTTCTGTTGTCCTGAAAATTCGGTGGTGTCTGATGAACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATg                                                                                                                 <  7:326448/149‑1 (MQ=38)
                                       gtggtgGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGAc                                                                                                             >  5:365262/1‑149 (MQ=255)
                                            ggCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                        >  8:50074/1‑149 (MQ=255)
                                               tttGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc                                                                                                     >  2:372211/1‑149 (MQ=255)
                                               tttGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                    >  4:225939/1‑149 (MQ=255)
                                                  gATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGt                                                                                                   <  7:405520/148‑1 (MQ=255)
                                                  gATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                                  >  1:25056/1‑149 (MQ=255)
                                                  gATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                                  >  1:17784/1‑149 (MQ=255)
                                                   aTCCTTCCCTGGACAAACCAGGACTGAGTCATGAGTTGCATGTTCTGTGTTCCGGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGTTAAGGCGCGCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                                  <  1:159884/148‑1 (MQ=37)
                                                   aTCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                 >  8:208501/1‑149 (MQ=255)
                                                    tCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                 >  8:460209/1‑148 (MQ=255)
                                                              gACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTa                                                                                      >  2:172566/1‑149 (MQ=255)
                                                                   aCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAAc                                                                                  >  3:305587/1‑148 (MQ=255)
                                                                     cAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcg                                                                                                      >  3:65578/1‑126 (MQ=255)
                                                                     cAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcg                                                                                                      <  4:65578/126‑1 (MQ=255)
                                                                      aGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAtt                                                                              >  5:375562/1‑149 (MQ=255)
                                                                      aGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAtt                                                                              <  8:5756/149‑1 (MQ=255)
                                                                        gACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGc                                                                            <  6:365262/149‑1 (MQ=255)
                                                                        gACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACCATTGc                                                                            <  5:46747/149‑1 (MQ=255)
                                                                           tGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGa                                                                         >  1:183393/1‑149 (MQ=255)
                                                                             aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                         >  7:44893/1‑131 (MQ=255)
                                                                             aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                         <  8:44893/131‑1 (MQ=255)
                                                                             aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAtt                                                                              >  4:266722/1‑142 (MQ=255)
                                                                             aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGAc                                                                        <  6:46745/148‑1 (MQ=255)
                                                                             aGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGAc                                                                        <  3:90685/148‑1 (MQ=255)
                                                                              gTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCt                                                                      >  7:162958/1‑149 (MQ=255)
                                                                               tAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGAt                                                                                                                  <  1:442897/104‑1 (MQ=255)
                                                                               tAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGAt                                                                                                                  >  2:442897/1‑104 (MQ=255)
                                                                               tAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTg                                                                     >  2:133278/1‑149 (MQ=255)
                                                                               tAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTg                                                                     >  8:276393/1‑149 (MQ=255)
                                                                               tAATGTTTTGCATGTTCTGTGTTCCTGTAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGt                                                                    <  3:225939/149‑1 (MQ=255)
                                                                                 aTGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTa                                                                   <  8:375205/149‑1 (MQ=255)
                                                                                   gTTTTGCATGTCCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGAt                                                                                                                  <  6:170307/100‑1 (MQ=255)
                                                                                   gTTTTGCATGTCCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGAt                                                                                                                  >  5:170307/1‑100 (MQ=255)
                                                                                      ttGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGcc                                                               <  7:62953/148‑1 (MQ=255)
                                                                                            gTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAATGCGATAGCGCCTCATCCGGCAATCAATGCCTGATGAGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAg                                                        >  1:154376/1‑149 (MQ=21)
                                                                                            gTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAg                                                        >  1:154374/1‑149 (MQ=255)
                                                                                                   gtTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTtg                                                 >  5:76048/1‑149 (MQ=255)
                                                                                                       cTGTAACTTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGcc                                             <  8:349576/149‑1 (MQ=255)
                                                                                                       cTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGc                                              <  2:17784/148‑1 (MQ=255)
                                                                                                        tGTAAATTCGGTGTTGTCGGATGCACGATCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGc                                                      <  7:72235/139‑1 (MQ=255)
                                                                                                        tGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCg                                            >  4:11767/1‑149 (MQ=255)
                                                                                                         gTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCg                                            >  1:256751/1‑148 (MQ=255)
                                                                                                            aaTTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATc                                        <  8:162958/149‑1 (MQ=255)
                                                                                                              ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGc                                                                            >  2:194075/1‑111 (MQ=255)
                                                                                                              ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGc                                                                            <  1:194075/111‑1 (MQ=255)
                                                                                                              ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCg                                      >  8:52603/1‑149 (MQ=255)
                                                                                                              ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCg                                      >  1:269877/1‑149 (MQ=255)
                                                                                                              ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCGACACCTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCg                                      >  5:110162/1‑149 (MQ=255)
                                                                                                                cGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                     <  4:305587/148‑1 (MQ=255)
                                                                                                                 ggTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCAGGca                                   <  1:191573/149‑1 (MQ=255)
                                                                                                                   tgttgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaat                                 <  3:11715/149‑1 (MQ=255)
                                                                                                                       gtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatcaat                             >  8:451020/1‑149 (MQ=255)
                                                                                                                        tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACt                                                                                 <  7:369385/96‑1 (MQ=255)
                                                                                                                        tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                            >  1:346548/1‑149 (MQ=255)
                                                                                                                          ggATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc                          <  2:256751/149‑1 (MQ=255)
                                                                                                                          ggATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTAGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc                          <  6:76048/149‑1 (MQ=255)
                                                                                                                           gATGCACGAACCTTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCt                         <  1:89756/149‑1 (MQ=255)
                                                                                                                                aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCAc                                                    <  1:101108/117‑1 (MQ=255)
                                                                                                                                aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCAc                                                    >  2:101108/1‑117 (MQ=255)
                                                                                                                                aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                    >  4:186891/1‑149 (MQ=255)
                                                                                                                                aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                    >  4:109985/1‑149 (MQ=255)
                                                                                                                                aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                    >  2:370035/1‑149 (MQ=255)
                                                                                                                                    cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACg                >  7:341077/1‑149 (MQ=255)
                                                                                                                                    cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACg                >  1:124802/1‑149 (MQ=255)
                                                                                                                                     ccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCGTACAACTATTGCCGACCTGTa                                                                   <  2:426641/97‑1 (MQ=255)
                                                                                                                                     ccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCGTACAACTATTGCCGACCTGTa                                                                   >  1:426641/1‑97 (MQ=255)
                                                                                                                                       gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTg             <  1:133278/149‑1 (MQ=39)
                                                                                                                                          ggCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcg          >  4:266843/1‑149 (MQ=39)
                                                                                                                                           gCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc         <  1:172566/149‑1 (MQ=39)
                                                                                                                                              ggATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc      <  7:460209/149‑1 (MQ=39)
                                                                                                                                               gATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt     >  4:162339/1‑149 (MQ=39)
                                                                                                                                                  aaGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt  >  7:447566/1‑149 (MQ=39)
                                                                                                                                                  |                                                                                                                                                    
CCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAT  >  NC_000913/4094181‑4094475

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: