Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 651,660 A→G 6.5% F335S (TTC→TCC)  dsbD‑I ← thiol/disulfide interchange protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947651,6600AG 6.5% 179.5 / 2.5 61F335S (TTC→TCC) dsbD‑Ithiol/disulfide interchange protein 1
Reads supporting (aligned to +/- strand):  ref base A (32/25);  new base G (2/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAGCCGCCGCGGGTGT  >  NC_002947/651516‑651807
                                                                                                                                                 |                                                                                                                                                   
cTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAg                                                                                                                                                 >  6:5772/1‑149 (MQ=255)
   aGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCtt                                                                                                                                              >  1:211014/1‑149 (MQ=255)
        gCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAgcg                                                                                                                                         >  7:314299/1‑149 (MQ=255)
          cGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAgcgag                                                                                                                                       >  4:156719/1‑147 (MQ=255)
                 cAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCAcc                                                                                                                                  >  5:88477/1‑147 (MQ=255)
                 cAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCa                                                                                                                                <  5:142464/149‑1 (MQ=255)
                   gCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAgg                                                                                                                              <  5:21596/149‑1 (MQ=255)
                              agTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCt                                                                                                                   >  6:279677/1‑149 (MQ=255)
                                    ccAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGTGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCt                                                                                                             <  6:270081/149‑1 (MQ=255)
                                     cAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTg                                                                                                            <  2:211014/149‑1 (MQ=255)
                                          gCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGata                                                                                                         >  6:145731/1‑147 (MQ=255)
                                          gCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGata                                                                                                         >  1:213659/1‑147 (MQ=255)
                                            gTTTTTCACCTTTTTCATCCATGGGCCCCTTTTGGGCTGCCAGGCCGGGCCACCGGTGGCTACCAGCAGCAGCGGTGAACCAATGCCCAAAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATag                                                                                                      <  1:103445/149‑1 (MQ=255)
                                            gTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGAtata                                                                                                       >  3:19915/1‑147 (MQ=255)
                                                            aTCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGat                                                                                                          <  2:310761/128‑1 (MQ=255)
                                                            aTCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAg                                                                                     <  3:156719/149‑1 (MQ=255)
                                                            aTCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAg                                                                                     >  2:222833/1‑149 (MQ=255)
                                                               cATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAgggg                                                                                  >  7:36344/1‑149 (MQ=255)
                                                               cATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAggg                                                                                   <  5:145731/148‑1 (MQ=255)
                                                                     ccGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCg                                                                             >  6:48346/1‑148 (MQ=255)
                                                                     ccGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGa                                                                            >  7:107147/1‑149 (MQ=255)
                                                                        cTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATac                                                                         <  1:215324/149‑1 (MQ=255)
                                                                           ttGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGCCGAACAGCTTTAGTGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATacac                                                                       <  1:293148/148‑1 (MQ=255)
                                                                                   ccAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGAACCGGCCAGGGGCGCCGATACACACGGAGaaa                                                              >  4:99311/1‑149 (MQ=255)
                                                                                   ccAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCCGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGaaa                                                              >  4:94221/1‑149 (MQ=255)
                                                                                        ccGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGgcgc                                                                               >  2:172919/1‑127 (MQ=255)
                                                                                        ccGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGgcgc                                                                               <  1:172919/127‑1 (MQ=255)
                                                                                        ccGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGc                                                         <  8:36344/149‑1 (MQ=255)
                                                                                          gcgcCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCaa                                                       <  5:5772/149‑1 (MQ=255)
                                                                                          gcgcCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCaa                                                       <  1:18130/149‑1 (MQ=255)
                                                                                              cACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTg                                                   <  6:88477/149‑1 (MQ=255)
                                                                                               aCCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGcc                                                                              <  8:337187/121‑1 (MQ=255)
                                                                                               aCCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGcc                                                                              >  7:337187/1‑121 (MQ=255)
                                                                                                ccGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGga                                                 <  5:48346/149‑1 (MQ=255)
                                                                                                  ggTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGgaga                                               <  4:19915/149‑1 (MQ=255)
                                                                                                      gCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGGACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATa                                                                          <  5:94220/118‑1 (MQ=255)
                                                                                                      gCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGGACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATa                                                                          <  5:94223/118‑1 (MQ=255)
                                                                                                      gCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGGACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATa                                                                          >  6:94220/1‑118 (MQ=255)
                                                                                                      gCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGGACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATa                                                                          >  6:94223/1‑118 (MQ=255)
                                                                                                      gCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGca                                            <  2:213659/148‑1 (MQ=255)
                                                                                                        tACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCacac                                         >  5:225729/1‑149 (MQ=255)
                                                                                                            atcagcatcGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCATGGCATCGCCGCTGGCGCTGATATATAGCAATTCACCGGCCAGGGGCGCCGATACACACGGAGTAACCTGCAAGCTGGAGAGCACACCGAg                                     <  7:53898/147‑1 (MQ=255)
                                                                                                            agcagcagcGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAg                                     >  1:15997/1‑149 (MQ=255)
                                                                                                            agcagcagcGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTATAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCCAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAg                                     >  7:255747/1‑149 (MQ=255)
                                                                                                                gcagcGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGca                                   >  1:158083/1‑147 (MQ=255)
                                                                                                                 cagcGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGga                                                 >  5:84786/1‑132 (MQ=255)
                                                                                                                 cagcGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGga                                                 <  6:84786/132‑1 (MQ=255)
                                                                                                                      gTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTgcg                           <  2:158083/149‑1 (MQ=255)
                                                                                                                          aCCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGAACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGcc                         >  5:317403/1‑147 (MQ=255)
                                                                                                                                gCCC‑AAACCCAGGGCGAACAGCTTTAGAGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGCGCACAGCTGCGCCCagcagc                  >  3:263331/1‑148 (MQ=255)
                                                                                                                                  cc‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAg               >  6:193146/1‑149 (MQ=255)
                                                                                                                                  cc‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAg               >  5:188668/1‑149 (MQ=255)
                                                                                                                                  cc‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAg               >  3:77796/1‑149 (MQ=255)
                                                                                                                                     aaaCCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCAGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAgcc             >  2:330771/1‑149 (MQ=255)
                                                                                                                                         ccAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAGccgccg         >  5:161423/1‑149 (MQ=255)
                                                                                                                                         ccAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAGccgccg         >  5:159249/1‑149 (MQ=255)
                                                                                                                                         ccAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAGccgccg         >  4:163792/1‑149 (MQ=255)
                                                                                                                                         ccAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAGccgccg         >  8:181365/1‑149 (MQ=255)
                                                                                                                                          cAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAGCcgccgc        <  6:317403/149‑1 (MQ=255)
                                                                                                                                             ggCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAGCCGCCGCggg     <  5:279677/149‑1 (MQ=255)
                                                                                                                                                gAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAGCCGCCGCGGgtgt  <  4:77796/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACCGTGTTCATCCATGGGCCGCTTTTGGGCAGCCAGGCCGCGCCACCGGTGGCTACCAGCAGCAGCGGTGCACCCATGCCC‑AAACCCAGGGCGAACAGCTTTAGCGCGCCACCCAGGGCATCGCCGCTGGCGCTGATATAGAGCAAGGCACCGGCCAGGGGCGCCGATACACACGGAGAAACCAGCAAGCTGGAGAGCACACCGAGCACAGCTGCGCCCAGCAGCGAGCCGCCGCGGGTGT  >  NC_002947/651516‑651807

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: ‑