Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,354 Δ1 bp 6.7% coding (121/639 nt) purN → phosphoribosylglycinamide formyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,622,3540T. 6.7% 673.4 / 38.3 193coding (121/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (88/92);  new base . (6/7);  total (94/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 195 total aligned reads displayed.

TTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCG  >  NC_000913/2622207‑2622502
                                                                                                                                                    |                                                                                                                                                    
 tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTc                                                                                                                                                    >  3:91551/1‑149 (MQ=255)
   cGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAATCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCgg                                                                                                                                                  >  3:1226953/1‑149 (MQ=255)
     aacaacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGcc                                                                                                                                                <  4:651871/149‑1 (MQ=255)
       caacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCtt                                                                                                                                              <  3:1107569/149‑1 (MQ=255)
       caacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCtt                                                                                                                                              <  1:800469/149‑1 (MQ=255)
        aacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCtta                                                                                                                                             <  3:275647/149‑2 (MQ=255)
         acGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttat                                                                                                                                            <  2:1349987/149‑3 (MQ=255)
         acGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttat                                                                                                                                            <  4:582957/149‑3 (MQ=255)
          cgcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGaa                                                                                                                                           <  2:1233726/149‑1 (MQ=255)
          cgcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGaa                                                                                                                                           <  2:1104708/149‑1 (MQ=255)
               ggTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGACGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACgcgc                                                                                                                                      <  3:1453024/149‑1 (MQ=255)
                gTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGcgcc                                                                                                                                     <  1:116040/149‑1 (MQ=255)
                   aTCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGcgcccgc                                                                                                                                  >  2:1497379/1‑149 (MQ=255)
                        ataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGc                                                                                                                             <  3:1001508/149‑1 (MQ=255)
                        ataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGc                                                                                                                             >  3:222184/1‑149 (MQ=255)
                          aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                           <  2:1046181/149‑1 (MQ=255)
                           atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCggg                                                                                                                          <  2:153901/149‑1 (MQ=255)
                           atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCAGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCggg                                                                                                                          >  3:975462/1‑149 (MQ=255)
                            tGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGt                                                                                                                         >  2:658012/1‑149 (MQ=255)
                             gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACTCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                        <  3:877654/149‑1 (MQ=255)
                             gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                        >  1:1363177/1‑149 (MQ=255)
                              aaTATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTAt                                                                                                                       <  2:957148/149‑1 (MQ=255)
                                tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttat                                                                                                                                            <  1:309204/126‑3 (MQ=255)
                                tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttat                                                                                                                                            >  2:309204/1‑124 (MQ=255)
                                tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                           >  1:1641138/1‑143 (MQ=255)
                                tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                           <  2:1641138/143‑1 (MQ=255)
                                tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTg                                                                                                                     >  3:435309/1‑149 (MQ=255)
                                 atTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                                    >  2:1528444/1‑149 (MQ=255)
                                        tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCAt                                                                                                             >  3:577716/1‑149 (MQ=255)
                                        tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCAt                                                                                                             >  4:704770/1‑149 (MQ=255)
                                           ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCt                                                                                                                                               <  4:454175/112‑1 (MQ=255)
                                               ttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAAACAACAAAACTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATa                                                                                                            <  3:1255242/143‑1 (MQ=255)
                                                   ggCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGc                                                                                                  >  1:202092/1‑149 (MQ=255)
                                                    gCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGcc                                                                                                 >  2:369709/1‑149 (MQ=255)
                                                     cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCa                                                                                                <  1:618801/149‑1 (MQ=255)
                                                     cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCa                                                                                                <  4:381116/149‑1 (MQ=255)
                                                     cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCa                                                                                                <  4:1251353/149‑1 (MQ=255)
                                                      aaCGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAAACAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAg                                                                                               >  3:679860/1‑149 (MQ=255)
                                                           aaGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGt                                                                                          <  4:1288651/149‑1 (MQ=255)
                                                           aaGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGt                                                                                          >  4:1135759/1‑149 (MQ=255)
                                                               aaTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                      <  1:1208895/149‑1 (MQ=255)
                                                                aTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGttt                                                                                        >  3:361747/1‑146 (MQ=255)
                                                                   tACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGt                                                                                  >  4:1572894/1‑149 (MQ=255)
                                                                      aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                               <  2:757726/149‑1 (MQ=255)
                                                                      aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                               >  4:266574/1‑149 (MQ=255)
                                                                        gCAATTATTGACGCCTGTAAACCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttatt                                                                                                                                           >  2:58568/1‑84 (MQ=38)
                                                                        gCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                      <  4:958947/140‑1 (MQ=255)
                                                                         cAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                            >  2:1613086/1‑149 (MQ=255)
                                                                                  aCGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAc                                                                   <  1:901803/149‑1 (MQ=255)
                                                                                  aCGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAc                                                                   <  3:1572894/149‑1 (MQ=255)
                                                                                     ccTGTAAAACCAACAAGATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgg                                                                <  4:552746/149‑1 (MQ=255)
                                                                                     ccTGAAAAACCAACAAATTTAAAGGCACATTACGGGCAGTTTACAGTAATAAGGCAGACGCGTTCGGCCTTGCACGCTGCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACACTCGTGAAGCCTATGCCCgg                                                                <  2:532650/149‑1 (MQ=255)
                                                                                      cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCtt                                                                                                                                              >  1:494494/1‑70 (MQ=255)
                                                                                         tAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTc                                                                                 <  2:96564/128‑1 (MQ=255)
                                                                                         tAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTc                                                                                 >  1:96564/1‑128 (MQ=255)
                                                                                               caacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTc                                                      <  1:385323/149‑1 (MQ=255)
                                                                                                aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                                     >  1:1031823/1‑149 (MQ=255)
                                                                                                aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                                     >  1:855492/1‑149 (MQ=255)
                                                                                                aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                                     <  1:658012/149‑1 (MQ=255)
                                                                                                 acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGt                                                            <  1:869250/141‑1 (MQ=255)
                                                                                                 acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGt                                                            >  2:869250/1‑141 (MQ=255)
                                                                                                 acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                                    >  4:1325388/1‑149 (MQ=255)
                                                                                                 acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                                    <  2:784633/149‑1 (MQ=255)
                                                                                                  caaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATg                                                   <  3:384493/149‑1 (MQ=255)
                                                                                                   aaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGa                                                  <  3:500762/149‑1 (MQ=255)
                                                                                                    aaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                               >  3:1160664/1‑119 (MQ=255)
                                                                                                     aaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTc                                                                                 >  1:44129/1‑116 (MQ=255)
                                                                                                     aaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTc                                                                                 <  2:44129/116‑1 (MQ=255)
                                                                                                       ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                              >  4:948366/1‑149 (MQ=255)
                                                                                                       ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                              <  4:1640225/149‑1 (MQ=255)
                                                                                                          aaGGCGCCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGACCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAGTCATGAAATCGat                                           <  4:1090867/149‑2 (MQ=255)
                                                                                                            ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGtt                                                           <  3:1432830/131‑1 (MQ=255)
                                                                                                            ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGtt                                                           >  4:1432830/1‑131 (MQ=255)
                                                                                                            ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACAt                                         >  2:664237/1‑149 (MQ=255)
                                                                                                                         cAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAt                            <  1:59160/149‑1 (MQ=255)
                                                                                                                         cAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAt                            <  3:1629710/149‑1 (MQ=255)
                                                                                                                            ttttCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAtgtg                         <  2:855492/149‑1 (MQ=255)
                                                                                                                               tCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTc                      <  1:1613086/149‑1 (MQ=255)
                                                                                                                                 aGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAtt                                                       <  4:870714/114‑1 (MQ=255)
                                                                                                                                 aGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGt                    >  4:162368/1‑149 (MQ=255)
                                                                                                                                   cAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgc                  <  4:1083533/149‑1 (MQ=255)
                                                                                                                                   cAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgc                  >  3:1646497/1‑149 (MQ=255)
                                                                                                                                    aaTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgct                 >  1:1011076/1‑149 (MQ=255)
                                                                                                                                    aaTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgct                 >  4:409931/1‑149 (MQ=255)
                                                                                                                                    aaTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgct                 >  2:247558/1‑149 (MQ=255)
                                                                                                                                     aTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg                >  4:147878/1‑149 (MQ=255)
                                                                                                                                     aTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg                >  3:808517/1‑149 (MQ=255)
                                                                                                                                      tAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg               <  1:996770/149‑1 (MQ=255)
                                                                                                                                      tAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg               >  3:1631496/1‑149 (MQ=255)
                                                                                                                                      caaGGCCGAACGCGTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg               >  4:432169/2‑149 (MQ=255)
                                                                                                                                       aaGGCCGAACGCGTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggc              >  3:858715/1‑149 (MQ=255)
                                                                                                                                        aGGCCGAACGCGTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct             >  4:31090/1‑149 (MQ=255)
                                                                                                                                        aGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct             >  3:1070983/1‑149 (MQ=255)
                                                                                                                                          gCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGctggctgg           <  4:1161136/149‑1 (MQ=255)
                                                                                                                                                   gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTa                                                                       >  2:1348227/1‑81 (MQ=255)
                                                                                                                                                   gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTa                                                                       <  1:1348227/81‑1 (MQ=255)
                                                                                                                                                   gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                      <  1:81946/82‑1 (MQ=255)
                                                                                                                                                   gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                      >  2:81946/1‑82 (MQ=255)
                                                                                                                                                   gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATgc   <  1:1504643/149‑1 (MQ=255)
                                                                                                                                                    ttCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATgcg  >  3:748559/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCG  >  NC_000913/2622207‑2622502

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