Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,354 Δ1 bp 7.1% 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. 7.1% 531.4 / 30.6 155coding (121/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (70/74);  new base . (5/6);  total (75/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 = 9.84e-01

Only 100 of 158 total aligned reads displayed.

TTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCG  >  NC_000913/2622207‑2622502
                                                                                                                                                    |                                                                                                                                                    
ttCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGtt                                                                                                                                                     <  1:907956/149‑1 (MQ=255)
 tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTc                                                                                                                                                    <  4:96819/149‑1 (MQ=255)
  ccGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCg                                                                                                                                                   >  3:1243674/1‑149 (MQ=255)
    gAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGc                                                                                                                                                 >  1:1355895/1‑149 (MQ=255)
     aacaacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGcc                                                                                                                                                >  4:119098/1‑149 (MQ=255)
      acaacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCt                                                                                                                                               <  3:576824/149‑1 (MQ=255)
        aacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTg                                                                                                                                             >  2:517523/1‑149 (MQ=255)
        aacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTg                                                                                                                                             >  2:470768/1‑149 (MQ=255)
           gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAAc                                                                                                                                          <  1:521118/149‑1 (MQ=255)
           gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAAc                                                                                                                                          <  3:76095/149‑1 (MQ=255)
           gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAAc                                                                                                                                          >  1:238508/1‑149 (MQ=255)
               ggTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTAACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACgcgc                                                                                                                                      >  3:136206/1‑149 (MQ=255)
                 ttATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGcgccc                                                                                                                                    >  3:310092/1‑149 (MQ=255)
                 ttATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGcgccc                                                                                                                                    >  1:742306/1‑149 (MQ=255)
                   aTCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGcgcccgc                                                                                                                                  >  4:172527/1‑149 (MQ=255)
                      gAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCtta                                                                                                                                             >  1:208420/1‑134 (MQ=255)
                      gAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCtta                                                                                                                                             <  2:208420/135‑2 (MQ=255)
                      gAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAg                                                                                                                               >  2:1573183/1‑149 (MQ=255)
                       aataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAgg                                                                                                                              >  3:371421/1‑149 (MQ=255)
                       aataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAgg                                                                                                                              >  3:853734/1‑149 (MQ=255)
                         taatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCg                                                                                                                            >  2:1034841/1‑149 (MQ=255)
                           atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAAcgc                                                                                                                                        >  4:1489266/1‑135 (MQ=255)
                           atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCggg                                                                                                                          >  4:199894/1‑149 (MQ=255)
                             gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                        >  1:575612/1‑149 (MQ=255)
                               atatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTAtt                                                                                                                      >  2:355381/1‑149 (MQ=255)
                                tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTg                                                                                                                     >  2:28685/1‑149 (MQ=255)
                                 atTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                                    <  1:1353619/149‑1 (MQ=255)
                                    gtggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                                 <  4:371421/149‑1 (MQ=255)
                                         gCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCtt                                                                                                                                              <  3:1234688/115‑1 (MQ=255)
                                          cTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                           <  2:1205678/133‑1 (MQ=255)
                                          cTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                           >  1:1205678/1‑133 (MQ=255)
                                           ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACGAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                        <  4:483216/135‑1 (MQ=255)
                                           ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACGAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                        >  3:483216/1‑135 (MQ=255)
                                           ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGc                                                                                                                             >  1:772361/1‑130 (MQ=255)
                                           ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGc                                                                                                                             <  2:772361/130‑1 (MQ=255)
                                            tATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGc                                                                                                         <  1:133850/149‑1 (MQ=255)
                                              tttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAATGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGTCGGGTATTGCAACGCATACGCTc                                                                                                       >  3:18513/1‑149 (MQ=255)
                                              tttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTc                                                                                                       >  4:1470825/1‑149 (MQ=255)
                                                    gCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGcc                                                                                                 <  4:364367/149‑1 (MQ=255)
                                                         ggAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgcgc                                                                                            <  1:687181/149‑1 (MQ=255)
                                                          gAAGTAATTTACAGGCAATTATTGACGCGTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGcgcg                                                                                           >  3:1347361/1‑149 (MQ=255)
                                                          gAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACATAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCtt                                                                                                                                              >  3:961730/1‑98 (MQ=255)
                                                          gAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCAGGTATTGCAACGCATACGCTCATCGCCAGcgcg                                                                                           >  2:622280/1‑149 (MQ=255)
                                                           aaGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGt                                                                                          >  1:1249570/1‑149 (MQ=255)
                                                             gTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGttt                                                                                        <  2:521351/149‑1 (MQ=255)
                                                                aTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGAc                                                                                     >  3:751599/1‑149 (MQ=255)
                                                                      aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                               <  2:1479651/149‑1 (MQ=255)
                                                                      aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                               <  2:125506/149‑1 (MQ=255)
                                                                          aaTTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGcgcccgc                                                                                                                                  <  4:233452/94‑1 (MQ=255)
                                                                          aaTTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGcgcccgc                                                                                                                                  >  3:233452/1‑94 (MQ=255)
                                                                          aaTTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                            <  2:1117491/148‑1 (MQ=255)
                                                                          aaTTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                            >  1:1117491/1‑148 (MQ=255)
                                                                            ttattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGcc                                                                         <  1:355381/149‑1 (MQ=255)
                                                                             tattGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                                        <  3:1062237/149‑1 (MQ=255)
                                                                               ttGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGTGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                      <  1:1042138/149‑1 (MQ=255)
                                                                               tagACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                      >  3:869055/3‑149 (MQ=255)
                                                                                     ccTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgg                                                                <  4:1055044/149‑1 (MQ=255)
                                                                                     ccTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgg                                                                <  1:961319/149‑1 (MQ=255)
                                                                                      cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCggg                                                               >  1:406810/1‑149 (MQ=255)
                                                                                       tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                              <  1:1093033/149‑1 (MQ=255)
                                                                                       tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTAAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                      <  1:1451069/141‑1 (MQ=255)
                                                                                        gTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                                 >  3:155990/1‑97 (MQ=255)
                                                                                        gTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                                 <  4:155990/97‑1 (MQ=255)
                                                                                           aaaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTg                                                          >  1:731903/1‑149 (MQ=255)
                                                                                            aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCTGGAGTTGa                                                         <  4:136206/149‑1 (MQ=255)
                                                                                            aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                         <  4:192754/149‑1 (MQ=255)
                                                                                            aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                         <  3:645943/149‑1 (MQ=255)
                                                                                             aCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAt                                                        <  1:1262861/149‑1 (MQ=255)
                                                                                                aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                                     >  4:69199/1‑149 (MQ=255)
                                                                                                 acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATGAGGCCGACGCGTTCGGCCTTGAACGCGGCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                                    >  3:976417/1‑149 (MQ=255)
                                                                                                  caaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATg                                                   >  4:1301749/1‑149 (MQ=255)
                                                                                                   aaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGa                                                  >  1:124321/1‑149 (MQ=255)
                                                                                                     aaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGaaa                                                <  4:271083/149‑1 (MQ=255)
                                                                                                      aTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAAt                                               <  2:575612/149‑1 (MQ=255)
                                                                                                          aaGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGAc                                           <  2:742306/149‑1 (MQ=255)
                                                                                                          aaGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGAc                                           <  2:1067794/149‑1 (MQ=255)
                                                                                                            ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACAt                                         <  4:751599/149‑1 (MQ=255)
                                                                                                            ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCATGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACAt                                         >  3:1032851/1‑149 (MQ=255)
                                                                                                                ccGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTAc                                     >  3:1483209/1‑149 (MQ=255)
                                                                                                                ccGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTAc                                     <  3:244187/149‑1 (MQ=255)
                                                                                                                      gggCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCAccc                               <  2:96186/149‑1 (MQ=255)
                                                                                                                       ggCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                                    <  4:1556501/127‑1 (MQ=255)
                                                                                                                       ggCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                                    >  3:1556501/1‑127 (MQ=255)
                                                                                                                       ggCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCg                              >  1:1113284/1‑149 (MQ=255)
                                                                                                                        gCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGtt                                                                                         <  4:81156/89‑1 (MQ=255)
                                                                                                                        gCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGtt                                                                                         >  3:81156/1‑89 (MQ=255)
                                                                                                                         cAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAt                            <  1:1387571/149‑1 (MQ=255)
                                                                                                                              ttCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGt                       <  2:124321/149‑1 (MQ=255)
                                                                                                                                 aGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGt                    <  1:496543/149‑1 (MQ=255)
                                                                                                                                 aGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGt                    >  2:261042/1‑149 (MQ=255)
                                                                                                                                  gCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTg                   >  2:807047/1‑149 (MQ=255)
                                                                                                                                   cAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgc                  <  3:389101/149‑1 (MQ=255)
                                                                                                                                   cAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgc                  <  3:1498320/149‑1 (MQ=255)
                                                                                                                                    ttcaaGGCCGAACGCGTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTAAGCACCCGATGTGGTCGTgct                 >  4:696907/4‑149 (MQ=255)
                                                                                                                                     aTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg                >  1:71955/1‑149 (MQ=255)
                                                                                                                                       aaGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggc              <  2:1355895/149‑1 (MQ=255)
                                                                                                                                        aGGCCGAACGCGTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct             >  3:1501427/1‑149 (MQ=255)
                                                                                                                                        aGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct             <  3:199894/149‑1 (MQ=255)
                                                                                                                                           ccGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGt          >  4:1427843/1‑149 (MQ=255)
                                                                                                                                             gACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGttt        <  1:1573183/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCG  >  NC_000913/2622207‑2622502

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