Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,630,918 A→T 6.5% intergenic (+53/+40) yfgI → / ← guaA nalidixic acid resistance protein, putative periplasmic protein/GMP synthetase (glutamine aminotransferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,630,9180AT 6.5% 377.3 / 9.2 106intergenic (+53/+40)yfgI/guaAnalidixic acid resistance protein, putative periplasmic protein/GMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base A (55/44);  new base T (5/2);  total (60/46)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAA  >  NC_000913/2630770‑2631063
                                                                                                                                                    |                                                                                                                                                 
gCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTgta                                                                                                                                                   <  2:456136/149‑1 (MQ=255)
gCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTgta                                                                                                                                                   <  4:255789/149‑1 (MQ=255)
    aCGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTaa                                                                                                                                               >  3:843263/1‑149 (MQ=255)
         aaCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACagag                                                                                                                                          >  2:384579/1‑149 (MQ=255)
         aaCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACagag                                                                                                                                          >  2:936011/1‑149 (MQ=255)
         aaCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACagag                                                                                                                                          >  4:282444/1‑149 (MQ=255)
               aaTATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGtttt                                                                                                                                    <  4:802772/149‑1 (MQ=255)
               aaTATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGtttt                                                                                                                                    <  1:25388/149‑1 (MQ=255)
                atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTg                                                                                                                                   >  4:377984/1‑149 (MQ=255)
                 tatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGt                                                                                                                                  <  4:213617/149‑1 (MQ=255)
                     aCTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcat                                                                                                                              <  2:804478/149‑1 (MQ=255)
                     aCTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcat                                                                                                                              <  3:34539/149‑1 (MQ=255)
                      cTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcatt                                                                                                                             <  4:796969/149‑1 (MQ=255)
                       tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAGGAGAACAAATCGGACAGCTGGTTTAATTAATCGGTGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcattc                                                                                                                            >  4:467440/1‑149 (MQ=255)
                       tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACaga                                                                                                                                           <  4:50633/134‑1 (MQ=255)
                       tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcattc                                                                                                                            >  1:494204/1‑149 (MQ=255)
                               gTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCa                                                                                                                    <  3:355323/149‑1 (MQ=255)
                               gTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCa                                                                                                                    >  1:253811/1‑149 (MQ=255)
                                tCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAg                                                                                                                   <  2:546972/149‑1 (MQ=255)
                                 cgcgATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAgg                                                                                                                  <  3:91772/149‑1 (MQ=255)
                                  gcgATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAggg                                                                                                                 >  4:341522/1‑149 (MQ=255)
                                    gATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTc                                                                                                               <  1:115165/149‑1 (MQ=255)
                                        tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                           >  3:832331/1‑149 (MQ=255)
                                        tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                           >  3:134388/1‑149 (MQ=255)
                                        tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                           >  2:883045/1‑149 (MQ=255)
                                        tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                           >  1:98949/1‑149 (MQ=255)
                                        tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGATTAATTAATCGATGTTAGTAACTTCAATCCTATAATGCTTGAAGACAAAAAACCTTCTGTAGTAACAGAGGGTTTTGTTCATTCGTAGTGCAGGGTCAAAt                                                                                                           >  4:551977/1‑149 (MQ=255)
                                         gatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATc                                                                                                          >  4:293243/1‑149 (MQ=255)
                                            gaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCAtt                                                                                                       <  3:718105/149‑1 (MQ=255)
                                            gaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCAtt                                                                                                       <  1:476105/149‑1 (MQ=255)
                                             aaaaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTc                                                                                                      <  3:523624/149‑1 (MQ=255)
                                               aaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTccc                                                                                                    <  2:481718/149‑1 (MQ=255)
                                                 gCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTc                                                                                                      <  2:813093/145‑1 (MQ=255)
                                                 gCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTc                                                                                                      >  1:813093/1‑145 (MQ=255)
                                                    tGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCa                                                                                               >  2:273276/1‑149 (MQ=255)
                                                       gCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATg                                                                                            >  4:867683/1‑149 (MQ=255)
                                                          aaGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTa                                                                                         <  3:775439/149‑1 (MQ=255)
                                                           aGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAACCATTCCAATTAAATgttgg                                                                                        >  2:794582/1‑145 (MQ=255)
                                                            gCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGc                                                                                       >  3:200036/1‑149 (MQ=255)
                                                               gaagaGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTgg                                                                                    >  2:157818/1‑149 (MQ=255)
                                                                 agagAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTggcg                                                                                  >  2:745170/1‑149 (MQ=255)
                                                                  gagaACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCATTCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTggcgg                                                                                 >  1:645411/1‑149 (MQ=255)
                                                                  gagaACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTggcgg                                                                                 >  2:290519/1‑149 (MQ=255)
                                                                    gaACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTAAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTGTAACAGAGGGTTTTGTTCATTCATATTACAGGGTCAAATCATTCACACTCAATGGTAGCTGTCGACt                                                                               >  2:916199/1‑149 (MQ=255)
                                                                      aCAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTg                                                                             >  1:178622/1‑149 (MQ=255)
                                                                        aaaTCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGcc                                                                           >  2:753255/1‑149 (MQ=255)
                                                                              gACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGAt                                                                     <  1:478804/149‑1 (MQ=255)
                                                                                  gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCa                                                                 <  2:494204/149‑1 (MQ=255)
                                                                                   cTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCAt                                                                >  1:934715/1‑149 (MQ=255)
                                                                                            ttaatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCAcgc                                                       <  3:467440/149‑1 (MQ=255)
                                                                                              aatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCgg                                                     <  2:178622/149‑1 (MQ=255)
                                                                                               atCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCggg                                                    >  4:268056/1‑149 (MQ=255)
                                                                                                tCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGa                                                   >  1:332768/1‑149 (MQ=255)
                                                                                                tCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGa                                                   >  3:670722/1‑149 (MQ=255)
                                                                                                  gATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGaaa                                                 >  1:357162/1‑149 (MQ=255)
                                                                                                      ttAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATAcc                                             <  1:273276/149‑1 (MQ=255)
                                                                                                      ttAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATAcc                                             <  2:937649/149‑1 (MQ=255)
                                                                                                        aGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGt                                           <  2:465739/149‑1 (MQ=255)
                                                                                                         gTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGtt                                          <  3:920137/149‑1 (MQ=255)
                                                                                                          tAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTc                                         >  4:827833/1‑149 (MQ=255)
                                                                                                           aaCTTCATTCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCa                                        <  2:645411/149‑1 (MQ=255)
                                                                                                            aCTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCAc                                       <  1:883045/149‑1 (MQ=255)
                                                                                                            aCTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCAc                                       <  1:228130/149‑1 (MQ=255)
                                                                                                              ttCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACtt                                     >  4:715533/1‑149 (MQ=255)
                                                                                                              ttCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACtt                                     >  3:164191/1‑149 (MQ=255)
                                                                                                                 aaTCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAt                                  >  1:496983/1‑149 (MQ=255)
                                                                                                                 aaTCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAt                                  <  3:827833/149‑1 (MQ=255)
                                                                                                                 aaTCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAt                                  <  1:157818/149‑1 (MQ=255)
                                                                                                                   tCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTg                                >  3:837570/1‑149 (MQ=255)
                                                                                                                      tataATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATa                             >  4:12008/1‑149 (MQ=255)
                                                                                                                       ataATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATaa                            <  2:300467/149‑1 (MQ=255)
                                                                                                                        taATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGt                                                                                                                                  <  1:354593/46‑1 (MQ=255)
                                                                                                                        taATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGt                                                                                                                                  >  2:354593/1‑46 (MQ=255)
                                                                                                                             cTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGt                      <  4:353945/149‑1 (MQ=255)
                                                                                                                              ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGATATACCGTTCACTTCATTGATAATGCGGtt                     <  3:551977/149‑1 (MQ=255)
                                                                                                                              ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                     >  1:654034/1‑149 (MQ=255)
                                                                                                                              ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                     <  3:867683/149‑1 (MQ=255)
                                                                                                                              ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                     <  1:531146/149‑1 (MQ=255)
                                                                                                                              ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                     <  3:372027/149‑1 (MQ=255)
                                                                                                                              ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                     <  1:290519/149‑1 (MQ=255)
                                                                                                                              ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATATCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                     >  3:580794/1‑149 (MQ=255)
                                                                                                                                gAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTgg                   >  1:299959/1‑149 (MQ=255)
                                                                                                                                  atgaacaaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCgg                       >  1:170546/7‑143 (MQ=39)
                                                                                                                                  atgaacaaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCgg                       <  2:170546/137‑1 (MQ=39)
                                                                                                                                   gATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                 <  2:687664/132‑1 (MQ=255)
                                                                                                                                   gATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                 >  1:687664/1‑132 (MQ=255)
                                                                                                                                    aTAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc               <  1:745170/149‑1 (MQ=255)
                                                                                                                                     aacaaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATAc                                              >  3:774191/4‑117 (MQ=39)
                                                                                                                                     aacaaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATAc                                              <  4:774191/114‑1 (MQ=38)
                                                                                                                                        aaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGa           >  1:691472/1‑149 (MQ=255)
                                                                                                                                         aaaCCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGAc          >  4:118896/1‑149 (MQ=255)
                                                                                                                                          aaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACGCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGAcc         >  2:846085/1‑149 (MQ=255)
                                                                                                                                          aaCCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGAcc         >  1:733330/1‑149 (MQ=255)
                                                                                                                                            ccccCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa       >  4:138136/1‑149 (MQ=255)
                                                                                                                                            cccTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa       >  1:596686/1‑149 (MQ=255)
                                                                                                                                             ccTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAg      >  3:545007/1‑149 (MQ=255)
                                                                                                                                               tctGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGa    >  2:722028/1‑149 (MQ=255)
                                                                                                                                               tctGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGa    >  4:830320/1‑149 (MQ=255)
                                                                                                                                                ctGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCTCTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAATGCGAACGAGGaa   <  4:580794/149‑1 (MQ=255)
                                                                                                                                                ctGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGaa   <  2:299959/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAA  >  NC_000913/2630770‑2631063

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