Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,665,249 T→G 5.1% I207L (ATC→CTC)  PP_1459 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,665,2490TG 5.1% 254.5 / 4.2 79I207L (ATC→CTC) PP_1459CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (34/41);  new base G (2/2);  total (36/43)
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.96e-01

AGGCGCTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAG  >  NC_002947/1665102‑1665391
                                                                                                                                                   |                                                                                                                                              
aGGCGCTGTGCAGCGCTTCCAGGGTCAGCTCAGCCTTGGGCAGCGGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATg                                                                                                                                               >  8:110234/1‑149 (MQ=255)
    gcTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGtt                                                                                                                                           <  8:146412/149‑1 (MQ=255)
     cTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTTGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGTCGGGTGTGGCGGCTTACGATCATCTGCTCGACTATTCGTTc                                                                                                                                          <  3:43962/149‑1 (MQ=255)
                ttCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGATTGGTGTTGAGGATCGCTTCTACCAGGTTGATGTCGTTGTGGTAGAGCGGCAGGCGTGTGTGGAGGCTGAAGATCAGATGCTCGACGATGCGTTCGATCGGCtcgt                                                                                                                               >  4:185182/1‑149 (MQ=255)
                 tCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCtcgtc                                                                                                                              >  5:229153/1‑149 (MQ=255)
                  ccAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGAt                                                                                                                                       >  2:2181/1‑139 (MQ=255)
                  ccAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGAt                                                                                                                                       <  1:2181/139‑1 (MQ=255)
                  ccAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGAt                                                                                                                                       >  4:204616/1‑139 (MQ=255)
                  ccAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGAt                                                                                                                                       <  3:204616/139‑1 (MQ=255)
                   cAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCtcgtc                                                                                                                              >  3:300239/1‑147 (MQ=255)
                    aGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAg                                                                                                                           <  6:229153/149‑1 (MQ=255)
                    aGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGACGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAg                                                                                                                           <  3:151707/149‑1 (MQ=255)
                        tCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGtt                                                                                                                        >  1:83448/1‑148 (MQ=255)
                          aGCTCGGGCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGAt                                                                                                                     <  3:185182/149‑1 (MQ=255)
                           gCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAgc                                                                                                                          <  1:323958/143‑2 (MQ=255)
                           gCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAgc                                                                                                                          >  2:323958/1‑142 (MQ=255)
                            cTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGc                                                                                                                   >  1:151282/1‑149 (MQ=255)
                                   ttGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCAc                                                                                                            >  7:78799/1‑149 (MQ=255)
                                     gggCAGCTGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACt                                                                                                           <  7:110234/148‑1 (MQ=255)
                                     gggCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCAc                                                                                                            <  3:66627/147‑1 (MQ=255)
                                       gcagcaGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCg                                                                                                        >  2:81985/1‑149 (MQ=255)
                                                  cTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCAc                                                                                             <  1:165317/149‑1 (MQ=255)
                                                    gATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCAc                                                                                             >  7:193742/1‑147 (MQ=255)
                                                           ttGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGa                                                                                                                      <  7:84244/115‑1 (MQ=255)
                                                           ttGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGa                                                                                                                      >  8:84244/1‑115 (MQ=255)
                                                             ggTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCg                                                                                   <  6:167254/148‑1 (MQ=255)
                                                              gtgtTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGtt                                                                                 <  8:193742/149‑1 (MQ=255)
                                                               tgtTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGtt                                                                                 >  3:210261/1‑148 (MQ=255)
                                                                 ttGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGAc                                                                              >  5:160941/1‑149 (MQ=255)
                                                                    aGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACgg                                                                            <  8:78799/148‑1 (MQ=255)
                                                                      gATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCgg                                                                                                                                    >  7:260104/1‑90 (MQ=255)
                                                                      gATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCgg                                                                                                                                    <  8:260104/90‑1 (MQ=255)
                                                                      gATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGt                                                                                  >  6:108496/1‑140 (MQ=255)
                                                                      gATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGt                                                                                  <  5:108496/140‑1 (MQ=255)
                                                                        tCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATc                                                                       <  5:351360/149‑1 (MQ=255)
                                                                        tCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATc                                                                       <  2:83448/149‑1 (MQ=255)
                                                                          gCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCtt                                                                     <  5:192428/149‑1 (MQ=255)
                                                                           cTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTg                                                                    >  4:44771/1‑149 (MQ=255)
                                                                            ttCGACCTGGTCGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGa                                                                               <  5:151697/137‑1 (MQ=255)
                                                                            ttCGACCTGGTCGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGa                                                                               >  6:151697/1‑137 (MQ=255)
                                                                              cGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTc                                                                  <  6:317928/148‑1 (MQ=255)
                                                                              cGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTc                                                                  <  6:328099/148‑1 (MQ=255)
                                                                                   tGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGa                                                            >  4:151618/1‑149 (MQ=255)
                                                                                     gTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAg                                                                                          >  4:55856/1‑117 (MQ=255)
                                                                                     gTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGa                                                            <  3:338144/147‑1 (MQ=255)
                                                                                      ttGATGTGGTTGGGGTAGACCGGCTGGCGGGTGTGGCGGCTGACGATGAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTg                                                         <  5:265365/149‑1 (MQ=255)
                                                                                          tGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAgg                                                     <  3:151618/149‑1 (MQ=255)
                                                                                          tGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAgg                                                     >  8:64658/1‑149 (MQ=255)
                                                                                            tCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGcc                                                   >  6:24526/1‑149 (MQ=255)
                                                                                              gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGAt                                                                                                                     <  6:229592/81‑1 (MQ=255)
                                                                                              gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGAt                                                                                                                     >  5:229592/1‑81 (MQ=255)
                                                                                              gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCtt                                                                     >  8:145946/1‑129 (MQ=255)
                                                                                              gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCtt                                                                     <  7:145946/129‑1 (MQ=255)
                                                                                              gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGt                                                 >  7:254703/1‑149 (MQ=255)
                                                                                                  tgGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGa                                              <  1:81985/148‑1 (MQ=255)
                                                                                                  tgGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCGCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGa                                              <  3:300198/148‑1 (MQ=255)
                                                                                                    gTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAg                                             >  7:196107/1‑147 (MQ=255)
                                                                                                       gACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCTTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGAc                                        <  2:255832/149‑1 (MQ=255)
                                                                                                         ccGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGAccc                                      >  7:138326/1‑149 (MQ=255)
                                                                                                                gCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAg                               >  4:150790/1‑149 (MQ=255)
                                                                                                                        ggcggcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCg                                 <  8:135158/139‑1 (MQ=255)
                                                                                                                        ggcggcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCg                                 >  7:135158/1‑139 (MQ=255)
                                                                                                                         gcggcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                              <  1:246252/109‑1 (MQ=255)
                                                                                                                         gcggcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                              >  2:246252/1‑109 (MQ=255)
                                                                                                                         gcggcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGt                      <  8:254703/149‑1 (MQ=255)
                                                                                                                            gcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCt                   <  7:135494/149‑1 (MQ=255)
                                                                                                                                aCGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGtt                >  6:215307/1‑148 (MQ=255)
                                                                                                                                aCGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGt                 <  6:160941/147‑1 (MQ=255)
                                                                                                                                  gATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCAt             <  7:64658/149‑1 (MQ=255)
                                                                                                                                  gATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCAt             >  3:186311/1‑149 (MQ=255)
                                                                                                                                    tCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATcc           <  4:186311/149‑1 (MQ=255)
                                                                                                                                          gctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCGTGGCTAGGCGCTTc                        <  8:172512/130‑1 (MQ=255)
                                                                                                                                          gctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCTAGGCGCTTc                        >  7:172512/1‑130 (MQ=255)
                                                                                                                                           ctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCg     <  3:44771/148‑1 (MQ=255)
                                                                                                                                            tCGACGAGCCGATCGAACGCATCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCg     >  5:55473/1‑147 (MQ=255)
                                                                                                                                             cgacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAg  <  5:24526/149‑1 (MQ=255)
                                                                                                                                              gacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGaa   <  5:215307/147‑1 (MQ=255)
                                                                                                                                                  aTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATc            <  5:18331/134‑1 (MQ=255)
                                                                                                                                                  aTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGCATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATc            >  6:18331/1‑134 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AGGCGCTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAG  >  NC_002947/1665102‑1665391

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: