Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,008,719 C→G 5.9% A148G (GCC→GGC)  yedE → UPF0394 family sulphur transport domain‑containing inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,008,7190CG 5.9% 226.1 / 5.8 68A148G (GCC→GGC) yedEUPF0394 family sulphur transport domain‑containing inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base G (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

ACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCAGGCGCGGCGTCGTTT  >  NC_000913/2008572‑2008866
                                                                                                                                                   |                                                                                                                                                   
acaaTGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGcc                                                                                                                                                    >  5:280926/1‑149 (MQ=255)
 caaTGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCt                                                                                                                                                   <  1:203281/149‑1 (MQ=255)
        aaaCTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCtt                                                                                                                                            >  4:232908/1‑149 (MQ=255)
         aaCTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCttt                                                                                                                                           <  6:232182/149‑1 (MQ=255)
         aaCTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCttt                                                                                                                                           <  6:280926/149‑1 (MQ=255)
          aCTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTg                                                                                                                                          <  5:297365/149‑1 (MQ=255)
                    ccGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCa                                                                                                                                 <  1:42131/148‑1 (MQ=255)
                              gTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGt                                                                                                                      <  4:162050/149‑1 (MQ=255)
                                 tCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGACGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCa                                                                                                                   >  1:353526/1‑149 (MQ=255)
                                   cGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCACTGGGCTGTACCCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTGCTTTGCCATCGCCACTGCCATTGGTTCATg                                                                                                                 <  2:427346/149‑1 (MQ=255)
                                   cGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCGGGTTCTTTGCCATCGCCACTGCCATTGGTTCATg                                                                                                                 <  2:453890/149‑1 (MQ=255)
                                      aTTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGtt                                                                                                              >  4:66489/1‑149 (MQ=255)
                                       ttATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGttt                                                                                                             >  7:102511/1‑149 (MQ=255)
                                       ttATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGttt                                                                                                             >  3:8404/1‑149 (MQ=255)
                                         aTGCCGTCCATCATTGGCGGCATTATCGCCGGTGTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGACGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTg                                                                                                            <  2:353526/148‑1 (MQ=255)
                                          tGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGgc                                                                                                          >  8:82292/1‑149 (MQ=255)
                                                 catcatTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCt                                                                                                   >  5:10771/1‑149 (MQ=255)
                                                         gcggcATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGgcgcgcg                                                                                                     <  7:97379/139‑1 (MQ=255)
                                                         gcggcATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGgcgcgcg                                                                                                     >  8:97379/1‑139 (MQ=255)
                                                              aTTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGcc                                                                                      >  7:307714/1‑149 (MQ=255)
                                                               ttATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGcc                                                                                      >  5:239617/1‑148 (MQ=255)
                                                                    gCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCttt                                                                                                                                           >  2:445282/1‑90 (MQ=255)
                                                                    gCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCttt                                                                                                                                           <  1:445282/90‑1 (MQ=255)
                                                                      cGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTcc                                                                              <  8:199081/149‑1 (MQ=255)
                                                                      cGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTcc                                                                              <  6:239617/149‑1 (MQ=255)
                                                                       ggTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGTCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCg                                                                             <  2:203080/149‑1 (MQ=255)
                                                                         ttttGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTa                                                                           >  7:84229/1‑149 (MQ=255)
                                                                             ggCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTcc                                                                       >  8:59832/1‑149 (MQ=255)
                                                                              tcgcgcgTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCATGTTTCTTTGACATCGCCACTGCCATTGGTTCATGGTTTg                                                                                                            <  3:282503/110‑1 (MQ=255)
                                                                              gcgcgcGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTg                                                                                                            >  4:282503/1‑111 (MQ=255)
                                                                                    gTCTGGCAATGGGCTGTACCCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCGTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGaa                                                                 >  2:124254/1‑148 (MQ=255)
                                                                                     tCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCttt                                                                                                                                           >  4:185256/1‑73 (MQ=255)
                                                                                     tCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCttt                                                                                                                                           <  3:185256/73‑1 (MQ=255)
                                                                                     tCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGaaaa                                                               <  7:82292/149‑1 (MQ=255)
                                                                                      cgggCCATGGGCTGGAACCGGGCGGCGTTCGTTACCGGTATTCCTCCGTTCGCGCGGCATGCCGGGGTCCTTGCCATCGCCACTGCCCTTGGTTCCTGGTTTGGGGCGCGCGTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAAt                                                              <  3:203766/147‑1 (MQ=255)
                                                                                       tGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGt                                                                                                               <  2:89096/99‑1 (MQ=255)
                                                                                       tGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATAGGTTCATGGt                                                                                                               >  1:89096/1‑99 (MQ=255)
                                                                                        ggCAATGGGCTGTAACCCGGCGGCGTTCTTTAGCTGTATTCCTCAGTTCTCGCTGCTTGCATTGTTCTTTGCCATCTCTTCTGCCATTGGTTCAGGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGc                                                            <  4:8404/149‑1 (MQ=255)
                                                                                         gCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATg                                                             <  7:59832/147‑1 (MQ=255)
                                                                                            aTGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGc                                                            >  7:92971/1‑145 (MQ=255)
                                                                                            aTGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGc                                                            <  8:92971/145‑1 (MQ=255)
                                                                                             tGGGCTGTAACCTGGCGGCGTTCTTGACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCAGTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCGTCCGTATTCCCGTGAAAATGCAGaac                                                       >  2:406174/1‑148 (MQ=255)
                                                                                                     aaCCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCt                                                 >  6:266794/1‑147 (MQ=255)
                                                                                                     aaCCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCt                                                 <  5:266794/147‑1 (MQ=255)
                                                                                                         tGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGt                                                                    >  5:280804/1‑124 (MQ=255)
                                                                                                         tGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGt                                                                    <  6:280804/124‑1 (MQ=255)
                                                                                                             ggcgTTCTTTACCGGGATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCAc                                       <  1:84244/149‑1 (MQ=255)
                                                                                                                    tttACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGAc                                >  1:417329/1‑149 (MQ=255)
                                                                                                                    tttACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCAGAGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGAc                                >  6:27304/1‑149 (MQ=255)
                                                                                                                      tACCGGTATTCCTCAGTTCTCGCTGCATGGCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACg                               >  8:265701/1‑148 (MQ=255)
                                                                                                                      tACCGGTATTCCTCAGTTCTCGCTGCATGGCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACg                               >  2:25595/1‑148 (MQ=255)
                                                                                                                              ttCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAAcc                       <  3:66489/148‑1 (MQ=255)
                                                                                                                               tCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCTGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCg                      >  4:47916/1‑148 (MQ=255)
                                                                                                                                 ctctggtctcGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGGTTGGCGCGCGCGTTACCCTGCTGCCCAGCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTATGCCGCCTGACCGCTGACGCAAAAACCGGAt                   <  1:406174/143‑1 (MQ=255)
                                                                                                                                 cTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGa                    <  3:232908/148‑1 (MQ=255)
                                                                                                                                 cTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGAt                   <  2:417329/149‑1 (MQ=255)
                                                                                                                                 cTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGAt                   <  8:307714/149‑1 (MQ=255)
                                                                                                                                    aGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCTGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCa                 <  3:47916/148‑1 (MQ=255)
                                                                                                                                    aGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATc                  <  1:394543/147‑1 (MQ=255)
                                                                                                                                    aGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCAg                <  6:10771/149‑1 (MQ=255)
                                                                                                                                      ttCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCAGgc              >  5:88372/1‑149 (MQ=255)
                                                                                                                                        ctcGCTGCGTGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCctgct                                                                                         >  1:91083/1‑72 (MQ=255)
                                                                                                                                        ctcGCTGCGTGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCctgct                                                                                         <  2:91083/72‑1 (MQ=255)
                                                                                                                                         tcGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCAGGcgcg           >  8:88479/1‑149 (MQ=255)
                                                                                                                                           gCTGCATGGCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCAGGcgcg           <  7:265701/147‑1 (MQ=255)
                                                                                                                                           gCTGCATGGCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCAGGcgcg           <  1:25595/147‑1 (MQ=255)
                                                                                                                                           gCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCAGGcgcg           <  8:102511/147‑1 (MQ=255)
                                                                                                                                             tGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCAGGCGCGGcgt       >  4:138672/1‑149 (MQ=255)
                                                                                                                                                  gCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATAAGGCGCGGCGTCGttt  >  1:197715/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTTGCCATCGCCACTGCCATTGGTTCATGGTTTGGCGCGCGCTTTACCCTGCTGCCCATCTTCCGTATTCCCGTGAAAATGCAGAAAGTTTCTGCCGCCTCACCGCTGACGCAAAAACCGGATCAGGCGCGGCGTCGTTT  >  NC_000913/2008572‑2008866

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