Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,194,887 T→G 12.5% intergenic (+29/+24) PP_4574 → / ← PP_4575 hypothetical protein/urea amidolyase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,194,8870TG12.5% 134.4 / 10.0 61intergenic (+29/+24)PP_4574/PP_4575hypothetical protein/urea amidolyase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base G (4/4);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCA  >  NC_002947/5194749‑5195032
                                                                                                                                          |                                                                                                                                                 
ccTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGaa                                                                                                                                         <  5:293978/149‑1 (MQ=255)
ccTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGa                                                                                                                                          >  8:223637/1‑148 (MQ=255)
      aGTACCCCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGcc                                                                                                                                   <  7:335462/149‑1 (MQ=255)
         acacCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCcagc                                                                                                                                >  2:175104/1‑149 (MQ=255)
         acacCGCCTCGACCAATCTGGTGCGCGACCTGATCGAGGACATGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCGTTTCGCGGGCTCGGCCGCGAAAGGGCcag                                                                                                                                 <  6:11169/148‑1 (MQ=255)
           acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCccgcag                                                                                                                              >  4:353392/1‑149 (MQ=255)
           acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCcagcag                                                                                                                              >  8:294670/1‑149 (MQ=255)
           acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGAAAg                                                                                                                                       >  6:221278/1‑140 (MQ=255)
           acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGAAAg                                                                                                                                       <  5:221278/140‑1 (MQ=255)
            ccGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGt                                                                                                                             <  1:49147/149‑1 (MQ=255)
               ccTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCtt                                                                                                                          <  6:159059/149‑1 (MQ=255)
                   gACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACg                                                                                                                       <  7:223637/148‑1 (MQ=255)
                    aCCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGAc                                                                                                                     >  8:91030/1‑149 (MQ=255)
                       cATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCaa                                                                                                                  >  7:96518/1‑149 (MQ=255)
                                cgcgACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAg                                                                                                          <  3:77131/148‑1 (MQ=255)
                                cgcgACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAg                                                                                                          <  7:285859/148‑1 (MQ=255)
                                     ccTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGCCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCCGTCTTACGACCAACGCCCCCGGTAAAg                                                                                                    >  5:9230/1‑149 (MQ=255)
                                             aGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctgct                                                                                            >  1:2073/1‑149 (MQ=255)
                                                    gATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCa                                                                                     <  7:200471/149‑1 (MQ=255)
                                                     aTGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCCACGCCCCAGGTAAAGCTGCTGAACTCGCCAg                                                                                    >  3:105992/1‑149 (MQ=255)
                                                           ccGCTGACCCTGTAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCg                                                                              >  2:263239/1‑149 (MQ=255)
                                                             gCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTc                                                                            >  5:365372/1‑149 (MQ=255)
                                                                gACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTg                                                                         >  5:22/1‑149 (MQ=255)
                                                                       cAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCaccacc                                                                  >  5:90110/1‑149 (MQ=255)
                                                                       cAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTAGCGGGCACGCCCGCGAAAGGGCCAGCAGTCGTACGACCAACGCCCCAGGTAAAGCGGCTGCTATCGCCAGGCCTCGTCCAGCaccacc                                                                  >  2:340414/1‑149 (MQ=255)
                                                                         tgggTTGGAGGACGGGATGCGCTAGCTGTCCATGTGAGGTGGCCTGTACCGGCCCTTGCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGc                                                                <  6:9230/148‑1 (MQ=255)
                                                                            gTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCGAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCGCTCGCCAGGCCTCGTCCTGCACCACCGCCCt                                                             >  3:244987/1‑149 (MQ=255)
                                                                             ttGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCt                                                             >  1:52402/1‑148 (MQ=255)
                                                                                gaggaCGGGATGCGCTACCTGTCCAAGTGAGGTGGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGaa                                                          <  5:106048/148‑1 (MQ=255)
                                                                                 aggaCGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAAcc                                                        <  2:2073/149‑1 (MQ=255)
                                                                                 aggaCGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAAcc                                                        >  3:189134/1‑149 (MQ=255)
                                                                                  ggaCGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCg                                                       <  6:22/149‑1 (MQ=255)
                                                                                   gaCGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTTCCGGCCCATTAGCGGGGTCGCCCGCGAAACGGCCCGCAGTCGTACGACCCAACCCCCACGTTAATATGCGGCGCTCGCCAGGCCGCCGGCTGCACCACCCCCCAGACCCGc                                                      >  8:311478/1‑149 (MQ=255)
                                                                                            cgcTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTc                                             <  7:294670/149‑1 (MQ=255)
                                                                                                 ccTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGTAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGccc                                                              >  4:21902/1‑127 (MQ=255)
                                                                                                 ccTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCg                                         <  6:90110/148‑1 (MQ=255)
                                                                                                 ccTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGCGCCCGCGAAAGGGCCAGTAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGccc                                                              <  3:21902/127‑1 (MQ=255)
                                                                                                     tCCAAGTGAGGTTGCCTGTACCGGCCCTTTGGCGGGCTCGCCCGCGAAAg                                                                                                                                       >  1:31774/1‑50 (MQ=255)
                                                                                                     tCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAg                                                                                                                                       <  2:31774/50‑1 (MQ=255)
                                                                                                      ccAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGc                                   >  2:294561/1‑149 (MQ=255)
                                                                                                       cAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTc                                                                            >  7:312403/1‑107 (MQ=255)
                                                                                                       cAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTc                                                                            <  8:312403/107‑1 (MQ=255)
                                                                                                       cAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCa                                  >  4:22709/1‑149 (MQ=255)
                                                                                                       cAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCa                                  >  1:95859/1‑149 (MQ=255)
                                                                                                        aaGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCAc                                 >  3:352492/1‑149 (MQ=255)
                                                                                                          gTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTg                               >  5:224317/1‑149 (MQ=255)
                                                                                                           tGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTgc                              <  1:340414/149‑1 (MQ=255)
                                                                                                                       tACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgcc                  >  7:273166/1‑149 (MQ=255)
                                                                                                                         ccGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccag                >  4:129157/1‑149 (MQ=255)
                                                                                                                         ccGGCCCTTTCGCGGGCGCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccag                <  1:152838/149‑1 (MQ=255)
                                                                                                                          cGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGc                                       <  4:363947/125‑1 (MQ=255)
                                                                                                                          cGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGc                                       >  3:363947/1‑125 (MQ=255)
                                                                                                                          cGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagc               <  8:96518/149‑1 (MQ=255)
                                                                                                                            gCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcgc             >  6:237465/1‑149 (MQ=255)
                                                                                                                            gCCCTTTCGCGGGCTAGCCCGCGAAAGGGCCAGCAGTCTTACGACGAACTCCCCAGGTAACGCTGCTGCTCTCGCCTGGCCTGGTCCTGCACCACCTCCCTGAACCGCACCACCGTCCCCGGCATGCTCTGGGCCAGTTgcgccagcgc             <  4:244987/149‑1 (MQ=255)
                                                                                                                                tttCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgc                     >  6:87117/1‑137 (MQ=255)
                                                                                                                                tttCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgc                     <  5:87117/137‑1 (MQ=255)
                                                                                                                                tttCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc         >  4:65473/1‑149 (MQ=255)
                                                                                                                                tttCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc         >  6:209528/1‑149 (MQ=255)
                                                                                                                                    gcgGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAgcggcg     >  6:371187/1‑149 (MQ=255)
                                                                                                                                       ggCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCa  <  7:91030/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
CCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCA  >  NC_002947/5194749‑5195032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: