Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,050,012 C→T 5.1% C219Y (TGT→TAT) ‡ 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,0120CT 5.1% 408.1 / 7.2 118C219Y (TGT→TAT) ‡gcvTaminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (57/55);  new base T (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.77e-01

Only 100 of 124 total aligned reads displayed.

TTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCC  >  NC_000913/3049864‑3050160
                                                                                                                                                    |                                                                                                                                                    
ttCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCcac                                                                                                                                                      >  3:1107163/1‑149 (MQ=255)
 tCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCcaca                                                                                                                                                     <  1:238963/149‑1 (MQ=255)
     gACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGc                                                                                                                                                 <  2:243836/149‑1 (MQ=255)
     gACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGc                                                                                                                                                 >  2:837871/1‑149 (MQ=255)
        cGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTa                                                                                                                                              >  4:251118/1‑149 (MQ=255)
         gATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTaa                                                                                                                                             >  1:940639/1‑149 (MQ=255)
         gATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTaa                                                                                                                                             >  4:362311/1‑149 (MQ=255)
         gATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTaa                                                                                                                                             >  4:1036207/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)
                  aCGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCtt                                                                                                                                    >  2:71517/1‑149 (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)
                                  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)
                                              cAGCCCATGTTGGCGGCTAGAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCtt                                                                                                                                               <  3:834130/110‑1 (MQ=255)
                                              cAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCtt                                                                                                                                               >  4:834130/1‑110 (MQ=255)
                                              cAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGc                                                                                                                      <  3:60676/135‑1 (MQ=255)
                                                  ccATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCtttt                                                                                                    <  1:675555/149‑1 (MQ=255)
                                                    aTGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTc                                                                                                  <  2:731604/149‑1 (MQ=255)
                                                       ttGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCAtt                                                                                               <  3:1039495/149‑1 (MQ=255)
                                                            ggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTg                                                                                              <  2:690400/145‑1 (MQ=255)
                                                            ggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTg                                                                                              >  1:690400/1‑145 (MQ=255)
                                                            ggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCa                                                                                          >  3:79072/1‑149 (MQ=255)
                                                              cTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAgc                                                                                        <  1:18541/149‑1 (MQ=255)
                                                               tAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAgcg                                                                                       >  3:518585/1‑149 (MQ=255)
                                                                  aGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCaa                                                                                    >  3:150127/1‑149 (MQ=255)
                                                                   ggAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAAt                                                                                   <  1:71517/149‑1 (MQ=255)
                                                                    gagagaTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAAtt                                                                                  >  1:740707/1‑149 (MQ=255)
                                                                    gagagaTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAACCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAAtt                                                                                  <  3:499640/149‑1 (MQ=255)
                                                                       agaTGGTTTCGTCCATCTGCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCATAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCa                                                                               >  1:333892/1‑149 (MQ=255)
                                                                       agaTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCa                                                                               >  1:179434/1‑149 (MQ=255)
                                                                        gaTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACAcc                                                                                                                                         <  1:1132288/90‑1 (MQ=255)
                                                                        gaTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACAcc                                                                                                                                         >  2:1132288/1‑90 (MQ=255)
                                                                         aTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATa                                                                             <  3:190213/149‑1 (MQ=255)
                                                                         aTGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATa                                                                             >  1:302805/1‑149 (MQ=255)
                                                                            gTTTCGTCCATCTCTTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGcc                                                                          >  3:85848/1‑149 (MQ=255)
                                                                               tCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGc                                                                       <  3:1036207/149‑1 (MQ=255)
                                                                                cGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAGGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCt                                                                      <  2:186344/149‑1 (MQ=255)
                                                                                     aTCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCAcc                                                                 >  4:401853/1‑149 (MQ=255)
                                                                                       ctcCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCgg                                                               >  3:871784/1‑149 (MQ=255)
                                                                                       ctcCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCgg                                                               <  2:940639/149‑1 (MQ=255)
                                                                                        tcCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGt                                                              >  2:513949/1‑149 (MQ=255)
                                                                                         ccTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGta                                                             >  4:153692/1‑149 (MQ=255)
                                                                                         ccTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGta                                                             <  4:803259/149‑1 (MQ=255)
                                                                                          cTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGtat                                                            >  1:926532/1‑149 (MQ=255)
                                                                                          cTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGtat                                                            >  4:425099/1‑149 (MQ=255)
                                                                                          cTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATAGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGtat                                                            >  1:838656/1‑149 (MQ=255)
                                                                                            gACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCAt                                                                                                <  3:479967/111‑1 (MQ=255)
                                                                                            gACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCAt                                                                                                >  4:479967/1‑111 (MQ=255)
                                                                                                aTAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACAccc                                                                                                                                        <  1:449206/67‑1 (MQ=255)
                                                                                                aTAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACAccc                                                                                                                                        >  2:449206/1‑67 (MQ=255)
                                                                                                aTAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAg                                                      >  3:948032/1‑149 (MQ=255)
                                                                                                aTAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAg                                                      >  2:316163/1‑149 (MQ=255)
                                                                                                  aaaGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtg                                                    >  3:308719/1‑149 (MQ=255)
                                                                                                  aaaGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtg                                                    >  1:777342/1‑149 (MQ=255)
                                                                                                     gATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtggtg                                                 <  1:1074767/149‑1 (MQ=255)
                                                                                                     gATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtggtg                                                 <  1:837871/149‑1 (MQ=255)
                                                                                                     gATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtggtg                                                 >  2:78082/1‑149 (MQ=255)
                                                                                                           tGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGccc                                                                         <  3:1104835/119‑1 (MQ=255)
                                                                                                           tGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGccc                                                                         >  4:1104835/1‑119 (MQ=255)
                                                                                                               cGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACa                                       <  2:179434/149‑1 (MQ=255)
                                                                                                                 cTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGa                                     <  4:518585/149‑1 (MQ=255)
                                                                                                                  ttCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGAt                                    <  3:1028226/149‑1 (MQ=255)
                                                                                                                  ttCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGAt                                    <  3:621583/149‑1 (MQ=255)
                                                                                                                   tCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATc                                   <  4:85848/149‑1 (MQ=255)
                                                                                                                   tCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATc                                   <  2:732973/149‑1 (MQ=255)
                                                                                                                      aGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGATGGCAATAAACAGATcgcc                                >  3:465359/1‑149 (MQ=255)
                                                                                                                          gCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACa                                       >  3:263560/1‑138 (MQ=255)
                                                                                                                          gCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACa                                       <  4:263560/138‑1 (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)
                                                                                                                                    aCGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGc                                                                       <  1:1121988/96‑1 (MQ=255)
                                                                                                                                      gcgcgcCCAAGCCACATGGCTTAACACCGGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAg                >  4:703960/1‑149 (MQ=255)
                                                                                                                                      gcgcgcCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAg                <  4:291470/149‑1 (MQ=255)
                                                                                                                                       cgcgcCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGa               >  4:503800/1‑149 (MQ=255)
                                                                                                                                        gcgcCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGaa              <  2:543793/149‑1 (MQ=255)
                                                                                                                                         gtgttaaGCCATGTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCt                           <  4:544255/130‑1 (MQ=255)
                                                                                                                                         gtgttaaGCCATGTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCt                           >  3:544255/6‑135 (MQ=255)
                                                                                                                                          gcCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACg            <  2:302805/149‑1 (MQ=255)
                                                                                                                                           cccAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACgg           <  4:150127/149‑1 (MQ=255)
                                                                                                                                            ccAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGt          <  4:1068868/149‑1 (MQ=255)
                                                                                                                                             cAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAACAACGGtt         <  2:620545/149‑1 (MQ=255)
                                                                                                                                                gCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAAc             <  2:81393/142‑1 (MQ=255)
                                                                                                                                                   acaTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATcc   >  1:737967/1‑149 (MQ=255)
                                                                                                                                                   acaTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATcc   <  1:208849/149‑1 (MQ=255)
                                                                                                                                                   acaTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATcc   <  3:503800/149‑1 (MQ=255)
                                                                                                                                                    caTGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATccc  >  3:277923/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCC  >  NC_000913/3049864‑3050160

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