Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 712,373 T→C 6.1% F724S (TTC→TCC)  ileS → isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947712,3730TC 6.1% 191.7 / 2.1 65F724S (TTC→TCC) ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base C (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GACCCGGCCCGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAGGCGCTGGTGCGCTGGATCGCGCCGATCCT  >  NC_002947/712228‑712520
                                                                                                                                                  |                                                                                                                                                   
gACCCGGCCCGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAGGGTGCACAACTTCt                                                                                                                                                   <  1:253644/149‑1 (MQ=255)
gACCCGGCCCGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCt                                                                                                                                                   <  2:220364/148‑1 (MQ=255)
gACCCGGCCCGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCt                                                                                                                                                   <  7:259099/148‑1 (MQ=255)
 aCCCGGCCCGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGc                                                                                                                                                 >  2:43290/1‑149 (MQ=255)
  cccgccccgcGACCTGATGGCACCGGAAGACATGCTGGCGCTGCACCGCTGGGCCGTGGACCGTACCCTCCTCCTGCAACGCGAACTGGAATAGCACTACAGCGAGTACAGATTCTGGAACGTCTACTCCAAGGTGAACAACttatacg                                                                                                                                                >  2:316028/1‑144 (MQ=255)
  cccGGCCCGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCg                                                                                                                                                >  6:30701/1‑149 (MQ=255)
    cGGCCCGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTa                                                                                                                                              <  1:43290/149‑1 (MQ=255)
    cGGCCCGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTa                                                                                                                                              >  7:205731/1‑149 (MQ=255)
       cccGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAg                                                                                                                                           <  6:240548/149‑1 (MQ=255)
         cgcgACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAg                                                                                                                                           >  3:108068/1‑147 (MQ=255)
               ctgctgGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTggg                                                                                                                                   >  7:132043/1‑149 (MQ=255)
                     gCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGggcggc                                                                                                                               <  6:7203/147‑1 (MQ=255)
                       aCCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGc                                                                                                                                                 <  6:208706/127‑1 (MQ=255)
                       aCCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGc                                                                                                                                                 <  8:200663/127‑1 (MQ=255)
                       aCCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGc                                                                                                                                                 >  7:200663/1‑127 (MQ=255)
                       aCCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGc                                                                                                                                                 >  5:208706/1‑127 (MQ=255)
                       aCCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCt                                                                                                                           <  1:316028/149‑1 (MQ=255)
                        ccGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACATCTGCGTACAGGAGCTGGGCGGCTTCTa                                                                                                                          <  1:195274/149‑1 (MQ=255)
                         cGGAAGAGATGCTGGCGCTGGCCCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGATCTGGAAGAGCACTACAGCGAGTACGGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGTGGCTTCTAt                                                                                                                         <  4:45641/149‑1 (MQ=255)
                            aaGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTAtctc                                                                                                                      <  2:113481/149‑1 (MQ=255)
                            aaGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTAtctc                                                                                                                      <  1:323302/149‑1 (MQ=255)
                              gACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGa                                                                                                                    <  1:188596/149‑1 (MQ=255)
                                 aTGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGAca                                                                                                                  >  3:235017/1‑148 (MQ=255)
                                             gACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCg                                                                                                     >  4:323987/1‑149 (MQ=255)
                                                   tGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTa                                                                                               >  6:122812/1‑149 (MQ=255)
                                                   tGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTa                                                                                               >  5:104357/1‑149 (MQ=255)
                                                     ggCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTaca                                                                                             >  6:178233/1‑149 (MQ=255)
                                                     ggCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTa                                                                                               <  4:108068/147‑1 (MQ=255)
                                                           ggACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCt                                                                                                                                                   >  6:93740/1‑89 (MQ=255)
                                                           ggACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCt                                                                                                                                                   <  5:93740/89‑1 (MQ=255)
                                                              ccGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTa                                                                                                                          <  6:51720/111‑1 (MQ=255)
                                                              ccGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTa                                                                                                                          >  5:51720/1‑111 (MQ=255)
                                                                    ccTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCaa                                                                                >  4:306944/1‑147 (MQ=255)
                                                                     ctgctgctgCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAg                                                                                                          >  8:91756/1‑120 (MQ=255)
                                                                        ctgctgcAACGCGAACTGGAAGAGTACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGgc                                                                                     >  2:76696/1‑138 (MQ=255)
                                                                        ctgctgcAACGCGAACTGGAAGAGTACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGCCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGgc                                                                                     <  1:76696/138‑1 (MQ=255)
                                                                           ctgcAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTcgc                                                                       <  3:323987/149‑1 (MQ=255)
                                                                           ctgcAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTcg                                                                        <  5:178233/148‑1 (MQ=255)
                                                                             gcAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAg                                                                             >  3:70581/1‑141 (MQ=255)
                                                                             gcAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAg                                                                             <  4:70581/141‑1 (MQ=255)
                                                                                aCGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTcgcccgcc                                                                  >  7:54816/1‑149 (MQ=255)
                                                                                aCGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGACAGTACACCACCGGCGCCAACAGTGTcgcccgcc                                                                  >  8:163737/1‑149 (MQ=255)
                                                                                 cgcgAACTGGAATAGCACTACAGCTAGTACCGTTTATGGAACTTCTACTCCAAGGTGCACAACTCCTGCGTACAGGAGCTGGGCGGCTTCTATCTCTACATCATCAAGGACCg                                                                                                     <  7:183230/113‑1 (MQ=255)
                                                                                 cgcgAACTGGAAGAGCACTACAGCGAGTTCCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTCCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCg                                                                                                     <  7:178297/113‑1 (MQ=255)
                                                                                 cgcgAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGAGCACAACTCCTGCGTACAGGAGCTCGGCGGCTTCTATATCGACATCATCAACGACCg                                                                                                     >  8:178297/1‑113 (MQ=255)
                                                                                 cgcgAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAAAGTCTACTACAAGGTGCACAACTCCTGCGTACAGGAGCTGGGCGGCTTATATATAGACATCATCAAGGAcag                                                                                                     >  8:183230/1‑111 (MQ=255)
                                                                                  gcgAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACAtcatca                                                                                                            <  3:151553/105‑1 (MQ=255)
                                                                                  gcgAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACAtcatca                                                                                                            >  4:151553/1‑105 (MQ=255)
                                                                                      aCTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCt                                                            <  8:205731/149‑1 (MQ=255)
                                                                                          gAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTca                                                             >  5:54921/1‑143 (MQ=255)
                                                                                          gAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTca                                                             <  6:54921/144‑2 (MQ=255)
                                                                                          gAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGc                                                          >  5:320114/1‑147 (MQ=255)
                                                                                          gAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGc                                                          <  8:14700/147‑1 (MQ=255)
                                                                                                 aCTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTgcg                                                 >  3:47111/1‑149 (MQ=255)
                                                                                                    tACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCt                                               >  3:138399/1‑149 (MQ=255)
                                                                                                      cAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTg                                              >  4:85400/1‑148 (MQ=255)
                                                                                                      cAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGt                                             <  4:138399/149‑1 (MQ=255)
                                                                                                       aGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTa                                            <  6:104357/149‑1 (MQ=255)
                                                                                                        gCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTAc                                           <  4:235017/149‑1 (MQ=255)
                                                                                                                cGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAg                                   >  4:272435/1‑149 (MQ=255)
                                                                                                                    tCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGa                                <  3:85400/148‑1 (MQ=255)
                                                                                                                     cTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAg                               <  5:122812/148‑1 (MQ=255)
                                                                                                                          aCGTCTACTCCAAGGGGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAGGCGCTg                         <  5:281578/149‑1 (MQ=255)
                                                                                                                                           cacaACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAGGCGCTGGTGCGCTGGATCgcgc         <  4:47111/148‑1 (MQ=255)
                                                                                                                                            acaACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAGGCGCTGGTGCGCTGGATCGCGcc        >  1:237577/1‑148 (MQ=255)
                                                                                                                                              aaCTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAGGCGCTGGTGCGCTGGATCGCGCCGAt     >  6:11774/1‑149 (MQ=255)
                                                                                                                                                 ttCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAGGCGCTGGTGCGCTGGATCGCGCCGATCCt  >  8:92865/1‑149 (MQ=255)
                                                                                                                                                  tCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAGGCGCTGGTg                      >  1:80306/1‑128 (MQ=255)
                                                                                                                                                  tCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAGGCGCTGGTg                      <  2:80306/128‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACCCGGCCCGCGACCTGCTGGCACCGGAAGACATGCTGGCGCTGGACCGCTGGGCCGTGGACCGTACCCTGCTGCTGCAACGCGAACTGGAAGAGCACTACAGCGAGTACCGTTTCTGGAACGTCTACTCCAAGGTGCACAACTTCTGCGTACAGGAGCTGGGCGGCTTCTATCTCGACATCATCAAGGACCGCCAGTACACCACCGGCGCCAACAGTGTCGCCCGCCGTTCCTGCCAGACTGCGCTGTACCACATCAGCGAGGCGCTGGTGCGCTGGATCGCGCCGATCCT  >  NC_002947/712228‑712520

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: