Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,432,202 C→G 6.1% L241L (CTG→CTC)  PP_3927 ← type 12 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,432,2020CG 6.1% 211.9 / 4.1 65L241L (CTG→CTC) PP_3927type 12 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base G (2/2);  total (32/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 = 9.30e-01

CTACACCACCGGCACCACGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGTTGATCGGTTTCGACGTAGGT  >  NC_002947/4432054‑4432348
                                                                                                                                                    |                                                                                                                                                  
cTACACCACCGGCACCACGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGc                                                                                                                                                    <  3:37199/149‑1 (MQ=255)
      caccGGCACCACGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAgcgc                                                                                                                                               >  4:104241/1‑148 (MQ=255)
      caccGGCACCACGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGcc                                                                                                                                              >  8:297799/1‑149 (MQ=255)
             accacGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAgg                                                                                                                                       <  4:69963/149‑1 (MQ=255)
              ccacGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAggg                                                                                                                                      <  7:297799/149‑1 (MQ=255)
               cacGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAggg                                                                                                                                      >  7:215453/1‑148 (MQ=255)
               cacGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGACTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAggg                                                                                                                                      >  3:331173/1‑148 (MQ=255)
                  gCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGGAGGCGCTGCCAgggcag                                                                                                                                   <  8:92983/148‑4 (MQ=255)
                  gCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGGAGGCGCTGCCAgggcag                                                                                                                                   >  7:92983/1‑145 (MQ=255)
                  gCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATg                                                                                                                                  <  3:104241/149‑1 (MQ=255)
                        tGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGGAGGCGCTGCCAgggcaga                                                                                                                                  <  5:48094/143‑5 (MQ=255)
                        tGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGGAGGCGCTGCCAgggcaga                                                                                                                                  >  6:48094/1‑139 (MQ=255)
                        tGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTgg                                                                                                                            >  6:44385/1‑149 (MQ=255)
                        tGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTgg                                                                                                                            >  6:210054/1‑149 (MQ=255)
                        tGTCGACGGTGGCTTCCTGTCCTGATCAGACTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTgg                                                                                                                            <  4:331173/149‑1 (MQ=255)
                             aCGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTc                                                                                                                       >  8:73176/1‑149 (MQ=255)
                                    cTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAgc                                                                                                                <  8:196334/149‑1 (MQ=255)
                                    cTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAgc                                                                                                                >  1:262631/1‑149 (MQ=255)
                                           tcctgATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCa                                                                                                         <  5:210054/149‑1 (MQ=255)
                                                aTCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCt                                                                                                      <  5:44385/147‑1 (MQ=255)
                                                 tCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGTGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCa                                                                                                   <  1:77556/149‑1 (MQ=255)
                                                 tCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTg                                                                                                     >  7:170346/1‑147 (MQ=255)
                                                 tCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTg                                                                                                     >  6:152280/1‑147 (MQ=255)
                                                  cAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAg                                                                                                  >  8:26635/1‑149 (MQ=255)
                                                       tttGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGt                                                                                             >  7:172168/1‑149 (MQ=255)
                                                         tGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGt                                                                                             <  2:232599/147‑1 (MQ=255)
                                                          gACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGa                                                                                           <  7:26635/148‑1 (MQ=255)
                                                          gACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGa                                                                                           <  5:152280/148‑1 (MQ=255)
                                                          gACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAg                                                                                          >  4:153943/1‑149 (MQ=255)
                                                           aCGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGc                                                                                         >  5:162377/1‑149 (MQ=255)
                                                             ggCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGc                                                                                         >  2:7903/1‑147 (MQ=255)
                                                             ggCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCa                                                                                       <  7:73176/149‑1 (MQ=255)
                                                                gACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCGCGACAAGCGCTGTCAGCTGCAGGCCGTGAGCCAGc                                                                                     >  4:207523/1‑148 (MQ=255)
                                                                   cAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGt                                                                                                                           <  6:248180/107‑1 (MQ=255)
                                                                    aaGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGc                                                                                <  2:262631/149‑1 (MQ=255)
                                                                     aGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCg                                                                               >  4:14698/1‑149 (MQ=255)
                                                                      ggTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCgg                                                                              >  4:47632/1‑149 (MQ=255)
                                                                        tgtAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTg                                                                            <  3:47632/149‑1 (MQ=255)
                                                                              tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGAt                                                                                                                                   <  5:255040/88‑1 (MQ=255)
                                                                              tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGAt                                                                                                                                   >  6:255040/1‑88 (MQ=255)
                                                                              tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGt                                                                                             <  5:52107/126‑1 (MQ=255)
                                                                              tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGt                                                                                             >  6:52107/1‑126 (MQ=255)
                                                                              tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATc                                                                         <  5:352994/146‑1 (MQ=255)
                                                                              tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATc                                                                         >  6:352994/1‑146 (MQ=255)
                                                                              tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACt                                                                      >  4:85971/1‑149 (MQ=255)
                                                                                             cAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTg                                                        <  3:153943/148‑1 (MQ=255)
                                                                                                gctgctCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTc                                                    <  6:162377/149‑1 (MQ=255)
                                                                                                    ctCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTTGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCAc                                                 <  3:207523/148‑1 (MQ=255)
                                                                                                            ggCGATATTCGCCGTCTTCCCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGtt                                        >  5:22846/1‑149 (MQ=255)
                                                                                                              cGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTc                                      <  6:145931/149‑1 (MQ=255)
                                                                                                                atatTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCgg                                    >  3:301512/1‑149 (MQ=255)
                                                                                                                      gCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTgg                               <  3:85971/148‑1 (MQ=255)
                                                                                                                         gTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCa                           <  8:170346/149‑1 (MQ=255)
                                                                                                                            ttCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTc                                               <  4:267240/126‑1 (MQ=255)
                                                                                                                            ttCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTc                                               >  3:267240/1‑126 (MQ=255)
                                                                                                                              cGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCt                                                                                                                        >  6:275707/1‑51 (MQ=255)
                                                                                                                              cGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCt                                                                                                                        <  5:275707/51‑1 (MQ=255)
                                                                                                                               gCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGtt                     <  4:301512/149‑1 (MQ=255)
                                                                                                                               ccccTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGtt                     <  6:22846/148‑1 (MQ=255)
                                                                                                                                 ccTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGt                                             >  7:55305/1‑123 (MQ=255)
                                                                                                                                 ccTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGt                                             <  8:55305/123‑1 (MQ=255)
                                                                                                                                        cATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGTTGATCGGttt            <  3:14698/149‑1 (MQ=255)
                                                                                                                                        cATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGTTGATCGGgtt            >  4:228681/1‑149 (MQ=255)
                                                                                                                                              cccTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCCTTGATCGCTTTCGACg       >  3:288458/1‑148 (MQ=255)
                                                                                                                                               ccTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGTTGATCGGTTTCGACg       >  7:232540/1‑147 (MQ=255)
                                                                                                                                                cTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGTTGATCGGTTTCGACGTAg    <  1:7903/149‑1 (MQ=255)
                                                                                                                                                  accAGAGCATCCCCTGGGAGGCGCTGGTGATCGACCAGCGCTGTCAGCTGCAGCCCGTGAGCCCGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCAACTCGTGGGTTTCGGAGAGGGTCAGGCGTTGATCGGTTTCGACGTAGGt  >  8:320563/3‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTACACCACCGGCACCACGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGTTGATCGGTTTCGACGTAGGT  >  NC_002947/4432054‑4432348

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