Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,354 Δ1 bp 5.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. 5.1% 891.6 / 36.8 253coding (121/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (119/119);  new base . (6/7);  total (125/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

Only 100 of 258 total aligned reads displayed.

CACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCGCATTCTCAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCG  >  NC_000913/2622084‑2622618
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                        
cACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTACTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGAGTTCg                                                                                                                                                                                                                                                                       >  2:113255/1‑274 (MQ=255)
                                                      cAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTCTTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACCCATACGCTCATCCCCAgcac                                                                                                                                                                                                                >  4:357449/1‑273 (MQ=255)
                                                             cTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                                                                                                                                          <  2:94348/274‑1 (MQ=255)
                                                                                aCGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAcc                                                                                                                                                                                      >  2:27934/1‑275 (MQ=255)
                                                                                     gAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                                                                                                                                                  <  2:484663/274‑1 (MQ=255)
                                                                                                             atcaAAGCCTCTGATTCCGAAAAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                                                                                                                                                          <  2:332900/274‑1 (MQ=255)
                                                                                                             atcaAAGCCTCTGATCCCTAACAATGCGTGGTTATCGAATATAGAATATTGTGGGGCTTATGTCCGGAAACGGAAGGAAGTAACAGGTAATTATTGACGCCTGTAAAAACAACAAAATTAAAGGCACCGTACGGGCAGGTTTCAGCAATAAGGCCGACGCGTTCGGCCTGGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGAGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                                                                                                                                                          <  4:234663/274‑1 (MQ=255)
                                                                                                                       ctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAtgt                                                                                                                                              >  1:159492/1‑276 (MQ=255)
                                                                                                                           ttCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGtt                                                                                                                                                                                                                                                                         >  5:1202395/1‑149 (MQ=255)
                                                                                                                            tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCa                                                                                                                                                      <  2:541559/263‑1 (MQ=255)
                                                                                                                              cGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCgg                                                                                                                                                                                                                                                                      <  5:468972/149‑1 (MQ=255)
                                                                                                                              cGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCgg                                                                                                                                                                                                                                                                      >  7:1251353/1‑149 (MQ=255)
                                                                                                                                 acaacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCt                                                                                                                                                                                                                                                                   <  6:1080718/149‑1 (MQ=255)
                                                                                                                                  caacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCtt                                                                                                                                                                                                                                                                  <  7:1107569/149‑1 (MQ=255)
                                                                                                                                   aacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCtta                                                                                                                                                                                                                                                                 <  7:275647/149‑2 (MQ=255)
                                                                                                                                     cgcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGaa                                                                                                                                                                                                                                                               <  6:1233726/149‑1 (MQ=255)
                                                                                                                                          ggTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGACGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACgcgc                                                                                                                                                                                                                                                          <  7:1453024/149‑1 (MQ=255)
                                                                                                                                             tATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGcgcccg                                                                                                                                                                                                                                                       >  6:1349763/1‑149 (MQ=255)
                                                                                                                                              aTCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                                                                                                                                                        >  2:37010/1‑227 (MQ=255)
                                                                                                                                                 gAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAACAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAg                                                                                                                                                                                                                                                   >  7:1161136/1‑149 (MQ=255)
                                                                                                                                                  aataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAgg                                                                                                                                                                                                                                                  >  7:789214/1‑149 (MQ=255)
                                                                                                                                                     aatGAATATTGTGGTGCTTATTTCCGGCAACGTAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                                                                                                                                               <  6:1033387/149‑1 (MQ=255)
                                                                                                                                                      atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCggg                                                                                                                                                                                                                                              <  6:153901/149‑1 (MQ=255)
                                                                                                                                                      atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCGGGCTGGTTTTATGCGCATtct                                                                                                                >  2:625039/1‑275 (MQ=255)
                                                                                                                                                       tGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGt                                                                                                                                                                                                                                             >  6:658012/1‑149 (MQ=255)
                                                                                                                                                        gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACTCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                                                                                                                                            <  7:877654/149‑1 (MQ=255)
                                                                                                                                                        gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                                                                                                                                            >  5:1363177/1‑149 (MQ=255)
                                                                                                                                                         aaTATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTAt                                                                                                                                                                                                                                           <  6:957148/149‑1 (MQ=255)
                                                                                                                                                          atatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTAtt                                                                                                                                                                                                                                          <  8:1226953/149‑1 (MQ=255)
                                                                                                                                                           tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttat                                                                                                                                                                                                                                                                >  6:309204/1‑124 (MQ=255)
                                                                                                                                                           tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttat                                                                                                                                                                                                                                                                <  5:309204/126‑3 (MQ=255)
                                                                                                                                                            atTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                                                                                                                                                        >  6:1528444/1‑149 (MQ=255)
                                                                                                                                                              tgtgGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTg                                                                                                                                                                                                                                                                 >  2:655941/1‑122 (MQ=255)
                                                                                                                                                              tgtgGTGCNTATTTCCGNNAACGGAAGTAATTTACAGGCAATNATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                                                                                                                                                        <  3:64881/147‑1 (MQ=255)
                                                                                                                                                               gtggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                                                                                                                                                     <  8:1021555/149‑1 (MQ=255)
                                                                                                                                                                   tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCAt                                                                                                                                                                                                                                 >  7:577716/1‑149 (MQ=255)
                                                                                                                                                                   tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCAt                                                                                                                                                                                                                                 >  8:704770/1‑149 (MQ=255)
                                                                                                                                                                     cTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATAc                                                                                                                                                                                                                               >  6:1622087/1‑149 (MQ=255)
                                                                                                                                                                      ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCt                                                                                                                                                                                                                                                                   >  7:454175/1‑112 (MQ=255)
                                                                                                                                                                      ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCt                                                                                                                                                                                                                                                                   <  8:454175/112‑1 (MQ=255)
                                                                                                                                                                          ttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCCACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATTTACGCACCCGATGTGTTCGGTCTGGCTGGTTTTATGCGCATTTTCAGCCCGGCGTTTGtccc                                                                                              >  4:433231/1‑271 (MQ=255)
                                                                                                                                                                            ccGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATc                                                                                                                                                                                                                        <  8:577716/149‑1 (MQ=255)
                                                                                                                                                                              ggCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCGCATTCTCAGCCCGGCGTTTGtc                                                                                                <  2:433235/267‑1 (MQ=255)
                                                                                                                                                                                cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATc                                                                                                                                                                                                                        >  2:689566/1‑145 (MQ=255)
                                                                                                                                                                                cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCa                                                                                                                                                                                                                    <  5:618801/149‑1 (MQ=255)
                                                                                                                                                                                 aaCGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgc                                                                                                                                                                                                                  >  6:996770/1‑149 (MQ=255)
                                                                                                                                                                                   cGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATa                                                                                                                                                                                                                                <  2:212459/134‑1 (MQ=255)
                                                                                                                                                                                    ggAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCa                                                                                                                                                                                                                    <  4:689561/145‑1 (MQ=255)
                                                                                                                                                                                      aaGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGt                                                                                                                                                                                                              <  8:1288651/149‑1 (MQ=255)
                                                                                                                                                                                         tAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTg                                                                                                                                                                                                           >  8:950132/1‑149 (MQ=255)
                                                                                                                                                                                           aTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGttt                                                                                                                                                                                                            >  7:361747/1‑146 (MQ=255)
                                                                                                                                                                                           aTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGAc                                                                                                                                                                                                         >  5:784633/1‑149 (MQ=255)
                                                                                                                                                                                            tttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACa                                                                                                                                                                                                        <  6:1363177/149‑1 (MQ=255)
                                                                                                                                                                                            tttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACa                                                                                                                                                                                                        <  5:1622087/149‑1 (MQ=255)
                                                                                                                                                                                                 aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                                                                                                                                                   >  8:266574/1‑149 (MQ=255)
                                                                                                                                                                                                 aGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                                                                                                                                                   <  6:757726/149‑1 (MQ=255)
                                                                                                                                                                                                   gCAATTATTGACGCCTGTAAACCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttatt                                                                                                                                                                                                                                                               >  6:58568/1‑84 (MQ=38)
                                                                                                                                                                                                   gCAATTATTGAAGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                                                                                                                                          >  7:958947/1‑140 (MQ=255)
                                                                                                                                                                                                       ttattGACGCCCGTAAAACCAACAAAACAAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTTGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACTCCCATCGCTCGCCCGTTTGCCAGTCGCGAAGCCTCAGACCGGGAGATGATCCATGAAATCGACCTGTACTCACCCGAAGGGATCGTGCCTGCGGTGTTTAGTCGCATTTCCAACCCGCGGTTTTCCTCCCACTATCCCGGGCTTTTGccgcccagtcac                                                               >  4:572746/1‑264 (MQ=255)
                                                                                                                                                                                                             aCGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGtt                                                                                                                                                                                                                                                                         >  1:44082/1‑67 (MQ=255)
                                                                                                                                                                                                             aCGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAc                                                                                                                                                                                       <  5:901803/149‑1 (MQ=255)
                                                                                                                                                                                                                ccTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgg                                                                                                                                                                                    <  8:605901/149‑1 (MQ=255)
                                                                                                                                                                                                                 cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCtt                                                                                                                                                                                                                                                                  <  6:494494/70‑1 (MQ=255)
                                                                                                                                                                                                                  tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                                                                                                                                                  <  6:833344/149‑1 (MQ=255)
                                                                                                                                                                                                                  tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                                                                                                                                                  <  7:421587/149‑1 (MQ=255)
                                                                                                                                                                                                                  tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                                                                                                                                                  >  8:421587/1‑149 (MQ=255)
                                                                                                                                                                                                                  tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACg                                                                                                                                                        <  1:185860/175‑1 (MQ=255)
                                                                                                                                                                                                                    tAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTc                                                                                                                                                                                                     >  5:96564/1‑128 (MQ=255)
                                                                                                                                                                                                                       aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGa                                                                                                                                                                             <  7:1002829/149‑1 (MQ=255)
                                                                                                                                                                                                                       aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATAGCCAGCGCGTTTGACAGTCGTGAAGCTTATGACCGGGAGTTGa                                                                                                                                                                             <  8:39637/149‑1 (MQ=255)
                                                                                                                                                                                                                           aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGc                                                                                                                                                                                              >  6:813101/1‑128 (MQ=255)
                                                                                                                                                                                                                           aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                                                                                                                                                         >  5:1031823/1‑149 (MQ=255)
                                                                                                                                                                                                                           aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                                                                                                                                                         >  5:855492/1‑149 (MQ=255)
                                                                                                                                                                                                                            acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATAc                                                                                                                                                                                                                               <  5:1301407/94‑1 (MQ=255)
                                                                                                                                                                                                                            acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGt                                                                                                                                                                                >  6:869250/1‑141 (MQ=255)
                                                                                                                                                                                                                               aaaTTAAAGGCACCGTACGGGCAGTTTTCATCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                                                                                                                                                   <  8:1160664/119‑1 (MQ=255)
                                                                                                                                                                                                                                aaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATa                                                                                                                                                                                                                                <  4:749726/89‑1 (MQ=255)
                                                                                                                                                                                                                                  ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                                                                                                                                                  >  7:719542/1‑149 (MQ=255)
                                                                                                                                                                                                                                   caagggtgccGTACGGGTAGTTTTTAGCAATAGNCCCGGCGCGTTCGGCCTTGACCGCGNCCGCCCAGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCNTGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCGCATTCTCAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTaca                               <  3:163085/271‑1 (MQ=255)
                                                                                                                                                                                                                                     aaGGCGCCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGACCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAGTCATGAAATCGat                                                                                                                                                               <  8:1090867/149‑2 (MQ=255)
                                                                                                                                                                                                                                       ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAg                                                                                                                                                                                 <  6:433599/129‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAg                                                                                                                                                                                 >  5:433599/1‑129 (MQ=255)
                                                                                                                                                                                                                                       ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGtt                                                                                                                                                                               <  7:1432830/131‑1 (MQ=255)
                                                                                                                                                                                                                                          aCCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa                                                                                                                                                                <  5:191199/143‑1 (MQ=255)
                                                                                                                                                                                                                                              tACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCa                                                                                                                                                      >  8:1585324/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   gCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGa                                                                                                                                                 >  8:1629710/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    cAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAt                                                                                                                                                <  5:59160/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAtt                                                                                                                                                                           <  8:870714/114‑1 (MQ=255)
                                                                                                                                                                                                                                                              cAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgc                                                                                                                                      >  7:1646497/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aaTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgct                                                                                                                                     >  6:247558/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg                                                                                                                                    >  8:147878/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg                                                                                                                                   >  7:1631496/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aaGGCCGAACGCGTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggc                                                                                                                                  >  7:858715/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aGGCCGAACGCGTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct                                                                                                                                 >  8:31090/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    ggCCCACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGCATTGCAACGCATACGCTCCTCGCCCGCGCCATTGACAGCCGTGAAACCTATGACCGGGAGTTGATTCATGACATCGACATGCACGCACCCGATGTGGTCGTGCGCGCTTGCTTTATGCGCATGCTCAGCACGGCGTCTGTCTCCCACCATACCGAGCGTTTGATGCACAGTCCCCCTTCTCGGCTCCCGACATCACCCGGTTTCCATACCCATCCTCATGCGCTGGAGATTGGca  >  4:517341/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                     gCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg                                                                                                                                   <  4:154902/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ccGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGt                                                                                                                              >  8:1367685/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                                                                                                                                          <  5:81946/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATgc                                                                                                                       <  5:1504643/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATgc                                                                                                                       <  7:162368/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                        
CACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCGCATTCTCAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCG  >  NC_000913/2622084‑2622618

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