Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,861,953 T→G 6.3% V14V (GTA→GTC)  purM ← phosphoribosylformylglycinamidine cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,861,9530TG 6.3% 198.9 / 3.7 64V14V (GTA→GTC) purMphosphoribosylformylglycinamidine cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base T (27/33);  new base G (2/2);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

AGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCTGTGAAGGCGC  >  NC_002947/1861811‑1862098
                                                                                                                                                |                                                                                                                                                 
agaCCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCg                                                                                                                                             <  1:222049/149‑1 (MQ=255)
   ccAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCg                                                                                                                                          >  3:298970/1‑149 (MQ=255)
     aGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGTGTAGACATCGACGCCGGCgaa                                                                                                                                        <  7:302491/149‑3 (MQ=255)
           ggCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTa                                                                                                                                  <  3:210870/149‑1 (MQ=255)
              tgctTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCt                                                                                                                               >  8:210180/1‑149 (MQ=255)
                 ttGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCag                                                                                                                            >  7:139923/1‑149 (MQ=255)
                    tAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCaggga                                                                                                                         <  7:210180/149‑1 (MQ=255)
                    tAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCaggga                                                                                                                         <  8:131697/149‑1 (MQ=255)
                      gccggccgGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCagggag                                                                                                                        <  8:258353/148‑1 (MQ=255)
                              ggATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgcttg                                                                                                               >  5:263867/1‑149 (MQ=255)
                               gATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCGGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgcttgc                                                                                                              <  7:12139/149‑1 (MQ=255)
                               gATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgcttgc                                                                                                              <  6:263867/149‑1 (MQ=255)
                                   tcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGTGTAGACATCGACGCCGGCgaa                                                                                                                                        >  7:283427/1‑117 (MQ=255)
                                   tcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGTGTAGACATCGACGCCGGCgaa                                                                                                                                        <  8:283427/119‑3 (MQ=255)
                                   tcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCAt                                                                                                          >  1:243217/1‑149 (MQ=255)
                                    cGCAGAGGGCGCCGAAGCCGCGCAGTCCACCCATGACTTCAGGGCGTGCTGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACGAGTGCTTCGCCGGCGTCGATGTGTACACCGGCGTCCTTGTATCTCAGGGATGGTTGCTTGCTCAt                                                                                                          <  5:202277/148‑1 (MQ=255)
                                    cGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCAtt                                                                                                         >  4:45664/1‑149 (MQ=255)
                                     gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTc                                                                                                            <  4:197643/145‑1 (MQ=255)
                                     gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCa                                                                                                           <  2:61952/146‑1 (MQ=255)
                                          gggCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCa                                                                                                   <  4:286122/149‑1 (MQ=255)
                                              gcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGcc                                                                                               >  3:157335/1‑149 (MQ=255)
                                              gcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGcc                                                                                               >  6:117730/1‑149 (MQ=255)
                                               ccGAATCCGCCCAGTGCACCCAAGTCTTCAGGGGGTGCGGTGCGGTTTTTCCTAAGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGc                                                                                                <  1:268394/149‑1 (MQ=255)
                                               ccGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTc                                                                                                                                        <  6:255320/107‑1 (MQ=255)
                                                 gAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCttt                                                                                            <  4:100914/149‑1 (MQ=255)
                                                  aaGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCtt                                                                                             >  1:125894/1‑147 (MQ=255)
                                                   aGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCG‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGTTGTCTACACCGGCGTCCTTGTAGCTCATGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTag                                                                                          <  4:298970/149‑1 (MQ=255)
                                                    gccgccCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTagg                                                                                         <  7:94692/149‑1 (MQ=255)
                                                       gccCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTaggagg                                                                                      <  3:25234/149‑1 (MQ=255)
                                                               cACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGAtt                                                                              <  1:282604/149‑1 (MQ=255)
                                                                     tGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTaaa                                                                        <  2:270054/149‑1 (MQ=255)
                                                                         ttCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGGCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAAc                                                                     >  1:185634/1‑148 (MQ=255)
                                                                                  tgcggtgc‑gTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTGAAAAACGACCATGACg                                                           >  7:237555/1‑149 (MQ=255)
                                                                                        tgtgTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTc                                                                             <  2:48535/123‑1 (MQ=255)
                                                                                           gTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGc                                                   >  5:125165/1‑149 (MQ=255)
                                                                                           gTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGc                                                   >  7:47352/1‑149 (MQ=255)
                                                                                           gTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGc                                                   >  7:7816/1‑149 (MQ=255)
                                                                                            ttttGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGcc                                                  <  4:157335/149‑1 (MQ=255)
                                                                                                  c‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGc                                               <  8:194889/146‑1 (MQ=255)
                                                                                                  c‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGc                                               >  7:194889/1‑146 (MQ=255)
                                                                                                    aCGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGAt                                           <  3:45664/149‑1 (MQ=255)
                                                                                                    aCGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGACGGCTGa                                            >  6:186728/1‑148 (MQ=255)
                                                                                                        ccTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAg                                                                                                  <  4:108143/90‑1 (MQ=255)
                                                                                                        ccTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAgc                                       <  2:243217/149‑1 (MQ=255)
                                                                                                             aTGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAgcggcgg                                  <  8:139923/149‑1 (MQ=255)
                                                                                                             aTGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAgcggcgg                                  >  5:18090/1‑149 (MQ=255)
                                                                                                              tGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAgcggcggc                                 <  2:185634/149‑1 (MQ=255)
                                                                                                                 gTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGt                              >  4:221048/1‑149 (MQ=255)
                                                                                                                 gTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGt                              <  5:23527/149‑1 (MQ=255)
                                                                                                                  ttCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTc                             >  3:254457/1‑149 (MQ=255)
                                                                                                                     gACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGa                          <  5:295979/149‑1 (MQ=255)
                                                                                                                       ccAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGaaa                        >  2:213782/1‑149 (MQ=255)
                                                                                                                          gTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACg                      >  2:267466/1‑148 (MQ=255)
                                                                                                                          gTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGt                     >  6:267918/1‑149 (MQ=255)
                                                                                                                             cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAg                  >  3:173949/1‑149 (MQ=255)
                                                                                                                                cGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTg                <  5:267918/148‑1 (MQ=255)
                                                                                                                                 gCCGGTGTAGACATCGACGCCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGt              >  7:130917/1‑149 (MQ=255)
                                                                                                                                   cGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCt            >  3:44855/1‑149 (MQ=255)
                                                                                                                                    ggCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTTCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCtg           >  2:240650/1‑149 (MQ=255)
                                                                                                                                          gATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCTGTGAAgg     >  2:284112/1‑149 (MQ=255)
                                                                                                                                            tGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCTGTGAAGgcg   <  6:18090/149‑1 (MQ=255)
                                                                                                                                              tCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGaaa                        <  4:53466/126‑1 (MQ=255)
                                                                                                                                              tCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGaaa                        >  3:53466/1‑126 (MQ=255)
                                                                                                                                              tCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCTGTGAAGgcgc  <  8:47/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGC‑GTTTTGCC‑ACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGTCTGAAACGTCAGTGGTCTGTGAAGGCGC  >  NC_002947/1861811‑1862098

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