Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 962,242 T→A 6.7% D58V (GAC→GTC)  selD ← selenide, water dikinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947962,2420TA 6.7% 178.7 / 2.8 60D58V (GAC→GTC) selDselenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base A (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

ACCGAGGATGGCGATGGCCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACT  >  NC_002947/962097‑962389
                                                                                                                                                 |                                                                                                                                                    
cccgctgatggccatcgcCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTcgtcgt                                                                                                                                                   >  8:158045/7‑149 (MQ=255)
aCCGAGGATGGCGATGGCCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTcgtcgt                                                                                                                                                   >  6:21480/1‑149 (MQ=255)
       atggcgatggcCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAgcgc                                                                                                                                            <  2:154592/149‑1 (MQ=255)
          gcgatggcCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTa                                                                                                                                         <  8:263827/149‑1 (MQ=255)
                gcCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACActgc                                                                                                                                   >  1:325117/1‑149 (MQ=255)
                  cATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACActgctg                                                                                                                                 >  1:57371/1‑149 (MQ=255)
                   aTCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACtgctgc                                                                                                                                <  7:158045/149‑1 (MQ=255)
                      aGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTGGACGACGAGcgcg                                                                                                                                           <  1:94298/135‑1 (MQ=255)
                      aGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTGGACGACGAGcgcg                                                                                                                                           >  2:94298/1‑135 (MQ=255)
                      aGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCAtc                                                                                                                             >  5:238955/1‑149 (MQ=255)
                      aGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCAtc                                                                                                                             >  2:92299/1‑149 (MQ=255)
                       gCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCAtcg                                                                                                                            >  7:217602/1‑149 (MQ=255)
                        cGGATCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCAtcgt                                                                                                                           >  3:154579/1‑149 (MQ=255)
                             cgccgcCCATGGCATAGATGTCGCTGAGGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTcgcgc                                                                                                                      <  6:283931/149‑1 (MQ=255)
                                     aTGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGc                                                                                                                 >  8:41190/1‑146 (MQ=255)
                                     aTGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGc                                                                                                                 <  7:41190/146‑1 (MQ=255)
                                          aTAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGa                                                                                                        >  5:378325/1‑149 (MQ=255)
                                           tAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGtt                                                                                                             <  2:190147/143‑1 (MQ=255)
                                           tAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGAATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGtt                                                                                                             >  1:190147/1‑143 (MQ=255)
                                                    cTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGttt                                                                                              >  6:247568/1‑149 (MQ=255)
                                                     tGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGAc                                                                                                       <  7:277188/139‑1 (MQ=255)
                                                     tGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGAc                                                                                                       >  8:277188/1‑139 (MQ=255)
                                                       aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTc                                                                                             <  5:21480/147‑1 (MQ=255)
                                                       aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCgg                                                                                           <  2:325117/149‑1 (MQ=255)
                                                       aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCg                                                                                            <  2:57371/148‑1 (MQ=255)
                                                           cGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTcc                                                                                       <  5:203284/149‑1 (MQ=255)
                                                             ttGGTGGCGTCAATGCGTCCGAAGTCGTAAGGGTCATCGACGATCGGAATGAAGAAGTCGGTGGTCGAGACCACGCCGTGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAg                                                                                     <  7:251246/149‑1 (MQ=255)
                                                               ggtggCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAgg                                                                                    <  6:378325/148‑1 (MQ=255)
                                                                 tggCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGcc                                                                                 >  8:216044/1‑149 (MQ=255)
                                                                     gtcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTggg                                                                             <  1:316876/147‑1 (MQ=255)
                                                                      gcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGt                                                                            >  8:33237/1‑149 (MQ=255)
                                                                      gcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGt                                                                            >  7:269530/1‑149 (MQ=255)
                                                                        aaTGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTTCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGc                                                                          <  4:279224/149‑1 (MQ=255)
                                                                         aTGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGcc                                                                         <  7:216044/149‑1 (MQ=255)
                                                                           gCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGc                                                                       >  5:106566/1‑149 (MQ=255)
                                                                                    gTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCaa                                                              <  8:217602/149‑1 (MQ=255)
                                                                                       gTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAg                                                             <  1:92299/148‑1 (MQ=255)
                                                                                       gTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGa                                                           >  5:327131/1‑149 (MQ=255)
                                                                                        tAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCAGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGAt                                                          >  4:239482/1‑149 (MQ=255)
                                                                                           gggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCAc                                                       <  7:33237/149‑1 (MQ=255)
                                                                                            ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTTGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACg                                                      >  6:271955/1‑149 (MQ=255)
                                                                                                aTCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTc                                                    >  6:316880/1‑147 (MQ=255)
                                                                                                   gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTgg                                                                              <  6:272040/118‑1 (MQ=255)
                                                                                                   gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTgg                                                                              >  5:272040/1‑118 (MQ=255)
                                                                                                        tCGGCATGAAGAAGTCGGTGGTTGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCTAGGATCACGTCCAGCACCTTg                                          <  5:271955/149‑1 (MQ=255)
                                                                                                           gCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                       <  4:154579/149‑1 (MQ=255)
                                                                                                           gCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                       <  6:106566/149‑1 (MQ=255)
                                                                                                           gCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGGCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                       <  6:238955/149‑1 (MQ=255)
                                                                                                            cATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTggggg                                      <  6:327131/149‑1 (MQ=255)
                                                                                                               gaagaaGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                       >  7:65701/1‑145 (MQ=255)
                                                                                                                   aagTCGGTTGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCtt                               <  4:149623/149‑1 (MQ=255)
                                                                                                                             tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACcaca                      <  1:173438/148‑1 (MQ=255)
                                                                                                                             tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACca                        >  2:173438/1‑146 (MQ=255)
                                                                                                                                  aCCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg                                          <  7:127572/123‑1 (MQ=255)
                                                                                                                                  aCCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg                                          >  8:127572/1‑123 (MQ=255)
                                                                                                                                    cACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGgcgc              >  6:209017/1‑149 (MQ=255)
                                                                                                                                    cACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGgcgc              >  5:50732/1‑149 (MQ=255)
                                                                                                                                     aCGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAgg                                                                                    <  6:158671/78‑1 (MQ=255)
                                                                                                                                     aCGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAgg                                                                                    >  5:158671/1‑78 (MQ=255)
                                                                                                                                       tacgcgCTGGACGACGAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTc                                                    >  6:64107/3‑108 (MQ=255)
                                                                                                                                       tacgcgCTGGACGACGAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTc                                                    <  5:64107/106‑1 (MQ=255)
                                                                                                                                         cgcgCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATgg         <  5:247568/149‑1 (MQ=255)
                                                                                                                                                gtcgtcCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCACCATGGCTGTACt  <  5:316880/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACCGAGGATGGCGATGGCCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACT  >  NC_002947/962097‑962389

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: