Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,227,391 C→T 6.8% H191Y (CAC→TAC)  PP_1072 → leucine‑rich repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,227,3910CT 6.8% 170.5 / 2.3 59H191Y (CAC→TAC) PP_1072leucine‑rich repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (32/23);  new base T (2/2);  total (34/25)
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.74e-01

ATGAGCTTCCACGCGTTTGCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCGGTGGTGCCTGCGCTGGTGGGCACTG  >  NC_002947/1227259‑1227538
                                                                                                                                     |                                                                                                                                                   
aTGAGCTTCCACGCGTTTGCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGc                                                                                                                                     <  5:27051/149‑1 (MQ=255)
    gCTTCCACGCGTTTGCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGt                                                                                                                                 >  5:221402/1‑149 (MQ=255)
      ttCCACGCGTTTGCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGa                                                                                                                                        <  8:285548/140‑1 (MQ=255)
      ttCCACGCGTTTGCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGa                                                                                                                                        >  7:285548/1‑140 (MQ=255)
        ccACGCGTTTGCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTc                                                                                                                              >  4:142827/1‑148 (MQ=255)
          aCGCGTTTGCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAg                                                                                                                           >  3:56899/1‑149 (MQ=255)
           cgcgTTTGCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAgg                                                                                                                          >  6:316029/1‑149 (MQ=255)
                  gCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCAtt                                                                                                                    >  8:187027/1‑148 (MQ=255)
                         cATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTg                                                                                                            >  5:17588/1‑149 (MQ=255)
                             tCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCCTTTCGACGTGGCGa                                                                                                        >  6:144183/1‑149 (MQ=255)
                               gcgAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATg                                                                                                      <  7:72137/149‑1 (MQ=255)
                               gcgAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATg                                                                                                      <  7:165124/149‑1 (MQ=255)
                                    gTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGACTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGt                                                                                                 >  1:75364/1‑149 (MQ=255)
                                         cGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTc                                                                                                                  >  1:14096/1‑127 (MQ=255)
                                         cGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTc                                                                                                                  <  2:14096/127‑1 (MQ=255)
                                           cAGGGGGCAAGTACCAGGCGTTGCTCAGCCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGt                                                                                                 >  1:13989/1‑143 (MQ=255)
                                           cAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGt                                                                                                 <  2:13989/142‑1 (MQ=255)
                                                   aaGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACa                                                                                  <  3:221909/149‑1 (MQ=255)
                                                   aaGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTAGCTGCGTTGaaa                                                                                    >  1:15842/1‑147 (MQ=255)
                                                    aGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAg                                                                                 <  5:316029/149‑1 (MQ=255)
                                                    aGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAg                                                                                 <  7:187027/149‑1 (MQ=255)
                                                        ccAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACa                                                                             >  3:201706/1‑149 (MQ=255)
                                                          aGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGTTGGCCCTGAGGTTGAGCCGGGCCTGGCCGTTTCACGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACAtt                                                                           <  7:117195/149‑1 (MQ=255)
                                                            gCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGa                                                                         >  4:191153/1‑149 (MQ=255)
                                                                    ca‑cccaccATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCa                                                                 <  4:56899/149‑1 (MQ=255)
                                                                               cTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTAtttttt                                                       >  1:58806/1‑149 (MQ=255)
                                                                                 aCAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTAttttttta                                                     >  3:20969/1‑148 (MQ=255)
                                                                                 aCAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTAtttttt                                                       <  4:20969/147‑1 (MQ=255)
                                                                                  cAGCCTGAAGATGGCCGGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTAtttttt                                                       <  2:15842/146‑1 (MQ=255)
                                                                                     ccTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGc                                                                                           <  1:161732/107‑1 (MQ=255)
                                                                                     ccTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGc                                                                                           >  2:161732/1‑107 (MQ=255)
                                                                                     ccTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGtt                                                 <  6:17588/149‑1 (MQ=255)
                                                                                        gAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGc                                               >  6:226421/1‑148 (MQ=255)
                                                                                              ggCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTAt                                        >  6:244964/1‑149 (MQ=255)
                                                                                              ggCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCGGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTAt                                        >  7:118694/1‑149 (MQ=255)
                                                                                                cccGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTg                                      >  4:250529/1‑149 (MQ=255)
                                                                                                 ccGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGc                                     <  8:118694/149‑1 (MQ=255)
                                                                                                  cGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGcc                                    <  2:58806/149‑1 (MQ=255)
                                                                                                       gaTGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTa                                         <  8:141117/139‑1 (MQ=255)
                                                                                                       gaTGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTa                                         <  8:136896/139‑1 (MQ=255)
                                                                                                       gaTGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTa                                         >  7:141117/1‑139 (MQ=255)
                                                                                                       gaTGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTa                                         >  7:136896/1‑139 (MQ=255)
                                                                                                       gaTGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGAt                                 >  2:304802/1‑147 (MQ=255)
                                                                                                       gaTGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGAt                                 >  4:5973/1‑147 (MQ=255)
                                                                                                        aTGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGTGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAg                              <  3:191153/149‑1 (MQ=255)
                                                                                                        aTGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAg                              >  8:39031/1‑149 (MQ=255)
                                                                                                            gCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCgg                          >  1:296772/1‑149 (MQ=255)
                                                                                                            gCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCgg                          >  1:291989/1‑149 (MQ=255)
                                                                                                            gCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCgg                          <  4:201706/149‑1 (MQ=255)
                                                                                                                         gTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCGGTGGTGCCTGCGCt             <  5:226421/149‑1 (MQ=255)
                                                                                                                           ttCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCGGTGGTGCCTGCGCtgg           >  6:123747/1‑149 (MQ=255)
                                                                                                                            tCCCGCTTGTACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCggtgg                       <  2:184801/136‑1 (MQ=255)
                                                                                                                            tCCCGCTTGTACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCggtgg                       <  6:25786/136‑1 (MQ=255)
                                                                                                                            tCCCGCTTGTACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCggtgg                       >  1:184801/1‑136 (MQ=255)
                                                                                                                            tCCCGCTTGTACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCggtgg                       >  5:25786/1‑136 (MQ=255)
                                                                                                                               cGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCGGTGGTGCCTGCGCTGGTggg       >  4:110845/1‑149 (MQ=255)
                                                                                                                                gCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCGGTGGTGCCTGCGCTGGTGGGc      <  6:221402/149‑1 (MQ=255)
                                                                                                                                gCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCGGTGGTGCCTGCGCTGGTGGGc      >  2:211378/1‑149 (MQ=255)
                                                                                                                                    gCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCGGTGGTGCCTGCGCTGGTGGGCACtg  >  4:262151/1‑149 (MQ=255)
                                                                                                                                     |                                                                                                                                                   
ATGAGCTTCCACGCGTTTGCCTCCGCATGTCGCGAGGTCGACGCAGGGGGCAAGTACCAGGCGTTGCTCA‑CCCACCATCTACAGCCTGAAGATGGCCCGGAGGATGAGCCGGGCCAGGCCGTTTCCCGCTTGCACGGTCTTTTCGAGGCGAGTTTCAAGGCCCATTTCGAGGTGGCGATGCGGGTCGCTGCGTTGAAACAGGACATTGACGAGCGCACGTATTTTTTTCTGTTGCCACTTATTGCCGATAAGCCGGTGGTGCCTGCGCTGGTGGGCACTG  >  NC_002947/1227259‑1227538

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