Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,586,530 C→T 7.1% K252K (AAG→AAA)  PP_4914 ← phytase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,586,5300CT 7.1% 210.9 / 4.0 70K252K (AAG→AAA) PP_4914phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (36/29);  new base T (3/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

CCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGG  >  NC_002947/5586385‑5586675
                                                                                                                                                 |                                                                                                                                                 
ccAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTg                                                                                                                                                >  2:183818/1‑149 (MQ=255)
   gTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCACCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGccc                                                                                                                                             >  7:254824/1‑149 (MQ=255)
     gCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCt                                                                                                                                            <  2:270600/148‑1 (MQ=255)
       ttCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCg                                                                                                                                          >  3:214694/1‑148 (MQ=255)
        tCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGaa                                                                                                                                        >  7:8449/1‑149 (MQ=255)
         cGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGa                                                                                                                                         <  6:316592/147‑1 (MQ=255)
            ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCtt                                                                                                                                                 <  5:290640/136‑1 (MQ=255)
            ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCtt                                                                                                                                                 >  6:290640/1‑136 (MQ=255)
            ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccag                                                                                                                                    >  4:298053/1‑149 (MQ=255)
                gccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                                  >  4:116931/1‑147 (MQ=255)
                 ccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggccag                                                                                                                               <  4:243502/149‑1 (MQ=255)
                    aaCGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTc                                                                                                                            <  3:116931/149‑1 (MQ=255)
                     aCGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGACATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCg                                                                                                                                          <  1:2331/134‑1 (MQ=255)
                     aCGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGACATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCg                                                                                                                                          >  2:2331/1‑134 (MQ=255)
                       ggCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGTGTAAACAGCTTGCCCTCGAACCAGGCCAGTTCAGc                                                                                                                         <  8:8449/149‑1 (MQ=255)
                         cTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGGCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGa                                                                                                                       <  7:340644/149‑1 (MQ=255)
                                 aGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTcgcgc                                                                                                               >  2:321852/1‑149 (MQ=255)
                                    accGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGAc                                                                                                            <  8:254824/148‑1 (MQ=255)
                                          aCTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCTGGCCAGGTCAGCGaa                                                                                                                      <  7:135638/133‑1 (MQ=255)
                                          aCTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTAGAGGGTAAACAGCTTGCCCTCGAACCTGGCCAGGTCAGCGaa                                                                                                                      >  8:135638/1‑133 (MQ=255)
                                            tCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCt                                                                                                    >  7:300622/1‑149 (MQ=255)
                                            tCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCt                                                                                                    >  7:46855/1‑149 (MQ=255)
                                                   gcgcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTg                                                                                             <  8:46855/149‑1 (MQ=255)
                                                   gcgcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTg                                                                                             <  1:183818/149‑1 (MQ=255)
                                                    cgcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTgg                                                                                            <  3:298053/149‑1 (MQ=255)
                                                       aTCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCa                                                                                         >  1:322715/1‑149 (MQ=255)
                                                                gAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                                <  1:66007/149‑1 (MQ=255)
                                                                      ccAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCt                                                                                      >  1:276968/1‑137 (MQ=255)
                                                                      ccAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCt                                                                                      <  2:276968/137‑1 (MQ=255)
                                                                       cAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAt                                                                         >  6:143198/1‑149 (MQ=255)
                                                                         gCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAttt                                                                       <  1:310175/149‑1 (MQ=255)
                                                                             cGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGAc                                                                   >  7:192777/1‑149 (MQ=255)
                                                                              gCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGAcc                                                                  >  7:151579/1‑149 (MQ=255)
                                                                                        gCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTc                                                          <  4:214694/147‑1 (MQ=255)
                                                                                        gCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGc                                                        >  6:103703/1‑149 (MQ=255)
                                                                                          cACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTGAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAttt                                                                       <  7:340383/132‑1 (MQ=255)
                                                                                          cACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAttt                                                                       >  8:340383/1‑132 (MQ=255)
                                                                                          cACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTc                                                      <  6:89981/149‑1 (MQ=255)
                                                                                               tCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagcag                                                 >  2:188267/1‑149 (MQ=255)
                                                                                               tCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagc                                                   <  5:143198/147‑1 (MQ=255)
                                                                                                 gATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCTCCCCCCGCCTCGCTCAgcagca                                               >  6:195565/1‑149 (MQ=255)
                                                                                                 gATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAgcagca                                               >  6:358616/1‑149 (MQ=255)
                                                                                                 gATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAgcagca                                               >  6:363788/1‑149 (MQ=255)
                                                                                                    atCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGAcc                                                                  >  8:162658/1‑127 (MQ=255)
                                                                                                    atCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGAcc                                                                  <  7:162658/127‑1 (MQ=255)
                                                                                                       aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCaa                                          <  1:321852/148‑1 (MQ=255)
                                                                                                       aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAAc                                         <  2:322715/149‑1 (MQ=255)
                                                                                                       aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAAc                                         >  8:121256/1‑149 (MQ=255)
                                                                                                             gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGcc                                   >  4:14544/1‑149 (MQ=255)
                                                                                                              cacaCCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGccc                                  <  3:49594/149‑1 (MQ=255)
                                                                                                               acacCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCg                                 >  5:179422/1‑149 (MQ=255)
                                                                                                               acacCCGGTACGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCg                                 >  5:85484/1‑149 (MQ=255)
                                                                                                                   ccGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAgg                             >  1:269792/1‑149 (MQ=255)
                                                                                                                      gTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTcc                          <  3:14544/149‑1 (MQ=255)
                                                                                                                            gTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCg                    <  5:103703/149‑1 (MQ=255)
                                                                                                                             ttACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGa                   >  1:215308/1‑149 (MQ=255)
                                                                                                                               tggTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGa                       <  7:243140/141‑1 (MQ=255)
                                                                                                                               tggTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGa                       >  8:243140/3‑143 (MQ=255)
                                                                                                                                cGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActg                >  6:14315/1‑149 (MQ=255)
                                                                                                                                 gTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActgc               >  2:235809/1‑149 (MQ=255)
                                                                                                                                  ttCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActgct              >  8:210399/1‑149 (MQ=255)
                                                                                                                                  ttCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCT‑GCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagca                                                  <  3:357987/112‑1 (MQ=38)
                                                                                                                                  ttCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCT‑GCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagca                                                  >  4:357987/1‑112 (MQ=38)
                                                                                                                                    cGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTgg            >  3:37614/1‑149 (MQ=255)
                                                                                                                                     gAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTgg            <  2:215308/148‑1 (MQ=255)
                                                                                                                                     gAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGgc           >  3:346172/1‑149 (MQ=255)
                                                                                                                                         gTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTc       <  5:14315/149‑1 (MQ=255)
                                                                                                                                         gTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACACGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCACCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTc       >  1:290025/1‑149 (MQ=255)
                                                                                                                                          tAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTc       >  5:38598/1‑148 (MQ=255)
                                                                                                                                              cAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATgg  <  8:300622/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGG  >  NC_002947/5586385‑5586675

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