Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,050,013 A→G 5.1% C219R (TGT→CGT) ‡ gcvT ← aminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,050,0130AG 5.1% 395.6 / 8.0 118C219R (TGT→CGT) ‡gcvTaminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base A (57/55);  new base G (3/3);  total (60/58)
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

Only 100 of 127 total aligned reads displayed.

TCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCC  >  NC_000913/3049865‑3050160
                                                                                                                                                    |                                                                                                                                                   
tCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCcaca                                                                                                                                                     <  1:238963/149‑1 (MQ=255)
    gACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGc                                                                                                                                                 <  2:243836/149‑1 (MQ=255)
    gACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGc                                                                                                                                                 >  2:837871/1‑149 (MQ=255)
      ccGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCtt                                                                                                                                               <  2:684792/149‑1 (MQ=255)
       cGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTa                                                                                                                                              >  4:251118/1‑149 (MQ=255)
        gATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTaa                                                                                                                                             >  4:1036207/1‑149 (MQ=255)
        gATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTaa                                                                                                                                             >  1:940639/1‑149 (MQ=255)
        gATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTaa                                                                                                                                             >  4:362311/1‑149 (MQ=255)
         aTAAAGTCACGATCTGCCGGTTCCCAGGAGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCATGTGGCttgg                                                                                                                                             >  3:785048/1‑146 (MQ=255)
         aTAAAGTCACGATCTGCCGGTTCCCAGGAGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCATGTGGCttcg                                                                                                                                             <  4:785048/148‑3 (MQ=255)
          tAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAaca                                                                                                                                           <  2:765688/149‑1 (MQ=255)
           aaaGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGTCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAacac                                                                                                                                          <  1:1107258/149‑1 (MQ=255)
               tCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAACCACATGGCTTAACACCCGc                                                                                                                                      >  4:499640/1‑149 (MQ=255)
                 aCGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCtt                                                                                                                                    >  2:71517/1‑149 (MQ=255)
                   gATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTcc                                                                                                                                  <  4:16019/149‑1 (MQ=255)
                   gATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTcc                                                                                                                                  >  4:529815/1‑149 (MQ=255)
                    aTCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTcca                                                                                                                                 >  1:620545/1‑149 (MQ=255)
                      cTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTccacc                                                                                                                               >  3:1068868/1‑149 (MQ=255)
                                ccAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGcc                                                                                                                     >  3:291470/1‑149 (MQ=255)
                                 cAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTcca                                                                                                                                 >  3:1050776/1‑136 (MQ=255)
                                 cAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCa                                                                                                                    >  4:1028226/1‑149 (MQ=255)
                                      gATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAAt                                                                                                               >  4:621583/1‑149 (MQ=255)
                                        tGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCg                                                                                                             <  4:826483/149‑1 (MQ=255)
                                         ggTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCgg                                                                                                            >  4:190213/1‑149 (MQ=255)
                                          gTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGc                                                                                                           <  1:1051392/149‑1 (MQ=255)
                                          gTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGc                                                                                                           <  1:302243/149‑1 (MQ=255)
                                           tCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGcc                                                                                                          <  4:1107163/149‑1 (MQ=255)
                                           gccAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCATGTGGCttgg                                                                                                                                             >  4:740255/2‑112 (MQ=255)
                                            ccAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCg                                                                                                         >  1:732973/1‑149 (MQ=255)
                                             cAGCCCATGTTGGCGGCTAGAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCtt                                                                                                                                               <  3:834130/110‑1 (MQ=255)
                                             cAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCtt                                                                                                                                               >  4:834130/1‑110 (MQ=255)
                                             cAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGc                                                                                                                      <  3:60676/135‑1 (MQ=255)
                                             cAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGc                                                                                                                      >  4:60676/1‑135 (MQ=255)
                                                      ttGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCAtt                                                                                               <  3:1039495/149‑1 (MQ=255)
                                                           ggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTg                                                                                              >  1:690400/1‑145 (MQ=255)
                                                           ggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTg                                                                                              <  2:690400/145‑1 (MQ=255)
                                                           ggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCa                                                                                          >  3:79072/1‑149 (MQ=255)
                                                             cTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAgc                                                                                        <  1:18541/149‑1 (MQ=255)
                                                                 aGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCaa                                                                                    >  3:150127/1‑149 (MQ=255)
                                                                  ggAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAAt                                                                                   <  1:71517/149‑1 (MQ=255)
                                                                   gagagaTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAACCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAAtt                                                                                  <  3:499640/149‑1 (MQ=255)
                                                                      agaTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCa                                                                               >  1:179434/1‑149 (MQ=255)
                                                                       gaTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACAcc                                                                                                                                         >  2:1132288/1‑90 (MQ=255)
                                                                       gaTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACAcc                                                                                                                                         <  1:1132288/90‑1 (MQ=255)
                                                                        aTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATa                                                                             <  3:190213/149‑1 (MQ=255)
                                                                        aTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATa                                                                             >  1:302805/1‑149 (MQ=255)
                                                                           gTTTCGTCCATCTCTTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGcc                                                                          >  3:85848/1‑149 (MQ=255)
                                                                                    aTCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCAcc                                                                 >  4:401853/1‑149 (MQ=255)
                                                                                     tctcCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCg                                                                >  2:208849/1‑149 (MQ=255)
                                                                                      ctcCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCgg                                                               <  2:940639/149‑1 (MQ=255)
                                                                                      ctcCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCgg                                                               >  3:871784/1‑149 (MQ=255)
                                                                                       tcCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGt                                                              >  2:513949/1‑149 (MQ=255)
                                                                                        ccTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGta                                                             >  4:153692/1‑149 (MQ=255)
                                                                                         cTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGtat                                                            >  1:926532/1‑149 (MQ=255)
                                                                                           gACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCAt                                                                                                >  4:479967/1‑111 (MQ=255)
                                                                                            aCCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAAc                                                         <  3:362311/149‑1 (MQ=255)
                                                                                               aTAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACAccc                                                                                                                                        <  1:449206/67‑1 (MQ=255)
                                                                                               aTAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACAccc                                                                                                                                        >  2:449206/1‑67 (MQ=255)
                                                                                               aTAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAg                                                      >  3:948032/1‑149 (MQ=255)
                                                                                                 aaaGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtg                                                    >  1:777342/1‑149 (MQ=255)
                                                                                                 aaaGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtg                                                    >  3:308719/1‑149 (MQ=255)
                                                                                                    gATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtggtg                                                 <  1:1074767/149‑1 (MQ=255)
                                                                                                    gATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtggtg                                                 <  1:837871/149‑1 (MQ=255)
                                                                                                    gATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtggtg                                                 <  2:1059953/149‑1 (MQ=255)
                                                                                                    gATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtggtg                                                 >  2:78082/1‑149 (MQ=255)
                                                                                                     aTTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGtggtgg                                                >  1:263358/1‑149 (MQ=255)
                                                                                                     aTTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAAAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGtggtgg                                                <  4:79072/149‑1 (MQ=255)
                                                                                                        cATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCaa                                             <  2:740707/149‑1 (MQ=255)
                                                                                                          tGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGccc                                                                         >  4:1104835/1‑119 (MQ=255)
                                                                                                              ctcTTCCAGACGCAGTGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGAAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACa                                       <  2:333892/147‑1 (MQ=255)
                                                                                                              cGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACa                                       <  2:179434/149‑1 (MQ=255)
                                                                                                                 ttCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGAt                                    <  3:621583/149‑1 (MQ=255)
                                                                                                                 ttCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGAt                                    <  3:1028226/149‑1 (MQ=255)
                                                                                                                  tCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATc                                   <  4:85848/149‑1 (MQ=255)
                                                                                                                  tCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATc                                   <  2:732973/149‑1 (MQ=255)
                                                                                                                    cAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGACAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATcga                                 >  1:543793/1‑148 (MQ=255)
                                                                                                                     aGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGATGGCAATAAACAGATcgcc                                >  3:465359/1‑149 (MQ=255)
                                                                                                                         gCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACa                                       >  3:263560/1‑138 (MQ=255)
                                                                                                                         gCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGAACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATcgcccgcc                            >  1:741289/1‑149 (MQ=255)
                                                                                                                            gCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGc                         >  2:273313/1‑149 (MQ=255)
                                                                                                                                 tCACGCGCGCCCAAGCCACATGGCGTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTTGTGGCAATAAACAGATCGCCCGCCTGCACGcc                    >  2:126875/1‑149 (MQ=255)
                                                                                                                                   aCGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGc                                                                       >  2:1121988/1‑96 (MQ=255)
                                                                                                                                    cgcgcgCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCaaa                 <  4:308719/149‑1 (MQ=255)
                                                                                                                                     gcgcgcCCAAGCCACATGGCTTAACACCGGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAg                >  4:703960/1‑149 (MQ=255)
                                                                                                                                     gcgcgcCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAg                <  4:291470/149‑1 (MQ=255)
                                                                                                                                     gcgcgcCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAg                >  1:165373/1‑149 (MQ=255)
                                                                                                                                      cgcgcCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGa               >  4:503800/1‑149 (MQ=255)
                                                                                                                                       gcgcCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGaa              <  2:543793/149‑1 (MQ=255)
                                                                                                                                        gtgttaaGCCATGTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCt                           >  3:544255/6‑135 (MQ=255)
                                                                                                                                        gtgttaaGCCATGTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCt                           <  4:544255/130‑1 (MQ=255)
                                                                                                                                         gcCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACg            <  2:302805/149‑1 (MQ=255)
                                                                                                                                           ccAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGt          <  4:1068868/149‑1 (MQ=255)
                                                                                                                                            cAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGtt         <  3:506817/149‑1 (MQ=255)
                                                                                                                                              aGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTc       <  3:251118/149‑1 (MQ=255)
                                                                                                                                               gCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCTGCCTGCACGCCAAAGAACGGTTTCa      <  2:838656/149‑1 (MQ=255)
                                                                                                                                               gCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAAc             >  1:81393/1‑142 (MQ=255)
                                                                                                                                                  acaTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATcc   >  1:737967/1‑149 (MQ=255)
                                                                                                                                                  acaTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATcc   <  3:503800/149‑1 (MQ=255)
                                                                                                                                                    aTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgt                                                     <  1:1094997/97‑1 (MQ=255)
                                                                                                                                                    aTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAAc             >  3:639014/1‑137 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCC  >  NC_000913/3049865‑3050160

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