Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,742,917 C→T 9.7% H5Y (CAC→TAC) ‡ 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,9170CT 9.7% 166.2 / 6.9 62H5Y (CAC→TAC) ‡PP_2397EF hand domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (23/33);  new base T (4/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

TGCCAGCGAAAAGGCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGAGCT  >  NC_002947/2742772‑2743064
                                                                                                                                                  |                                                                                                                                                   
tgccAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCCAGGccc                                                                                                                                                   >  1:257310/1‑147 (MQ=255)
tgccAGCGAAAAGGCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCAcc                                                                                                                                                  >  5:28895/1‑149 (MQ=255)
  ccAGCGAAAAGGCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCtt                                                                                                                                                <  4:63970/149‑1 (MQ=255)
   cAGCGAAAAGGCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTg                                                                                                                                               <  3:194267/149‑1 (MQ=255)
   cAGCGAAAAGGCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCAGGGCCCCAACGTGCCTTg                                                                                                                                               >  2:205699/1‑149 (MQ=255)
         aaaGGCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGccc                                                                                                                                         <  6:289424/149‑1 (MQ=255)
            ggCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGtt                                                                                                                                      <  4:132759/149‑1 (MQ=255)
            ggCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGtt                                                                                                                                      <  5:326530/149‑1 (MQ=255)
             gCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGccc                                                                                                                                         <  3:196021/145‑1 (MQ=255)
             gCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGccc                                                                                                                                         >  4:196021/1‑145 (MQ=255)
               ggTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGcc                                                                                                                                   >  7:242360/1‑149 (MQ=255)
                gTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCCGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGc                                                                                                                                           >  1:65015/1‑140 (MQ=255)
                gTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGc                                                                                                                                           <  2:65015/140‑1 (MQ=255)
                gTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCg                                                                                                                                  <  6:320559/149‑1 (MQ=255)
                   gcggcACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGggt                                                                                                                               <  7:271273/149‑1 (MQ=255)
                    cggcACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGGTGCCGggtg                                                                                                                              <  4:266859/149‑1 (MQ=255)
                     ggcACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGggtgg                                                                                                                             >  6:362702/1‑149 (MQ=255)
                      gcACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGggtggt                                                                                                                            >  6:258670/1‑149 (MQ=255)
                       cACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAATGAAAAGTACAGAAAGCCGGCCAAGCTCTATAATTTATAGGACCCGCCAAAACCTCTTTAAGACTGATCACGCCTCATGGCCCCAAGGCACGTTGGGGCCCGTTGCCGGgtggtg                                                                                                                           <  6:219565/149‑1 (MQ=255)
                        aCCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGGAAGCCGGCCCAGCCCCCCCAGCCAGCCGCCCCGCCCCAACCGCGGGAAGAATTATAACGCCTAATTGCCCCAAGGCACCATGGGGCCCCGTGGCGGGTGgggt                                                                                                                          >  1:278470/1‑149 (MQ=255)
                        aCCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGgtgt                                                                                                                          <  5:114774/149‑1 (MQ=255)
                          ccGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGccat                                                                                                                                        <  4:359078/133‑3 (MQ=255)
                          ccGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTCGGGccat                                                                                                                                        >  3:359078/1‑131 (MQ=255)
                              gCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGGGGTGTGTATCg                                                                                                                    <  8:242360/149‑1 (MQ=255)
                                           ggcCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGc                                                                                                        <  6:28895/149‑1 (MQ=255)
                                           ggcCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGc                                                                                                        <  8:175248/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                                                                                                                  >  4:267837/1‑97 (MQ=255)
                                                                                     cGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCa                                                                                                                  <  3:267837/97‑1 (MQ=255)
                                                                                                   ccccGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGa                                                 >  7:137405/1‑148 (MQ=255)
                                                                                                   ccccGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGa                                                 >  5:143192/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)
                                                                                                                                                  |                                                                                                                                                   
TGCCAGCGAAAAGGCGGTGGCGGCACCCGAGCCC‑GTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGAGCT  >  NC_002947/2742772‑2743064

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