Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,348:1 +A 8.5% 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 8.5% 274.2 / 21.7 71coding (115/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base . (33/32);  new base A (3/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

CCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTT  >  NC_000913/2622200‑2622496
                                                                                                                                                     |                                                                                                                                                    
ccTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCg                                                                                                                                                       >  3:428526/1‑149 (MQ=255)
  tctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCATCGTACGGGCAGTTTTCAGCAATAAGGCCG‑Ac                                                                                                                                                    >  2:452184/1‑149 (MQ=255)
  tctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑aa                                                                                                                                                    >  2:196257/1‑148 (MQ=255)
  tctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑Ac                                                                                                                                                    >  3:32827/1‑149 (MQ=255)
       ttCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGtt                                                                                                                                               >  3:330868/1‑149 (MQ=255)
       ttCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGtt                                                                                                                                               >  4:578296/1‑149 (MQ=255)
       ttCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGtt                                                                                                                                               >  3:186518/1‑149 (MQ=255)
           gAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGc                                                                                                                                           >  2:66299/1‑149 (MQ=255)
           gAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGc                                                                                                                                           <  2:408332/149‑1 (MQ=255)
              caacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCtt                                                                                                                                        <  1:499014/149‑1 (MQ=255)
              caacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                        <  4:431232/149‑1 (MQ=255)
                 cgcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttatt                                                                                                                                     <  1:196257/149‑4 (MQ=255)
                 cgcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGaa                                                                                                                                     <  1:181731/149‑1 (MQ=255)
                 agcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGaa                                                                                                                                     <  3:20852/148‑1 (MQ=255)
                  gcgcGGTTATCGAATAATGAATATTGTTGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAc                                                                                                                                    <  4:19890/149‑1 (MQ=255)
                  gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttatt                                                                                                                                     >  2:85294/1‑145 (MQ=255)
                  gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttatt                                                                                                                                     <  1:85294/148‑4 (MQ=255)
                        ttATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGcgccc                                                                                                                              >  3:249178/1‑149 (MQ=255)
                              aataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAgg                                                                                                                        >  4:511401/1‑149 (MQ=255)
                               ataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGc                                                                                                                       <  4:368081/149‑1 (MQ=255)
                                 aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                     <  3:254144/149‑1 (MQ=255)
                                           gtggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                           <  1:284949/149‑1 (MQ=255)
                                             ggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGc                                                                                                         >  2:456865/1‑149 (MQ=255)
                                              gtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCa                                                                                                        >  2:539673/1‑149 (MQ=255)
                                              gtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCa                                                                                                        >  1:416803/1‑149 (MQ=255)
                                               tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTc                                                                                                                                              <  3:142461/110‑1 (MQ=255)
                                               tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGCACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTc                                                                                                                                              >  4:142461/1‑110 (MQ=255)
                                                gCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATa                                                                                                      <  3:352591/149‑1 (MQ=255)
                                                    aTTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCt                                                                                                  <  4:67680/149‑1 (MQ=255)
                                                     tttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTc                                                                                                 >  1:386784/1‑149 (MQ=255)
                                                                ggAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgcgc                                                                                      >  1:505447/1‑149 (MQ=255)
                                                                  aaGTAATTTACAGGCAATTATTGACGCCTGTAAAACCGACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGt                                                                                    >  3:579633/1‑149 (MQ=255)
                                                                      aaTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                <  1:539673/149‑1 (MQ=255)
                                                                          tACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGt                                                                            <  1:533090/149‑1 (MQ=255)
                                                                             aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACgcgc                                                                                                                                <  2:474762/94‑1 (MQ=255)
                                                                             aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACgcgc                                                                                                                                >  1:474762/1‑94 (MQ=255)
                                                                             aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                         <  2:416803/149‑1 (MQ=255)
                                                                                cAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                      <  3:60903/149‑1 (MQ=255)
                                                                                      ttGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCGGCGCGTTTGACAGTCGTGAAGCCTAt                                                                <  2:386784/149‑1 (MQ=255)
                                                                                        gACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGa                                                              >  3:581864/1‑149 (MQ=255)
                                                                                        gACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGa                                                              <  2:221286/149‑1 (MQ=255)
                                                                                         aCGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAc                                                             >  4:46485/1‑149 (MQ=255)
                                                                                           gCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGc                                                                    <  1:8334/140‑1 (MQ=255)
                                                                                           gCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGc                                                                    >  2:8334/1‑140 (MQ=255)
                                                                                           gCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCg                                                           >  2:228234/1‑149 (MQ=255)
                                                                                              tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                        <  3:339217/149‑1 (MQ=255)
                                                                                              tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                        >  1:292296/1‑149 (MQ=255)
                                                                                                 aaaaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCAGGAGtt                                                     <  4:72397/149‑1 (MQ=255)
                                                                                                  aaaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTg                                                    >  2:125582/1‑149 (MQ=255)
                                                                                                     ccAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAtt                                                 <  1:199680/149‑1 (MQ=255)
                                                                                                       aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgg                                                          >  3:288399/1‑138 (MQ=255)
                                                                                                       aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgg                                                          <  4:288399/138‑1 (MQ=255)
                                                                                                        acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                              <  3:194713/149‑1 (MQ=255)
                                                                                                        acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                              <  4:579633/149‑1 (MQ=255)
                                                                                                             aTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAAt                                         <  4:298833/149‑1 (MQ=255)
                                                                                                              ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTc                                                                           >  4:397734/1‑114 (MQ=255)
                                                                                                              ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTc                                                                           <  3:397734/114‑1 (MQ=255)
                                                                                                              ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                        >  2:454940/1‑149 (MQ=255)
                                                                                                              ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                        >  1:213226/1‑149 (MQ=255)
                                                                                                                aaaGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa                                      <  1:189326/149‑1 (MQ=255)
                                                                                                                aaaGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa                                      <  4:330868/149‑1 (MQ=255)
                                                                                                                  aGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACa                                    >  4:91871/1‑149 (MQ=255)
                                                                                                                    gCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATg                                  >  2:381657/1‑149 (MQ=255)
                                                                                                                      aCCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                >  4:284776/1‑101 (MQ=255)
                                                                                                                      aCCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                <  3:284776/101‑1 (MQ=255)
                                                                                                                      aCCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                                >  3:139814/1‑149 (MQ=255)
                                                                                                                         gTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGc                             >  2:345733/1‑149 (MQ=255)
                                                                                                                               gCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGa                       >  3:501574/1‑149 (MQ=255)
                                                                                                                                  gTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAtgt                    <  4:32827/149‑1 (MQ=255)
                                                                                                                                            tcaaGGCCGAACGCG‑TCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg          >  4:197869/3‑149 (MQ=255)
                                                                                                                                              aaGGCCGAACGCG‑TCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggc        >  4:447439/1‑149 (MQ=255)
                                                                                                                                                ggCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggctg      >  3:583048/1‑149 (MQ=255)
                                                                                                                                                ggCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggctg      <  3:46485/149‑1 (MQ=255)
                                                                                                                                                    g‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGGTTGACAGGCGTGAAGCCGATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTCGTCGTGCTGGCTGGttt  <  4:207990/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTT  >  NC_000913/2622200‑2622496

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