Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,742,918 A→G 9.7% H5R (CAC→CGC) ‡ PP_2397 → EF hand domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,742,9180AG 9.7% 158.9 / 6.7 62H5R (CAC→CGC) ‡PP_2397EF hand domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (23/32);  new base G (4/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

TGCCAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGAGCTGC  >  NC_002947/2742772‑2743066
                                                                                                                                                  |                                                                                                                                                    
tgccAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCAcc                                                                                                                                                    >  5:28895/1‑149 (MQ=255)
  ccAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCtt                                                                                                                                                  <  4:63970/149‑1 (MQ=255)
   cAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTg                                                                                                                                                 <  3:194267/149‑1 (MQ=255)
   cAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCAGGGCCCCAACGTGCCTTg                                                                                                                                                 >  2:205699/1‑149 (MQ=255)
         aaaGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGccc                                                                                                                                           <  6:289424/149‑1 (MQ=255)
            ggCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGtt                                                                                                                                        <  5:326530/149‑1 (MQ=255)
            ggCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGtt                                                                                                                                        <  4:132759/149‑1 (MQ=255)
             gCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGccc                                                                                                                                           >  4:196021/1‑145 (MQ=255)
             gCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGccc                                                                                                                                           <  3:196021/145‑1 (MQ=255)
               ggTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGcc                                                                                                                                     >  7:242360/1‑149 (MQ=255)
                gTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCCGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGc                                                                                                                                             >  1:65015/1‑140 (MQ=255)
                gTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGc                                                                                                                                             <  2:65015/140‑1 (MQ=255)
                gTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCg                                                                                                                                    <  6:320559/149‑1 (MQ=255)
                   gcggcACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGggt                                                                                                                                 <  7:271273/149‑1 (MQ=255)
                    cggcACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGGTGCCGggtg                                                                                                                                <  4:266859/149‑1 (MQ=255)
                     ggcACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGggtgg                                                                                                                               >  6:362702/1‑149 (MQ=255)
                      gcACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGggtggt                                                                                                                              >  6:258670/1‑149 (MQ=255)
                       cACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAATGAAAAGTACAGAAAGCCGGCCAAGCTCTATAATTTATAGGACCCGCCAAAACCTCTTTAAGACTGATCACGCCTCATGGCCCCAAGGCACGTTGGGGCCCGTTGCCGGgtggtg                                                                                                                             <  6:219565/149‑1 (MQ=255)
                        aCCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGGAAGCCGGCCCAGCCCCCCCAGCCAGCCGCCCCGCCCCAACCGCGGGAAGAATTATAACGCCTAATTGCCCCAAGGCACCATGGGGCCCCGTGGCGGGTGgggt                                                                                                                            >  1:278470/1‑149 (MQ=255)
                        aCCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGgtgt                                                                                                                            <  5:114774/149‑1 (MQ=255)
                          ccGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGccat                                                                                                                                          <  4:359078/133‑3 (MQ=255)
                          ccGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTCGGGccat                                                                                                                                          >  3:359078/1‑131 (MQ=255)
                              gCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGGGGTGTGTATCg                                                                                                                      <  8:242360/149‑1 (MQ=255)
                                          ggcCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGc                                                                                                          <  8:175248/149‑1 (MQ=255)
                                          ggcCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGc                                                                                                          <  6:28895/149‑1 (MQ=255)
                                            ccccGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCt                                                                                                        <  2:278470/149‑1 (MQ=255)
                                                   aGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGtgctgc                                                                                                 <  5:258670/149‑1 (MQ=255)
                                                    gCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCGAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCa                                                                                                <  5:165515/149‑1 (MQ=255)
                                                    gCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCa                                                                                                <  8:351453/149‑1 (MQ=255)
                                                       aaGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGAc                                                                                             >  5:319791/1‑149 (MQ=255)
                                                          gAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTa                                                                                          >  2:50324/1‑149 (MQ=255)
                                                          gAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTa                                                                                          <  7:49129/149‑1 (MQ=255)
                                                           aaaaGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTAc                                                                                         <  2:282905/149‑1 (MQ=255)
                                                                 aCAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCAt                                                                                   <  6:226542/149‑1 (MQ=255)
                                                                    gCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATAc                                                                                 >  7:309925/1‑148 (MQ=255)
                                                                       aGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACaaaa                                                                             <  1:66134/149‑1 (MQ=255)
                                                                               aaGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCt                                                                     <  5:362702/149‑1 (MQ=255)
                                                                                    cGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCa                                                                                                                    <  3:267837/97‑1 (MQ=255)
                                                                                    cGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCa                                                                                                                    >  4:267837/1‑97 (MQ=255)
                                                                                                  ccccGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGa                                                   >  5:143192/1‑148 (MQ=255)
                                                                                                  ccccGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGa                                                   >  7:137405/1‑148 (MQ=255)
                                                                                                     cGCCCAAACCGCGGGCAGACTGATCACGCCTCATGGCCCCCAGCCCCCTCGGGGCCCGTTGCCGGGTGGTGGGTAGCGCACTGTCTCTGCCTGTGCTGCAGACCTACTGCCATCCAAAAGCGCGCCTGGCACGATACATCCCCGACaaa                                               <  2:264041/149‑1 (MQ=255)
                                                                                                          aaaCCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTg                                                                    <  7:323/123‑1 (MQ=255)
                                                                                                          aaaCCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTg                                                                    >  8:323/1‑123 (MQ=255)
                                                                                                            aCCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGAGGTGGGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAAc                                        <  2:257310/149‑1 (MQ=255)
                                                                                                               gCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGc                                     >  4:76137/1‑149 (MQ=255)
                                                                                                                  ggAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCt                                   >  2:274153/1‑148 (MQ=255)
                                                                                                                   gAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTaa                                 >  2:102020/1‑149 (MQ=255)
                                                                                                                    aaGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCCTACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAAt                                <  1:102020/149‑1 (MQ=255)
                                                                                                                     aGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATg                               <  1:50324/149‑1 (MQ=255)
                                                                                                                      gACTGATGACGCCTCATGGCCCCAAGGCAACTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTCTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATgg                              >  2:83882/1‑149 (MQ=255)
                                                                                                                           aTCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCa                         >  5:228416/1‑149 (MQ=255)
                                                                                                                            tCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCAc                        >  1:281882/1‑149 (MQ=255)
                                                                                                                                gCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGATCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTc                    >  5:82452/1‑149 (MQ=255)
                                                                                                                                gCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTc                    >  1:12401/1‑149 (MQ=255)
                                                                                                                                gCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACCATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTc                    <  1:274153/149‑1 (MQ=255)
                                                                                                                                     aTGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCaa                 <  6:143192/147‑1 (MQ=255)
                                                                                                                                      tGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACttt                     >  6:109983/1‑142 (MQ=255)
                                                                                                                                      tGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACttt                     <  5:109983/142‑1 (MQ=255)
                                                                                                                                      gggCCCCAAGGTGCCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAAccc              >  7:207120/2‑149 (MQ=255)
                                                                                                                                           ccAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTg         <  1:157744/149‑1 (MQ=255)
                                                                                                                                                gCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGAGCt    >  4:254133/1‑149 (MQ=255)
                                                                                                                                                  aCCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGGGCTGc  <  8:137405/149‑1 (MQ=255)
                                                                                                                                                  aCCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGAGCTg   >  1:132775/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGCCAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGAGCTGC  >  NC_002947/2742772‑2743066

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