Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,599,894 C→G 8.0% P1196P (CCG→CCC PP_4924 ← subtilase family serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,599,8940CG 8.0% 178.7 / 4.3 59P1196P (CCG→CCCPP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base G (2/3);  total (26/33)
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.31e-01

ATCATGCCAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATAAAACATGGCTCCCCGGTG  >  NC_002947/5599748‑5600042
                                                                                                                                                   |                                                                                                                                                      
aTCATGCCAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTCGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCgg                                                                                                                                                      <  4:13985/149‑1 (MQ=255)
 tCATGCCAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCggg                                                                                                                                                     >  3:205349/1‑149 (MQ=255)
  cATGCCAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGt                                                                                                                                                    >  5:19047/1‑149 (MQ=255)
    tGCCAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTg                                                                                                                                                  >  6:221595/1‑149 (MQ=255)
      ccAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACCTCCCGGTAACCCAACCCCGGGTTGGc                                                                                                                                                <  1:175338/149‑1 (MQ=255)
      ccAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGc                                                                                                                                                <  4:205349/149‑1 (MQ=255)
      ccAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGc                                                                                                                                                >  3:2772/1‑149 (MQ=255)
       cAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGCCAGTGCGAGGATCTGCTGACCATCCCGGTAACCCAACCCCGGGTTGGCt                                                                                                                                               >  8:238228/1‑149 (MQ=255)
       cAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAAGATCCCGGTAACCCAACCCCGGGTTGGCt                                                                                                                                               <  1:182343/149‑1 (MQ=255)
       cAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCt                                                                                                                                               <  4:92042/149‑1 (MQ=255)
          aGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGc                                                                                                                                            <  5:268026/149‑1 (MQ=255)
              aTAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACa                                                                                                                                        >  3:184446/1‑149 (MQ=255)
                       gctggcGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGcc                                                                                                                                <  2:222472/148‑1 (MQ=255)
                                 aTGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTga                                                                                                                     <  3:239494/149‑1 (MQ=255)
                                   gCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCGgggttgggttaccg                                                                                                                                          >  6:85255/1‑122 (MQ=255)
                                   gCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCGgggttgggttaccg                                                                                                                                          <  5:85255/126‑5 (MQ=255)
                                    ccACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTgaga                                                                                                                   <  6:42239/148‑1 (MQ=255)
                                       ccGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGAt                                                                                                               <  5:221595/149‑1 (MQ=255)
                                                aaaTTCCGCCCTGCGTTGTAAGCCAATACGGTTGATGACTCAACCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGATGCCTTGCAACATCAACGCAACCACGCATGAGAAGAAcggcgcgg                                                                                                       >  3:233490/2‑148 (MQ=255)
                                                      cgcCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGc                                                                                                 <  3:303425/148‑1 (MQ=255)
                                                         ccTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGt                                                                                             <  4:2772/149‑1 (MQ=255)
                                                          cTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTg                                                                                            <  6:265384/149‑1 (MQ=255)
                                                             cGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGcc                                                                                          <  6:19047/148‑1 (MQ=255)
                                                              gTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCCACCCCGGGTTGGCTTGCAACATCAACGCCACCACTCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCt                                                                                       >  2:200997/1‑149 (MQ=255)
                                                               ttGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCt                                                                                       >  6:213565/1‑148 (MQ=255)
                                                               ttGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCt                                                                                       >  2:277324/1‑148 (MQ=255)
                                                                tGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCt                                                                                       >  6:305797/1‑147 (MQ=255)
                                                                    aGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCg                                                                                <  6:206983/149‑1 (MQ=255)
                                                                       ccATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGc                                                                             <  6:273553/149‑1 (MQ=255)
                                                                          ttGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGt                                                                          <  5:46264/149‑1 (MQ=255)
                                                                          ttGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGt                                                                          <  5:118585/149‑1 (MQ=255)
                                                                            gggTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTa                                                                         <  6:23602/148‑1 (MQ=255)
                                                                             ggTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAg                                                                        <  1:258716/148‑1 (MQ=255)
                                                                               ttGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTgaga                                                                                                                   <  6:41014/105‑1 (MQ=255)
                                                                               ttGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTgaga                                                                                                                   >  5:41014/1‑105 (MQ=255)
                                                                                  aTGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTg                                                                  <  4:269239/149‑1 (MQ=255)
                                                                                             aCTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGc                                                       >  8:137374/1‑149 (MQ=255)
                                                                                              cTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAAACCCAACCCCGGGTTGGCTTGCACCATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGACGTGCTGTAGGCGCTGCCGACCGACGc                                                       >  1:67957/1‑149 (MQ=255)
                                                                                               tttGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCAc                                                     >  2:260978/1‑149 (MQ=255)
                                                                                               tttGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCAc                                                     >  2:287971/1‑149 (MQ=255)
                                                                                                           cAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCtt                                         <  1:326997/149‑1 (MQ=255)
                                                                                                           cAGTGCGAGGATATGCTGAACATCCCGGTTACCCTCCCGCGGGTTGGCTTGCTACATCAACGCCACGACGGCTGATACGTTGGGCGCGGCGAAGCTGGTGGTCTGTGTCGTGCTGTTGGCGCTGCCGACCGACGCACGGCGCTCCGCtt                                         <  4:233490/149‑1 (MQ=255)
                                                                                                            aGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCtt                                         >  4:106174/1‑148 (MQ=255)
                                                                                                              tGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTc                                        >  6:165421/1‑147 (MQ=255)
                                                                                                              tGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTc                                        >  7:128800/1‑147 (MQ=255)
                                                                                                              tGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCaa                                      <  5:305797/149‑1 (MQ=255)
                                                                                                                 gAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGc                                   <  5:165421/149‑1 (MQ=255)
                                                                                                                 gAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGc                                   >  3:206749/1‑149 (MQ=255)
                                                                                                                       ctgctgAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGc                             >  4:103319/1‑149 (MQ=255)
                                                                                                                          ctgAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCAc                          >  3:104697/1‑149 (MQ=255)
                                                                                                                               cATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGAt                     <  4:108247/149‑1 (MQ=255)
                                                                                                                                aTCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATa                    <  5:213565/149‑1 (MQ=255)
                                                                                                                                    cGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATaaaa                 >  8:52168/1‑148 (MQ=255)
                                                                                                                                    cGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATAAAAc                <  8:128800/149‑1 (MQ=255)
                                                                                                                                       tAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATAAAACATg             <  7:52168/149‑1 (MQ=255)
                                                                                                                                         aagccaacccGGGGTTGGGTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACcgcg                                                 <  2:191077/109‑1 (MQ=255)
                                                                                                                                         aagccaacccGGGGTTGGGTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACcgcg                                                 >  1:191077/4‑112 (MQ=255)
                                                                                                                                         aaccaaacccGGGGTTGGCTTGAAACAACAAGGCCACAAGGCCTTTTAGTTGCGGGGCGGGTATGCTGTTGCCATTTTCTCCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATAAAACa               <  1:260978/148‑1 (MQ=255)
                                                                                                                                            ccaacccCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATAAAACATGGCTc         >  3:26200/1‑149 (MQ=255)
                                                                                                                                                 cccGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATaaaa                 <  3:156671/136‑1 (MQ=255)
                                                                                                                                                 cccGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATaaaa                 >  4:156671/1‑136 (MQ=255)
                                                                                                                                                 cccGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATAAAACATGGCTCCCCgg    >  6:24344/1‑149 (MQ=255)
                                                                                                                                                   cGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATAAAACATGGCTCCCCGGTg  >  3:120584/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ATCATGCCAAAGCCATAGTCGTGGCTGGCGTGCATGCCACCGCCATTCCAATTGCGCCCTGCGTTGTAAGCCCATTGGGTTGATGCGTCATCCACTTTGCGCGCAGACAGTGCGAGGATCTGCTGAACATCCCGGTAACCCAACCCCGGGTTGGCTTGCAACATCAACGCCACCACGCCTGAGACGATCGGCGCGGCGAAGCTGGTGCCCTGAGTCGTGCTGTAGGCGCTGCCGACCGACGCACCGCGCTCCGCTTCAAGGCTGACGCCACTGGATAAAACATGGCTCCCCGGTG  >  NC_002947/5599748‑5600042

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: