Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,356,815 A→C 7.4% V649G (GTT→GGT)  infB ← translation initiation factor IF‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,356,8150AC 7.4% 179.4 / 7.0 68V649G (GTT→GGT) infBtranslation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base C (2/3);  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.43e-01

TGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTTTCTTCTCGTCGGCAACCA  >  NC_002947/5356667‑5356962
                                                                                                                                                     |                                                                                                                                                   
tGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGaa                                                                                                                                                     >  5:227006/1‑149 (MQ=255)
  cATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAAc                                                                                                                                                    >  3:337576/1‑148 (MQ=255)
    ttCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGAt                                                                                                                                                 <  3:339543/149‑1 (MQ=255)
       gAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGTTCTCAAGACCGATGt                                                                                                                                               <  4:270295/148‑1 (MQ=255)
       gAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGTTCTCAAGACCGATGt                                                                                                                                               >  3:270295/1‑148 (MQ=255)
          gCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                           <  4:31704/149‑1 (MQ=255)
          gCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                           >  3:31704/1‑149 (MQ=255)
          gCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                           >  8:196701/1‑149 (MQ=255)
          gCCAGCGCCAGGTTGGCATAGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                           <  5:122187/149‑1 (MQ=255)
           ccAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAgg                                                                                                                                          <  5:329814/149‑1 (MQ=255)
            cAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAggg                                                                                                                                         <  7:196701/149‑1 (MQ=255)
             aGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCCCCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGCGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGCGt                                                                                                                                        >  7:338468/1‑149 (MQ=255)
                  cAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTttct                                                                                                                                   >  3:272745/1‑149 (MQ=255)
                        cgCAGCGCACTCGGGGCTTCCCCCGACGCCACCACCCATCACGCGCACTTGCTCTTCGTGGTGTGCCCGGGCGCGGAGTGAACCCTGCAGAGCTTGCATATATACACGCACATATGTATTGATAACGATGTTGAGGGTTTtcttctctt                                                                                                                             <  8:338468/148‑1 (MQ=255)
                          cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTcc                                                                                                                           >  7:34054/1‑149 (MQ=255)
                            tCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCt                                                                                                                          >  5:348304/1‑148 (MQ=255)
                               ctctCGGTGATACCACCGACGCCACCCCCAATCACGCGCACTTGCACTTCGTCGTTTCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGAcc                                                                                                                      <  6:306970/149‑1 (MQ=255)
                                    ggTGATACCACCGACGCCACCACCACTCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATgg                                                                                                                 <  2:107969/149‑1 (MQ=255)
                                    ggTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                                  <  6:213839/148‑1 (MQ=255)
                                     gTGATACCTCCGACGCCACCCCGAATCACGCGCTCTTGCTCTCCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGACCCACGCACATCGGGCTTGAGCACGATGTTGAGGGTTTTCTTCTCTACCTGACCGATGGt                                                                                                                <  8:303320/149‑1 (MQ=255)
                                           ccaccGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGaa                                                                                                          <  4:272745/149‑1 (MQ=255)
                                             accGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGTTCTCAAGACCGAt                                                                                                                                                 <  8:43684/108‑1 (MQ=255)
                                             accGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGa                                                                                                        <  4:27496/149‑1 (MQ=255)
                                             accGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGa                                                                                                        <  2:53488/149‑1 (MQ=255)
                                              ccGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGAt                                                                                                       >  5:311812/1‑149 (MQ=255)
                                                  cGCCAACACCACTCACGCGCACTTGCACTTCCTCGTTGCCCAGGCCGCCTAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACCTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTATTCCTGACCCATGGTCTCGAAGATGttt                                                                                                   <  8:172554/149‑1 (MQ=255)
                                                    ccaccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTc                                                                                                 <  7:283979/149‑1 (MQ=255)
                                                    ccaccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTc                                                                                                 <  7:315858/149‑1 (MQ=255)
                                                      accaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTc                                                                                                 >  4:152249/1‑147 (MQ=255)
                                                                aCGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGc                                                                                     >  5:272893/1‑149 (MQ=255)
                                                                  gcgcACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGt                                                                                   <  6:227006/149‑1 (MQ=255)
                                                                  gcgcACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGt                                                                                   >  8:312761/1‑149 (MQ=255)
                                                                  gcgcACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGt                                                                                   >  2:309152/1‑149 (MQ=255)
                                                                    gcacttgcacttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGa                                                                                 <  6:311812/149‑1 (MQ=255)
                                                                      acttgcacttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAg                                                                                >  5:254324/1‑148 (MQ=255)
                                                                            acttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGc                                                                         >  4:177900/1‑149 (MQ=255)
                                                                              ttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCGCATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGc                                                                         >  4:99071/1‑147 (MQ=255)
                                                                              ttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCa                                                                       >  5:142297/1‑149 (MQ=255)
                                                                                cgtcgtTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAg                                                                      <  6:142297/148‑1 (MQ=255)
                                                                                   cgtTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                                  >  8:144347/1‑101 (MQ=255)
                                                                                   cgtTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                                  <  7:144347/101‑1 (MQ=255)
                                                                                   cgtTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCtt                                                                   >  3:326320/1‑148 (MQ=255)
                                                                                    gtTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGa                                                                 >  7:58536/1‑149 (MQ=255)
                                                                                          ccAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACgg                                                                           <  8:326908/134‑1 (MQ=255)
                                                                                          ccAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACgg                                                                           >  7:326908/1‑134 (MQ=255)
                                                                                           cAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTcgc                                                           >  1:178469/1‑149 (MQ=255)
                                                                                           cAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCATCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTcgc                                                           <  6:348304/149‑1 (MQ=255)
                                                                                           cAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCATCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTcgc                                                           <  8:58536/149‑1 (MQ=255)
                                                                                                    gAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGc                                                  <  7:312761/149‑1 (MQ=255)
                                                                                                    gAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGc                                                  <  2:178469/149‑1 (MQ=255)
                                                                                                          aCCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGa                                             <  6:254324/148‑1 (MQ=255)
                                                                                                             cTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACgg                                          >  1:316213/1‑148 (MQ=255)
                                                                                                              tGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGc                                                  >  6:346980/1‑139 (MQ=255)
                                                                                                              tGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGc                                                  <  5:346980/139‑1 (MQ=255)
                                                                                                                   cgcTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCt                                                                    >  8:45652/1‑116 (MQ=255)
                                                                                                                   cgcTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCt                                                                    <  7:45652/116‑1 (MQ=255)
                                                                                                                      ttCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTCGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACgg                                                                           <  7:260352/106‑1 (MQ=255)
                                                                                                                      ttCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTAGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACgg                                                                           >  8:260352/1‑106 (MQ=255)
                                                                                                                        ccAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCaa                                <  2:316213/147‑1 (MQ=255)
                                                                                                                                  ctcaacaTCGTTCTCAAGACCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGcc                                                                        >  1:27434/5‑97 (MQ=38)
                                                                                                                                  ctcaacaTCGTTCTCAAGACCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGcc                                                                        <  2:27434/93‑1 (MQ=38)
                                                                                                                                   cGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCttt                   <  4:326320/149‑1 (MQ=255)
                                                                                                                                    gCACATCGGGCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttc                  <  3:152249/149‑1 (MQ=255)
                                                                                                                                      agaTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttctt                <  1:309152/147‑1 (MQ=255)
                                                                                                                                      acaTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTATCGAAGAGGTTATACAGCTTGACGGCGTGAGCACGGGCAAGCTTCACCTCTCCATACTTTCCTTGACGGAACAGGGCAACTTCAAGGTATttct                 >  7:221675/1‑148 (MQ=255)
                                                                                                                                                  agaACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTATTCCAGCTTGCAGGCGTGAGTACCGGCCAGCTTCACCTCGCGATACTTCCAATGACTCAACATGTCAACTTCACGGGCTTTCGTCTCGTCGGCAAc    >  4:68870/1‑149 (MQ=255)
                                                                                                                                                   gaACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTTTCTTCTCGTCGGCAAc    >  4:35744/1‑148 (MQ=255)
                                                                                                                                                   gaACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGACGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAgg                                    <  6:308877/116‑1 (MQ=255)
                                                                                                                                                   gaACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGACGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAgg                                    >  5:308877/1‑116 (MQ=255)
                                                                                                                                                    aaCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCACCTTCACGGGCTTTCTTCTCGTCGGCAACCa  >  8:67542/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTTTCTTCTCGTCGGCAACCA  >  NC_002947/5356667‑5356962

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: