Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,223,854 T→C 5.2% R9R (CGT→CGC)  metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,8540TC 5.2% 257.9 / 3.8 77R9R (CGT→CGC) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base C (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

CGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTG  >  NC_000913/4223706‑4223995
                                                                                                                                                    |                                                                                                                                             
cGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgt                                                                                                                                               <  3:317635/149‑1 (MQ=255)
 gaggAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTCGAATTTTTATTAAATCGGTGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAGAGTGGAACAACtgcgtg                                                                                                                                              >  2:205152/4‑149 (MQ=255)
 gATGAGTTAATGTTGGACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                                                                                                                              >  2:166792/1‑149 (MQ=255)
 gATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                                                                                                                                              >  2:49449/1‑149 (MQ=255)
     aGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                                                                                                                                          <  4:362002/149‑1 (MQ=255)
            gTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt                                                                                                                                   >  4:535214/1‑149 (MQ=255)
             ttGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATg                                                                                                                                  >  4:7468/1‑149 (MQ=255)
               gAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa                                                                                                                                >  4:31357/1‑149 (MQ=255)
                aaCAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGAGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc                                                                                                                               >  2:34390/1‑149 (MQ=255)
                 aCAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                                                                                                                              >  4:160349/1‑149 (MQ=255)
                      tctcATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgt                                                                                                                                     >  1:106342/1‑135 (MQ=255)
                      tctcATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgt                                                                                                                                     <  2:106342/137‑3 (MQ=255)
                          aTGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt                                                                                                                     >  2:615362/1‑149 (MQ=255)
                           tGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTg                                                                                                                    <  3:73621/149‑1 (MQ=255)
                             ttGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCt                                                                                                                  >  3:104794/1‑149 (MQ=255)
                               gCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGTAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTgg                                                                                                                >  4:352184/1‑149 (MQ=255)
                                      ggtCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTAAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGAcggcgg                                                                                                         <  3:53991/149‑1 (MQ=255)
                                       gtCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGt                                                                                                        <  4:104794/149‑1 (MQ=255)
                                        tCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTa                                                                                                       >  4:529813/1‑149 (MQ=255)
                                           cTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa                                                                                                                                >  2:238875/1‑121 (MQ=255)
                                           cTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa                                                                                                                                <  1:238875/121‑1 (MQ=255)
                                           cTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATgg                                                                                                    <  2:526702/149‑1 (MQ=255)
                                            ttttACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGTAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATggg                                                                                                   >  1:392595/1‑149 (MQ=255)
                                                 ccACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCa                                                                                              >  3:339947/1‑149 (MQ=255)
                                                  cacaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCAt                                                                                             >  3:613949/1‑149 (MQ=255)
                                                  cacaGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCAt                                                                                             <  2:68620/149‑1 (MQ=255)
                                                    caGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGa                                                                                           >  2:253821/1‑149 (MQ=255)
                                                        tGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCa                                                                                       <  1:165334/149‑1 (MQ=255)
                                                          cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCaga                                                                                     <  1:34390/149‑1 (MQ=255)
                                                          cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCaga                                                                                     <  1:59850/149‑1 (MQ=255)
                                                          cGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTAAGCGTGTCCGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCaga                                                                                     <  3:535214/149‑1 (MQ=255)
                                                           gTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                                                                    >  4:497936/1‑149 (MQ=255)
                                                           gTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCagag                                                                                    >  4:100638/1‑149 (MQ=255)
                                                            tttATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGt                                                                                   >  4:489260/1‑149 (MQ=255)
                                                                  ccAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGa                                                                             >  3:29951/1‑149 (MQ=255)
                                                                   cAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGAc                                                                            <  4:289097/149‑1 (MQ=255)
                                                                     gTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTg                                                                          <  1:112836/149‑1 (MQ=255)
                                                                     gTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTg                                                                          >  1:575445/1‑149 (MQ=255)
                                                                       aTGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGaa                                                                        <  2:181323/149‑1 (MQ=255)
                                                                       aTGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTAGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGaa                                                                        <  4:351729/149‑1 (MQ=255)
                                                                         ggTTTGTTGAATTTTTATTAAATCTGGGGTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATAATCCAGAGTTATCGACTGAACg                                                                      <  1:379025/149‑1 (MQ=255)
                                                                            ttgttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAg                                                                   <  4:613949/149‑1 (MQ=255)
                                                                              gttgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGcc                                                                 <  1:615362/149‑1 (MQ=255)
                                                                                tgAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGa                                                               <  3:7468/149‑1 (MQ=255)
                                                                                 gAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAt                                                              >  4:164707/1‑149 (MQ=255)
                                                                                    tttttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGAtttt                                                           >  1:204204/1‑149 (MQ=255)
                                                                                     ttttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgt                                                                                                                                     <  4:551394/74‑3 (MQ=255)
                                                                                     ttttATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGCACGCAGttgt                                                                                                                                     >  3:551394/1‑72 (MQ=255)
                                                                                        tattaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtg                                                       <  3:489260/149‑1 (MQ=255)
                                                                                         attaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg                                                      >  3:332616/1‑149 (MQ=255)
                                                                                         attaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtgg                                                      <  1:166792/149‑1 (MQ=255)
                                                                                          ttaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAAc                                                                       <  4:535053/131‑1 (MQ=255)
                                                                                          ttaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCgtggt                                                     <  3:497936/149‑1 (MQ=255)
                                                                                          ttaAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGAGGTATGGGTACCATGATCCAGAGTTATCGACTGAAc                                                                       >  3:535053/1‑131 (MQ=255)
                                                                                             aaTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGaa                                                  <  2:575445/149‑1 (MQ=255)
                                                                                               tCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACg                                                >  4:487551/1‑149 (MQ=255)
                                                                                                    gTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTg                                           >  1:339645/1‑149 (MQ=255)
                                                                                                       gAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCg                                        <  3:160349/149‑1 (MQ=255)
                                                                                                           gtgtCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTg                                    <  4:29212/149‑1 (MQ=255)
                                                                                                                 ggAGTAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGGTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGCGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATg                              >  2:571067/1‑149 (MQ=255)
                                                                                                                  gAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGc                             <  3:31357/149‑1 (MQ=255)
                                                                                                                    gCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGa                           >  1:132918/1‑149 (MQ=255)
                                                                                                                           tgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCaaa                    <  2:108333/149‑1 (MQ=255)
                                                                                                                           tgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCaaa                    <  1:95148/149‑1 (MQ=255)
                                                                                                                           tgAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCaaa                    <  2:92058/149‑1 (MQ=255)
                                                                                                                            gAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAg                   <  3:487551/149‑1 (MQ=255)
                                                                                                                            cagcagcAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTGTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAg                   >  2:528984/2‑149 (MQ=255)
                                                                                                                              gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGc                 >  4:122990/1‑149 (MQ=255)
                                                                                                                              gcagcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGc                 <  3:529813/149‑1 (MQ=255)
                                                                                                                                 gcaAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaac              <  4:290876/149‑1 (MQ=255)
                                                                                                                                  caAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaaca             >  2:602975/1‑149 (MQ=255)
                                                                                                                                   aaaGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaacaa            <  3:164707/149‑1 (MQ=255)
                                                                                                                                   aaaGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaa               <  1:83687/146‑1 (MQ=255)
                                                                                                                                   aaaGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaa               >  2:83687/1‑146 (MQ=255)
                                                                                                                                   aaaGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAACGcaacaa            <  2:339645/149‑1 (MQ=255)
                                                                                                                                    aaGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCaacaac           <  1:253821/149‑1 (MQ=255)
                                                                                                                                           acaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCaga                                                                                     >  4:338253/1‑68 (MQ=255)
                                                                                                                                           acaacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCaga                                                                                     <  3:338253/68‑1 (MQ=255)
                                                                                                                                             aacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGAc                          <  1:402237/125‑1 (MQ=255)
                                                                                                                                             aacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGAc                          >  2:402237/1‑125 (MQ=255)
                                                                                                                                             aacTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACctgctg  >  3:506916/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTG  >  NC_000913/4223706‑4223995

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑