Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,432,210 C→G 6.0% E239Q (GAG→CAG)  PP_3927 ← type 12 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,432,2100CG 6.0% 216.9 / 4.0 67E239Q (GAG→CAG) PP_3927type 12 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base G (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

ACCACGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGTTGATCGGTTTCGACGTAGGTCTGGTC  >  NC_002947/4432067‑4432354
                                                                                                                                               |                                                                                                                                                
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                                                                                                                                         >  6:255040/1‑88 (MQ=255)
                                                                 tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGAt                                                                                                                                         <  5:255040/88‑1 (MQ=255)
                                                                 tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGt                                                                                                   >  6:52107/1‑126 (MQ=255)
                                                                 tCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGt                                                                                                   <  5:52107/126‑1 (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                                                   <  8:55305/123‑1 (MQ=255)
                                                                                                                    ccTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGt                                                   >  7:55305/1‑123 (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)
                                                                                                                                           gcCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGTTGATCGGTTTCGACGTAGGTCTGGTc  >  4:53996/1‑149 (MQ=255)
                                                                                                                                             cTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGt                                                   <  3:206688/98‑1 (MQ=255)
                                                                                                                                             cTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGt                                                   >  4:206688/1‑98 (MQ=255)
                                                                                                                                               |                                                                                                                                                
ACCACGCTCAATGTCGACGGTGGCTTCCTGTCCTGATCAGGCTTTGACGGCGACCAAGGTGTAACTCAGGGGCAGGCGTGCAGGCTGCTCGCGCAGGCGATATTCGCCGTCTTCGCCCTTGACCATCTGCCCTGGCAGCGCCTCCCAGGGGATGCTCTGGTGCTCGACCAGCGCTGTCAGCTGCAGGCCGTGAGCCAGCAAGGCGGTGATCACTTCACCCAGGCCGTGGTTCCACTCGTGGGTTTCGGTGTGGGTCAGGCGTTGATCGGTTTCGACGTAGGTCTGGTC  >  NC_002947/4432067‑4432354

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