Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 639,841 C→G 6.4% G302R (GGC→CGC)  PP_0551 ← MORN domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947639,8410CG 6.4% 202.2 / 3.6 62G302R (GGC→CGC) PP_0551MORN domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base G (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

ACCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGCCATCAACGCT  >  NC_002947/639696‑639988
                                                                                                                                                    |                                                                                                                                                   
aCCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATc                                                                                                                                                  <  8:241592/149‑1 (MQ=255)
aCCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCAGCGAtg                                                                                                                                                  >  1:259604/1‑148 (MQ=255)
   ttGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGGGACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGc                                                                                                                                                >  1:111907/1‑149 (MQ=255)
          tGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAg                                                                                                                                         >  5:1984/1‑148 (MQ=255)
            ggCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCAGCGATGGCAGccactaca                                                                                                                                      <  2:259604/149‑7 (MQ=255)
             gCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAAc                                                                                                                                     >  3:133017/1‑149 (MQ=255)
               aGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTc                                                                                                                                   >  6:64018/1‑149 (MQ=255)
                gAATCAGCCGGCCATGGCCATGGTTGGCATCGTTGAGGAAACCGCCCTTGTACTTGGTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACATTCCGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTc                                                                                                                                   <  1:97565/149‑1 (MQ=255)
                   tCAGCCGCCCATGGTCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCAGCGATGGCAGcca                                                                                                                                           >  4:329657/1‑136 (MQ=255)
                   tCAGCCGCCCATGGTCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCAGCGATGGCAGcca                                                                                                                                           <  3:329657/137‑2 (MQ=255)
                     aGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTc                                                                                                                             >  4:45511/1‑149 (MQ=255)
                        cgccCATGGCCATTGTAGGCCTCTGTGAGGAAACAGCCGATGTACTTGCTGCCGTCGGCCAGTTTCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGcc                                                                                                                          <  7:131596/149‑1 (MQ=255)
                                tccatggTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCGCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGc                                                                                                                   <  3:309565/147‑1 (MQ=255)
                                 ccatggTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGa                                                                                                                 <  7:91795/149‑1 (MQ=255)
                                    tggTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGcc                                                                                                              >  2:192247/1‑149 (MQ=255)
                                           cATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGAAGTATGTTACCTTGAAGTGGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCtt                                                                                                        >  1:121645/1‑149 (MQ=255)
                                            aTCGTTGAGTAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGa                                                                                                      <  2:136908/149‑1 (MQ=255)
                                            aTCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTg                                                                                                       >  2:360726/1‑148 (MQ=255)
                                               gTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAAtt                                                                                                   >  4:246383/1‑149 (MQ=255)
                                                   aGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCAc                                                                                                <  3:45511/148‑1 (MQ=255)
                                                     gAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGc                                                                                                                           <  8:197028/119‑1 (MQ=255)
                                                     gAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGc                                                                                                                           >  7:197028/1‑119 (MQ=255)
                                                        aCCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGc                                                                                                                   >  1:86153/1‑124 (MQ=255)
                                                        aCCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGc                                                                                                                   <  2:86153/124‑1 (MQ=255)
                                                          cgccgATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCcac                                                                                        <  2:70303/149‑1 (MQ=255)
                                                          cgccgATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCcac                                                                                        <  5:348950/149‑1 (MQ=255)
                                                          cgccgATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCcac                                                                                        <  2:111907/149‑1 (MQ=255)
                                                          cgccgATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCcac                                                                                        <  1:251674/149‑1 (MQ=255)
                                                            ccgATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCacac                                                                                      <  2:116792/149‑1 (MQ=255)
                                                                    aCTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCa                                                                               >  6:361347/1‑149 (MQ=255)
                                                                      ttGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATc                                                                             >  4:46768/1‑149 (MQ=255)
                                                                        gctgcCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCgg                                                                           >  1:81079/1‑149 (MQ=255)
                                                                         ctgcCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACACACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCgg                                                                           >  5:306466/1‑148 (MQ=255)
                                                                          tgcCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCtt                                                                                                        <  7:192067/118‑1 (MQ=255)
                                                                          tgcCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCtt                                                                                                        >  8:192067/1‑118 (MQ=255)
                                                                                 ggCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGt                                                                  <  1:192247/149‑1 (MQ=255)
                                                                                    cAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTc                                                                                   >  7:207871/1‑129 (MQ=255)
                                                                                    cAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTc                                                                                   <  8:207871/129‑1 (MQ=255)
                                                                                    cAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGc                                                               >  2:209431/1‑149 (MQ=255)
                                                                                       ttGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGcggc                                                            <  3:188102/149‑1 (MQ=255)
                                                                                       ttGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGcggc                                                            >  6:205648/1‑149 (MQ=255)
                                                                                                 ccTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGcggcctcggcct                                                    <  3:46768/147‑1 (MQ=255)
                                                                                                 ccTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTc                                                  <  2:121645/149‑1 (MQ=255)
                                                                                                  cTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTc                                                  >  7:336574/1‑148 (MQ=255)
                                                                                                   ttGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTc                                                  >  4:262007/1‑147 (MQ=255)
                                                                                                        ccGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCgg                                           <  3:246383/149‑1 (MQ=255)
                                                                                                        ccGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCgg                                           <  6:1984/149‑1 (MQ=255)
                                                                                                        ccGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCgg                                           <  1:269461/149‑1 (MQ=255)
                                                                                                             aGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTc                                      >  1:368538/1‑149 (MQ=255)
                                                                                                             aGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTc                                      <  5:361347/149‑1 (MQ=255)
                                                                                                                 cGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCtt                                   <  4:224659/149‑1 (MQ=255)
                                                                                                                 cGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCtt                                   <  7:90261/149‑1 (MQ=255)
                                                                                                                         ggCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCCGCCTTCCAGGCGGTTGTCCTTGAATTCg                            <  1:209431/148‑1 (MQ=255)
                                                                                                                          gCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGcc                          >  6:354403/1‑149 (MQ=255)
                                                                                                                            gAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGa                        <  5:64018/149‑1 (MQ=255)
                                                                                                                                    ccTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGc                >  3:95065/1‑149 (MQ=255)
                                                                                                                                      ttGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCAGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCTTTTTCGTAGCGTCCTCGGCATTCCAGGTGGTTGTCCTTGAATTCGCCGATGTACTGCTc              >  5:272959/1‑149 (MQ=255)
                                                                                                                                        gTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTc              <  2:185143/147‑1 (MQ=255)
                                                                                                                                        gTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCg             <  3:262007/148‑1 (MQ=255)
                                                                                                                                        gTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGc            >  4:259396/1‑149 (MQ=255)
                                                                                                                                           gTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGCCAt         >  4:5767/1‑149 (MQ=255)
                                                                                                                                                 tGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGCCATCAACg    >  4:353873/1‑149 (MQ=255)
                                                                                                                                                   ccATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTGCTGGAAGTCACGGAGCCGCACGTCGCCATCGGCATTGTCGTAGCGGCCTCGGCCGGCCAGGCGGTTGTCCGTGAATTCGCCGAGGTCCTGCTCGCCTGCACCGCt  >  8:160000/1‑149 (MQ=255)
                                                                                                                                                   ccATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGCCATCAACg    <  6:306466/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACAGACGCCAGTCGGCGAACTCGCCCTTGTAGTGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGCCATCAACGCT  >  NC_002947/639696‑639988

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: