Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,354 Δ1 bp 5.5% 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. 5.5% 621.3 / 25.8 182coding (121/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (77/93);  new base . (3/7);  total (81/101)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

Only 100 of 184 total aligned reads displayed.

ATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGC  >  NC_000913/2622206‑2622501
                                                                                                                                                     |                                                                                                                                                   
aTTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGt                                                                                                                                                     >  3:102605/1‑149 (MQ=255)
aTTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGt                                                                                                                                                     >  3:333175/1‑149 (MQ=255)
  tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTc                                                                                                                                                   <  2:1639259/149‑1 (MQ=255)
    cGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCgg                                                                                                                                                 >  3:780312/1‑149 (MQ=255)
     gAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGc                                                                                                                                                >  4:500109/1‑149 (MQ=255)
      aacaacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGcc                                                                                                                                               <  2:323756/149‑1 (MQ=255)
       acaacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCt                                                                                                                                              <  3:1312573/149‑1 (MQ=255)
        caacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                             <  1:513107/149‑1 (MQ=255)
         aacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTg                                                                                                                                            >  1:1703869/1‑149 (MQ=255)
          acGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttat                                                                                                                                           <  3:1324128/149‑3 (MQ=255)
             cgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttattgc                                                                                                                                        <  2:873876/149‑6 (MQ=255)
              gTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAcgc                                                                                                                                       <  3:807989/149‑1 (MQ=255)
                     tCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgcccgcc                                                                                                                                >  1:1505610/1‑149 (MQ=255)
                     tCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgcccgcc                                                                                                                                >  1:24275/1‑149 (MQ=255)
                      cGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCa                                                                                                                               <  2:776551/149‑1 (MQ=255)
                       gAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAg                                                                                                                              >  3:693850/1‑149 (MQ=255)
                          taatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCg                                                                                                                           >  1:1298557/1‑149 (MQ=255)
                           aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                          <  1:302748/149‑1 (MQ=255)
                           aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                          <  1:997053/149‑1 (MQ=255)
                           aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                          <  4:1010806/149‑1 (MQ=255)
                            atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCggg                                                                                                                         >  2:45154/1‑149 (MQ=255)
                            atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCggg                                                                                                                         >  2:528296/1‑149 (MQ=255)
                             tGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgtgt                                                                                                                        <  2:1089441/149‑1 (MQ=255)
                              gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                       >  1:12060/1‑149 (MQ=255)
                               aaTATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTAt                                                                                                                      >  1:802375/1‑149 (MQ=255)
                                atatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGCTCGGCCTTGAACGCGCCCGCCAGGCGGGTAtt                                                                                                                     <  2:1495728/149‑1 (MQ=255)
                                 tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTg                                                                                                                    >  1:66033/1‑149 (MQ=255)
                                   ttGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCa                                                                                                                  <  1:1351066/149‑1 (MQ=255)
                                     gtggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                                <  4:658663/149‑1 (MQ=255)
                                           cTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATAc                                                                                                          <  2:24275/149‑1 (MQ=255)
                                               tttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTc                                                                                                      >  4:1612616/1‑149 (MQ=255)
                                                  ccGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATc                                                                                                   <  1:855947/149‑1 (MQ=255)
                                                          ggAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgcgc                                                                                           <  3:1145649/149‑1 (MQ=255)
                                                          ggAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgcgc                                                                                           <  4:570680/149‑1 (MQ=255)
                                                          ggAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgcgc                                                                                           >  4:991179/1‑149 (MQ=255)
                                                           gAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGcgcg                                                                                          >  2:809111/1‑149 (MQ=255)
                                                           gAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGcgcg                                                                                          >  1:1477390/1‑149 (MQ=255)
                                                            aaGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGt                                                                                         >  4:1316093/1‑149 (MQ=255)
                                                             aGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGtt                                                                                        <  2:518108/149‑1 (MQ=255)
                                                               tAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttat                                                                                                                                           >  4:495553/1‑94 (MQ=255)
                                                                aaTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                     <  3:689055/149‑1 (MQ=255)
                                                                  tttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACa                                                                                   <  3:1612616/149‑1 (MQ=255)
                                                                   ttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAg                                                                                  >  2:392656/1‑149 (MQ=255)
                                                                    tACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGt                                                                                 <  2:1314977/149‑1 (MQ=255)
                                                                     aCAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTc                                                                                <  1:1359876/149‑1 (MQ=255)
                                                                       aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                              >  3:1503849/1‑149 (MQ=255)
                                                                         gCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGa                                                                            >  3:288769/1‑149 (MQ=255)
                                                                         gCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGa                                                                            >  2:62214/1‑149 (MQ=255)
                                                                            attattGACGCCTGTAACACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCACTAAGGCCG‑ACGTGTTCGCCCTTGACCGCGCCGGCCAGGCGCGTACTGCAACGCATACGCTCATCGCCAGCGCCTTTGACCGTCGTGAAGc                                                                         <  4:941399/149‑1 (MQ=255)
                                                                            attattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGc                                                                         <  4:6695/149‑1 (MQ=255)
                                                                             ttattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGcc                                                                        <  3:1244621/149‑1 (MQ=255)
                                                                             ttattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTAATGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGcc                                                                        >  4:1436037/1‑149 (MQ=255)
                                                                              tattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCTCCAGCGCGTTTGACAGTCGTGAAGCCt                                                                       <  3:1155154/149‑1 (MQ=255)
                                                                               attGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGACCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTa                                                                      <  2:919256/149‑1 (MQ=255)
                                                                               attGACGCCTGTAAAACCAACAAAATTAAAGGCACAGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTa                                                                      >  4:1624689/1‑149 (MQ=255)
                                                                                 tGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATg                                                                    >  1:1003895/1‑149 (MQ=255)
                                                                                     gCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCg                                                                >  3:31604/1‑149 (MQ=255)
                                                                                      ccTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgg                                                               <  3:1693074/149‑1 (MQ=255)
                                                                                        tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                             >  3:462002/1‑149 (MQ=255)
                                                                                           aaaaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCTCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGtt                                                          <  4:333175/149‑1 (MQ=255)
                                                                                             aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                        >  4:183887/1‑149 (MQ=255)
                                                                                             aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                        >  2:1531563/1‑149 (MQ=255)
                                                                                             aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                        <  4:1441988/149‑1 (MQ=255)
                                                                                             aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                        <  4:31604/149‑1 (MQ=255)
                                                                                             aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                        <  3:883329/149‑1 (MQ=255)
                                                                                              aCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                                >  2:796725/1‑92 (MQ=255)
                                                                                               ccAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAtt                                                      <  1:1740209/149‑1 (MQ=255)
                                                                                                caacaaAATTAAAGGCACCGTACAGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTc                                                     >  2:572931/1‑149 (MQ=255)
                                                                                                 aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                           >  1:1041355/1‑126 (MQ=255)
                                                                                                 aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                           <  2:1041355/126‑1 (MQ=255)
                                                                                                  actaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATACGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGGTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATGCAt                                                   <  4:127185/146‑1 (MQ=255)
                                                                                                  acaaAATTAAAGGAACGGTACGCGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGCCGTGAAGCCAATGACCGGGAGTTGATTCAt                                                   >  2:110284/1‑149 (MQ=255)
                                                                                                    aaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGc                                                                         <  3:366970/125‑1 (MQ=255)
                                                                                                       aTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAAt                                              >  1:1071230/1‑149 (MQ=255)
                                                                                                        ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                             <  4:1678691/149‑1 (MQ=255)
                                                                                                        ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                             <  2:1298557/149‑1 (MQ=255)
                                                                                                        ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                             >  4:148378/1‑149 (MQ=255)
                                                                                                          aaaGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa                                           <  1:45154/149‑1 (MQ=255)
                                                                                                            aGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACa                                         <  2:1467983/149‑1 (MQ=255)
                                                                                                             ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACAt                                        >  4:433958/1‑149 (MQ=255)
                                                                                                                   gTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGc                                  >  3:498120/1‑149 (MQ=255)
                                                                                                                    tACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCa                                 <  2:889710/149‑1 (MQ=255)
                                                                                                                        ggCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCggg                                                              <  1:1742171/116‑1 (MQ=255)
                                                                                                                        ggCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCggg                                                              >  2:1742171/1‑116 (MQ=255)
                                                                                                                        ggCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGAc                                          <  3:1192387/136‑1 (MQ=255)
                                                                                                                        ggCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGAc                                          >  4:1192387/1‑136 (MQ=255)
                                                                                                                         gCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGa                            >  3:1458129/1‑149 (MQ=255)
                                                                                                                         gCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGa                            >  3:1353829/1‑149 (MQ=255)
                                                                                                                            gTTTTCAGCGATAAGGCCT‑ACGCGTTCGGCCTGGAACGCGCCCGCCAGGCAGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAAGGGTCATGTACGCACCCGAtgt                         <  2:1693717/149‑1 (MQ=255)
                                                                                                                             ttttCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAtgtg                        >  2:1730517/1‑149 (MQ=255)
                                                                                                                             ttttCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAtgtg                        >  1:72073/1‑149 (MQ=255)
                                                                                                                                tCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTc                     <  1:1572361/149‑1 (MQ=255)
                                                                                                                                 cAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCAACATGTACGCACCCGATGTGGTCg                    >  2:976313/1‑149 (MQ=255)
                                                                                                                                  aGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGt                   <  2:1477390/149‑1 (MQ=255)
                                                                                                                                   gCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTg                  >  4:1676899/1‑149 (MQ=255)
                                                                                                                                             cG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGtt        <  3:1408872/149‑1 (MQ=255)
                                                                                                                                             cG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGtt        >  1:110317/1‑149 (MQ=255)
                                                                                                                                                    gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATgc  <  2:12060/149‑1 (MQ=255)
                                                                                                                                                    gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATgc  <  1:483302/149‑1 (MQ=255)
                                                                                                                                                    gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCAGCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTACGGGCTGGTTTTATgc  <  2:125220/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGC  >  NC_000913/2622206‑2622501

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