Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,630,919 G→A 5.6% intergenic (+54/+39) 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,9190GA 5.6% 356.9 / 7.8 106intergenic (+54/+39)yfgI/guaAnalidixic acid resistance protein, putative periplasmic protein/GMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base G (56/44);  new base A (4/2);  total (60/46)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 106 total aligned reads displayed.

ACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAA  >  NC_000913/2630774‑2631063
                                                                                                                                                 |                                                                                                                                                
aCGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTaa                                                                                                                                               >  3:843263/1‑149 (MQ=255)
     aaCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACagag                                                                                                                                          >  2:384579/1‑149 (MQ=255)
     aaCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACagag                                                                                                                                          >  2:936011/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                                                                                                                              <  3:34539/149‑1 (MQ=255)
                 aCTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcat                                                                                                                              <  2:804478/149‑1 (MQ=255)
                  cTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcatt                                                                                                                             <  4:796969/149‑1 (MQ=255)
                   tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAGGAGAACAAATCGGACAGCTGGTTTAATTAATCGGTGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcattc                                                                                                                            >  4:467440/1‑149 (MQ=255)
                   tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACaga                                                                                                                                           >  3:50633/1‑134 (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                                                                                                           >  1:98949/1‑149 (MQ=255)
                                    tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                           >  2:883045/1‑149 (MQ=255)
                                    tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                           >  3:832331/1‑149 (MQ=255)
                                    tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                           >  3:134388/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                                                                                                      >  1:813093/1‑145 (MQ=255)
                                             gCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTc                                                                                                      <  2:813093/145‑1 (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)
                                                                gaACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTAAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTGTAACAGAGGGTTTTGTTCATTCATATTACAGGGTCAAATCATTCACACTCAATGGTAGCTGTCGACt                                                                               >  2:916199/1‑149 (MQ=255)
                                                                  aCAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTg                                                                             >  1:178622/1‑149 (MQ=255)
                                                                    aaaTCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGcc                                                                           >  2:753255/1‑149 (MQ=255)
                                                                          gACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGAt                                                                     >  3:863612/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                                                   >  4:372027/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                                     >  3:164191/1‑149 (MQ=255)
                                                                                                          ttCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACtt                                     >  4:715533/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                                                                                                                                  >  2:354593/1‑46 (MQ=255)
                                                                                                                    taATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGt                                                                                                                                  <  1:354593/46‑1 (MQ=255)
                                                                                                                         cTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGt                      <  2:389045/149‑1 (MQ=255)
                                                                                                                         cTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGt                      <  4:353945/149‑1 (MQ=255)
                                                                                                                          ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGATATACCGTTCACTTCATTGATAATGCGGtt                     <  3:551977/149‑1 (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)
                                                                                                                          ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                     >  1:784974/1‑149 (MQ=255)
                                                                                                                          ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                     >  1:654034/1‑149 (MQ=255)
                                                                                                                            gAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTgg                   >  1:299959/1‑149 (MQ=255)
                                                                                                                              atgaacaaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCgg                       <  2:170546/137‑1 (MQ=39)
                                                                                                                              atgaacaaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCgg                       >  1:170546/7‑143 (MQ=39)
                                                                                                                               gATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                 >  1:687664/1‑132 (MQ=255)
                                                                                                                               gATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                 <  2:687664/132‑1 (MQ=255)
                                                                                                                                aTAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc               <  1:745170/149‑1 (MQ=255)
                                                                                                                                 aacaaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATAc                                              <  4:774191/114‑1 (MQ=38)
                                                                                                                                 aacaaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATAc                                              >  3:774191/4‑117 (MQ=39)
                                                                                                                                    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)
                                                                                                                                             tGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGaaa  >  1:410927/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAA  >  NC_000913/2630774‑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: