Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,931,506 T→C 5.8% V209A (GTC→GCC)  purT → phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,931,5060TC 5.8% 204.1 / 2.7 69V209A (GTC→GCC) purTphosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base C (2/2);  total (32/37)
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.01e-01

CTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGC  >  NC_000913/1931359‑1931646
                                                                                                                                                   |                                                                                                                                            
ctTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGt                                                                                                                                              >  7:132422/1‑148 (MQ=255)
 ttCCGGCAAGGGGCAGAGGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGt                                                                                                                                              >  2:42783/1‑147 (MQ=255)
 ttCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGt                                                                                                                                              >  4:293461/1‑147 (MQ=255)
     ggCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAgcgc                                                                                                                                        >  5:131820/1‑149 (MQ=255)
     ggCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAgcgc                                                                                                                                        <  5:165306/149‑1 (MQ=255)
       cAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCgg                                                                                                                                      >  2:126107/1‑149 (MQ=255)
        aaGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGt                                                                                                                                     <  4:46359/149‑1 (MQ=255)
            ggCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGAt                                                                                                                                 <  7:279691/149‑1 (MQ=255)
               aGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGc                                                                                                                              >  7:300331/1‑149 (MQ=255)
                  cGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTc                                                                                                                           <  6:131820/149‑1 (MQ=255)
                   gTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTc                                                                                                                           >  1:135011/1‑148 (MQ=255)
                   gTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTc                                                                                                                           >  4:149371/1‑148 (MQ=255)
                    tttATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCa                                                                                                                         <  8:142028/149‑1 (MQ=255)
                      tATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCAt                                                                                                                        <  4:32031/148‑1 (MQ=255)
                      tATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCAt                                                                                                                        <  3:291721/148‑1 (MQ=255)
                       aTTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCAttt                                                                                                                      >  1:167843/1‑149 (MQ=255)
                        ttCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCTGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCAttt                                                                                                                      <  4:156844/148‑1 (MQ=255)
                          cGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCt                                                                                                                    <  3:149371/148‑1 (MQ=255)
                            ttCTGCTGATGAACTTGCTCAGGCATGGAAGTACGCTAAGAAATGCTTTAGAGCCTTAGCGGGCCGCGTTATTGTTTATGGCGTCGTTTAGTTTGACTTCGAAATTTCCCTGCTAAGCGTCTGCGAGGTGGATGGATTCCATTTCtgtg                                                                                                                 <  6:37187/149‑1 (MQ=255)
                                gCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCAc                                                                                                              >  4:44878/1‑148 (MQ=255)
                                gCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCAc                                                                                                              >  4:288847/1‑148 (MQ=255)
                                   gagCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGt                                                                                                          <  6:138718/149‑1 (MQ=255)
                                              tCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGc                                                                                                >  6:85213/1‑148 (MQ=255)
                                              tCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGc                                                                                                >  6:86336/1‑148 (MQ=255)
                                              tCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGc                                                                                                >  4:304591/1‑148 (MQ=255)
                                              tCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGc                                                                                                >  2:120153/1‑148 (MQ=255)
                                              tCAGGCATCGAAGTACGCTCAGCAAGGCGGTCGCGAAGGAGCGGGCCGCGTAATTGTAGAAGGCGTCGTTAAGTTTGACTTCGAATTTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCACTAGGTCATCGc                                                                                                >  2:19827/1‑148 (MQ=255)
                                                aGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAg                                                                                             >  5:289653/1‑149 (MQ=255)
                                                    aTGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTACACCAGTAGGTCATCGCCAgg                                                                                            >  7:168066/1‑146 (MQ=255)
                                                     tGGAAGTACGCTCAGCAAGGCGGACGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAAACGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCAAgg                                                                                            >  5:152395/1‑145 (MQ=255)
                                                           tACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGa                                                                                  <  1:126107/149‑1 (MQ=255)
                                                              gCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTa                                                                               >  5:139865/1‑149 (MQ=255)
                                                              gCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACGTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGCAGATGGCGACTa                                                                               >  5:92445/1‑149 (MQ=255)
                                                                  aGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGt                                                                           <  8:328518/149‑1 (MQ=255)
                                                                    cAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGa                                                                         <  3:293461/149‑1 (MQ=255)
                                                                    cAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGa                                                                         <  6:289653/149‑1 (MQ=255)
                                                                     aaGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGaa                                                                        <  6:187796/149‑1 (MQ=255)
                                                                     aaGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGaa                                                                        >  2:245044/1‑149 (MQ=255)
                                                                        gCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATcc                                                                     <  8:132422/149‑1 (MQ=255)
                                                                         cGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCt                                                                    <  3:19986/149‑1 (MQ=255)
                                                                          ggTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGACATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGACGATGGCGACTACCGTGAATCCTg                                                                   <  1:42783/149‑1 (MQ=255)
                                                                          ggTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTg                                                                   <  3:288847/149‑1 (MQ=255)
                                                                          ggTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTg                                                                   <  1:180066/149‑1 (MQ=255)
                                                                          ggTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTg                                                                   <  3:44878/149‑1 (MQ=255)
                                                                                  cGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACcac                                                           <  3:52478/149‑1 (MQ=255)
                                                                                        gggCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCaaa                                                     <  3:304591/149‑1 (MQ=255)
                                                                                          gCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGc                                                                 >  4:271275/1‑135 (MQ=255)
                                                                                          gCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGc                                                                 <  3:271275/135‑1 (MQ=255)
                                                                                           ccGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCCGCGCGGTGGATGGCGTCAATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATgc                                                  <  2:167843/149‑2 (MQ=255)
                                                                                           ccGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGa                                                  <  5:86336/149‑1 (MQ=255)
                                                                                           ccGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGa                                                  <  2:212842/149‑1 (MQ=255)
                                                                                           ccGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGCCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAAt                                                    >  7:163023/1‑147 (MQ=255)
                                                                                           ccGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGCCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAAt                                                    >  7:158185/1‑147 (MQ=255)
                                                                                              cgTAATTGTTGAAGGCGTCGTTAAGTTTGACTTAGAAATTACCCTGCTAACCGTCAGCGCGCTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGCATATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAgac                                               >  4:102123/1‑147 (MQ=255)
                                                                                                      ttgAAGGCGTCGTTTAGTTTGACTTCGAAATTACCCTGATAACCGTCAGCGCGGATGATGTCGTCCATTTATGTGCACCAGTAGGTCAACGCCACGAAGATCGCGACGACCGTGAATCATGGCAACCACATCAAATCAGCCCGCTTGcc                                       >  6:78914/1‑149 (MQ=255)
                                                                                                      ttgAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGcc                                       >  3:46200/1‑149 (MQ=255)
                                                                                                      ttgAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGcc                                       >  6:28567/1‑149 (MQ=255)
                                                                                                             agtcgtTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTACACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAAc                                <  8:168066/148‑1 (MQ=255)
                                                                                                                 gtTAAGTATGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTgc                            >  1:190918/1‑149 (MQ=255)
                                                                                                                  ttAAGTTTGACTTCGAAATTAACCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACAAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTgcg                           >  4:325165/1‑149 (MQ=255)
                                                                                                                     aGTTTGACTTCGAAATTACCCTGCTACACGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCa                         <  2:296295/148‑1 (MQ=255)
                                                                                                                            cTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTg                 >  1:252930/1‑149 (MQ=255)
                                                                                                                             ttCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGc                >  5:256999/1‑149 (MQ=255)
                                                                                                                               cGATATTTCCCTTGTAACAGGGAGCGCTGTTGATGTCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGTTGGCGACTACCGTGAATCCTTGTAACCACATCATATGAGCCCGCTTGCCCTGGATGGTGGGCAGGATATTGca               <  2:135011/148‑2 (MQ=255)
                                                                                                                                    ttACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTaaa         <  5:28567/149‑1 (MQ=255)
                                                                                                                                      aCCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAg        <  2:252930/148‑1 (MQ=255)
                                                                                                                                       cccTGCTAACCGCCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAgt       <  8:158185/148‑1 (MQ=255)
                                                                                                                                       cccTGCTAACCGCCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAgt       <  8:163023/148‑1 (MQ=255)
                                                                                                                                         cTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAgtggt    >  1:2418/1‑149 (MQ=255)
                                                                                                                                         cTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGGGCGCAGGAGATTGCCCGTAAAgtggt    <  6:92445/149‑1 (MQ=255)
                                                                                                                                           gCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGc  >  1:182348/1‑149 (MQ=255)
                                                                                                                                            cTAACCGTCAGCGCGGGGGATGGCGTCCATTTCTGTGCCCCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAgtggtg   <  3:102123/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGC  >  NC_000913/1931359‑1931646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: