Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,348:1 +A 7.0% coding (115/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,3481.A 7.0% 645.6 / 44.0 158coding (115/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base . (73/74);  new base A (5/6);  total (78/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 161 total aligned reads displayed.

CTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTT  >  NC_000913/2622201‑2622496
                                                                                                                                                    |                                                                                                                                                    
 tctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑aa                                                                                                                                                    >  3:446266/1‑148 (MQ=255)
      ttCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGtt                                                                                                                                               <  1:907956/149‑1 (MQ=255)
      ttCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGtt                                                                                                                                               >  1:1593805/1‑149 (MQ=255)
       tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTc                                                                                                                                              <  4:96819/149‑1 (MQ=255)
        ccGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCg                                                                                                                                             >  3:1243674/1‑149 (MQ=255)
          gAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAACTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGc                                                                                                                                           >  2:961319/1‑149 (MQ=255)
          gAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGc                                                                                                                                           >  1:1355895/1‑149 (MQ=255)
           aacaacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGcc                                                                                                                                          >  4:119098/1‑149 (MQ=255)
            acaacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCt                                                                                                                                         <  3:576824/149‑1 (MQ=255)
             caacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                        <  1:485927/149‑1 (MQ=255)
              aacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTg                                                                                                                                       >  2:470768/1‑149 (MQ=255)
              aacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTg                                                                                                                                       >  2:517523/1‑149 (MQ=255)
               acGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGa                                                                                                                                      >  4:508998/1‑149 (MQ=255)
                 gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAc                                                                                                                                    <  3:508998/149‑1 (MQ=255)
                 gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAc                                                                                                                                    <  1:722329/149‑1 (MQ=255)
                 gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAc                                                                                                                                    <  1:1530500/149‑1 (MQ=255)
                 gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAc                                                                                                                                    <  1:521118/149‑1 (MQ=255)
                   gTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAcgc                                                                                                                                  <  2:1593805/149‑1 (MQ=255)
                     ggTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTAACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACgcgc                                                                                                                                >  3:136206/1‑149 (MQ=255)
                       ttATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgccc                                                                                                                              >  1:742306/1‑149 (MQ=255)
                       ttATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgccc                                                                                                                              >  3:310092/1‑149 (MQ=255)
                         aTCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgcccgc                                                                                                                            >  4:172527/1‑149 (MQ=255)
                         aTCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgcccgc                                                                                                                            >  1:96186/1‑149 (MQ=255)
                            gAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCtta                                                                                                                                       >  1:208420/1‑134 (MQ=255)
                             aataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAgg                                                                                                                        >  3:853734/1‑149 (MQ=255)
                             aataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAgg                                                                                                                        >  3:371421/1‑149 (MQ=255)
                               taatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCg                                                                                                                      >  2:1034841/1‑149 (MQ=255)
                                 atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAcgc                                                                                                                                  <  3:1489266/135‑1 (MQ=255)
                                 atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAcgc                                                                                                                                  >  4:1489266/1‑135 (MQ=255)
                                 atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCggg                                                                                                                    >  4:199894/1‑149 (MQ=255)
                                       atTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                              <  1:1353619/149‑1 (MQ=255)
                                          gtggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                           <  4:371421/149‑1 (MQ=255)
                                          gtggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                           <  2:238508/149‑1 (MQ=255)
                                            ggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGc                                                                                                         <  1:745882/149‑1 (MQ=255)
                                               gCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCtt                                                                                                                                        >  4:1234688/1‑115 (MQ=255)
                                                cTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                     >  1:1205678/1‑133 (MQ=255)
                                                 ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACGAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                  <  4:483216/135‑1 (MQ=255)
                                                 ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACg                                                                                                    >  1:125506/1‑149 (MQ=255)
                                                 ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACg                                                                                                    >  4:204644/1‑149 (MQ=255)
                                                    tttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAATGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGTCGGGTATTGCAACGCATACGCTc                                                                                                 >  3:18513/1‑149 (MQ=255)
                                                          gCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGcc                                                                                           <  4:364367/149‑1 (MQ=255)
                                                           cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTGCGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCa                                                                                          <  1:1341716/149‑1 (MQ=255)
                                                           cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCa                                                                                          <  4:18513/149‑1 (MQ=255)
                                                                gAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACATAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                        >  3:961730/1‑98 (MQ=255)
                                                                gAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCAGGTATTGCAACGCATACGCTCATCGCCAGcgcg                                                                                     >  2:622280/1‑149 (MQ=255)
                                                                      aTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGAc                                                                               >  3:751599/1‑149 (MQ=255)
                                                                            aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                         <  2:1479651/149‑1 (MQ=255)
                                                                            aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                         <  2:125506/149‑1 (MQ=255)
                                                                            aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                         <  1:519231/149‑1 (MQ=255)
                                                                               cAATTATTGACGCCTGTAAAACCAACAACATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                      >  3:1004095/1‑149 (MQ=255)
                                                                               cAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                      <  2:767043/149‑1 (MQ=255)
                                                                               cAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                      <  1:1248325/149‑1 (MQ=255)
                                                                                aaTTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgcccgc                                                                                                                            >  3:233452/1‑94 (MQ=255)
                                                                                aaTTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                      <  2:1117491/148‑1 (MQ=255)
                                                                                aaTTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                      >  1:1117491/1‑148 (MQ=255)
                                                                                aaTTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAg                                                                     >  4:1439080/1‑149 (MQ=255)
                                                                                  ttattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGcc                                                                   <  1:355381/149‑1 (MQ=255)
                                                                                   tattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                                  <  3:1062237/149‑1 (MQ=255)
                                                                                   tattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                                  <  2:1249570/149‑1 (MQ=255)
                                                                                   tattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                                  <  3:106849/149‑1 (MQ=255)
                                                                                     ttGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGTGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                <  1:1042138/149‑1 (MQ=255)
                                                                                        aCGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAc                                                             <  2:1117315/149‑1 (MQ=255)
                                                                                          gCCTGTAAAACCAACACAATTAAAGGCACCGTACGGGCAGTTTTCAGCCATACGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCg                                                           >  3:933512/1‑149 (MQ=255)
                                                                                          gCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCg                                                           >  2:555411/1‑149 (MQ=255)
                                                                                            cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTa                                                                 >  1:3337/1‑141 (MQ=255)
                                                                                            cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTa                                                                 <  2:3337/141‑1 (MQ=255)
                                                                                             tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                >  2:1451069/1‑141 (MQ=255)
                                                                                             tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                        <  1:1093033/149‑1 (MQ=255)
                                                                                             tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTAAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                <  1:1451069/141‑1 (MQ=255)
                                                                                              gTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                           <  4:155990/97‑1 (MQ=255)
                                                                                              gTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                           >  3:155990/1‑97 (MQ=255)
                                                                                                  aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCTGGAGTTGa                                                   <  4:136206/149‑1 (MQ=255)
                                                                                                   aCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAt                                                  <  1:1262861/149‑1 (MQ=255)
                                                                                                   aCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAt                                                  >  4:1018422/1‑149 (MQ=255)
                                                                                                         aaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGa                                            >  1:124321/1‑149 (MQ=255)
                                                                                                           aaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGaaa                                          <  4:271083/149‑1 (MQ=255)
                                                                                                            aTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAAt                                         <  2:575612/149‑1 (MQ=255)
                                                                                                             ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGACCGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                        <  3:696632/149‑1 (MQ=255)
                                                                                                             ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                        >  3:1506842/1‑149 (MQ=255)
                                                                                                                aaGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGAc                                     <  2:742306/149‑1 (MQ=255)
                                                                                                                 aGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACa                                    <  3:864307/149‑1 (MQ=255)
                                                                                                                  ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACAt                                   <  4:751599/149‑1 (MQ=255)
                                                                                                                   gCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATg                                  >  3:81127/1‑149 (MQ=255)
                                                                                                                     aCCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                                <  1:608143/149‑1 (MQ=255)
                                                                                                                      ccGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTAc                               >  3:1483209/1‑149 (MQ=255)
                                                                                                                      ccGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTAc                               <  3:244187/149‑1 (MQ=255)
                                                                                                                             ggCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCg                        >  1:1113284/1‑149 (MQ=255)
                                                                                                                              gCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGtt                                                                                   <  4:81156/89‑1 (MQ=255)
                                                                                                                              gCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGa                       >  2:123500/1‑149 (MQ=255)
                                                                                                                               cAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAt                      <  1:1387571/149‑1 (MQ=255)
                                                                                                                                       aGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGt              >  2:261042/1‑149 (MQ=255)
                                                                                                                                        gCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTg             >  4:1498320/1‑149 (MQ=255)
                                                                                                                                        gCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTg             >  2:807047/1‑149 (MQ=255)
                                                                                                                                         cAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgc            <  3:389101/149‑1 (MQ=255)
                                                                                                                                          ttcaaGGCCGAACGCG‑TCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTAAGCACCCGATGTGGTCGTgct           >  4:696907/4‑149 (MQ=255)
                                                                                                                                           aTAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg          >  1:71955/1‑149 (MQ=255)
                                                                                                                                             aaGGCCGAACGCG‑TCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggc        >  3:1308314/1‑149 (MQ=255)
                                                                                                                                              aGGCCGAACGCG‑TCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct       >  3:1501427/1‑149 (MQ=255)
                                                                                                                                              aGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct       <  3:199894/149‑1 (MQ=255)
                                                                                                                                                 ccG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGt    >  4:1427843/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTT  >  NC_000913/2622201‑2622496

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: ‑