Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,630,920 T→C 18.1% intergenic (+55/+38) 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,9200TC18.1% 214.1 / 44.4 105intergenic (+55/+38)yfgI/guaAnalidixic acid resistance protein, putative periplasmic protein/GMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base T (50/36);  new base C (8/11);  total (58/47)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

Only 100 of 113 total aligned reads displayed.

CTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAC  >  NC_000913/2630771‑2631069
                                                                                                                                                      |                                                                                                                                                     
cTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTg‑tagt                                                                                                                                                       <  1:61495/150‑1 (MQ=255)
cTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATAAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTg‑tagt                                                                                                                                                       >  1:1685913/1‑150 (MQ=255)
  gACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTaa                                                                                                                                                     >  2:1149281/1‑150 (MQ=255)
    cGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACa                                                                                                                                                   <  2:1513811/150‑1 (MQ=255)
    cGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGAAGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACa                                                                                                                                                   <  2:907431/150‑1 (MQ=255)
      ccAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCATATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACaga                                                                                                                                                 >  1:1142464/1‑150 (MQ=255)
        aaCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAgg                                                                                                                                               >  2:1174041/1‑150 (MQ=255)
         aCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAggg                                                                                                                                              <  1:480417/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAATCAGAAGAGAACAAATCGGACAGCTTGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  1:369892/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  1:222482/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  2:588583/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  2:65775/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGt                                                                                                                                             <  2:1632635/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGGCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  1:203040/150‑1 (MQ=255)
               atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGtttttt                                                                                                                                        <  2:693626/150‑3 (MQ=255)
               atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGt                                                                                                                                        <  2:941929/150‑1 (MQ=255)
               atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGt                                                                                                                                        <  2:948060/150‑1 (MQ=255)
                tatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCAATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGtt                                                                                                                                       <  2:91066/150‑1 (MQ=255)
                tatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAACAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGtttttta                                                                                                                                       <  1:197818/150‑4 (MQ=255)
                 atGTCTGAATGCGTTCTCTATCTGATGAAAAGAGTGGGCAAATCAGAAGAGATCAAATGGGAGAGGTGTTTTAATTAATCGATGTTAGTAACTTGAATCCTTTAATTCTTTAAGATAAAATACCCTCTGTTAGT‑ACAGAGGGttttttat                                                                                                                                      <  2:1551642/150‑5 (MQ=255)
                 atGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGttttttat                                                                                                                                      <  2:436221/150‑5 (MQ=255)
                  tGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGttca                                                                                                                                     >  2:678090/1‑150 (MQ=255)
                   gACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGttcat                                                                                                                                    >  1:890571/1‑150 (MQ=255)
                   gACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGttcat                                                                                                                                    <  2:602840/150‑1 (MQ=255)
                      tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGttcattca                                                                                                                                 >  1:839316/1‑150 (MQ=255)
                       gAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTtcattcat                                                                                                                                >  2:1270497/1‑150 (MQ=255)
                       gAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTtcattcat                                                                                                                                >  1:13579/1‑150 (MQ=255)
                             ggTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCa                                                                                                                          >  2:115367/1‑150 (MQ=255)
                              gTCGCGATCTGATGAAAAGCGTGGGCAAAGTAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAg                                                                                                                         >  2:1692660/1‑150 (MQ=255)
                              gTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAg                                                                                                                         <  1:854173/150‑1 (MQ=255)
                                cgcgATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAggg                                                                                                                       <  2:458825/150‑1 (MQ=255)
                                      cTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                                 >  1:1525423/1‑150 (MQ=255)
                                       tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATc                                                                                                                <  1:1270497/150‑1 (MQ=255)
                                       tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACACTCTG‑TTGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATc                                                                                                                >  1:846972/1‑150 (MQ=255)
                                        gatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCa                                                                                                               <  1:450416/150‑1 (MQ=255)
                                           gaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTc                                                                                                            >  1:20863/1‑150 (MQ=255)
                                             aaaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAATATAAAAAACCCTCTG‑TAGTAAAAGAGGGTTTTGTTCATTCATAGTTCAGGGTCAAATCATTccc                                                                                                          <  1:738238/150‑1 (MQ=255)
                                                       cAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGt                                                                                                >  2:113301/1‑150 (MQ=255)
                                                        aaaGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTa                                                                                               <  1:471005/150‑1 (MQ=255)
                                                          aGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAg                                                                                              <  1:573670/149‑1 (MQ=255)
                                                            cAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACACTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTg                                                                                           >  1:1214951/1‑150 (MQ=255)
                                                             agaagaGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTgg                                                                                          >  2:187708/1‑150 (MQ=255)
                                                                agagAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAAAAGAGGGTTTTTTTCATTCATATTGAAGGGAAAAATAATTCACAATAAATGGTAGATtgcgg                                                                                       >  1:18529/1‑150 (MQ=255)
                                                                  agaACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCt                                                                                     >  2:43229/1‑150 (MQ=255)
                                                                       acaTCGGACAGCTGGTTTACTTACTCGATGGTAGTCACTTCACTCCTCGACTGCTTGCAGATACACAACCCTCTG‑TCGTCCCAGAGGGTTTTGTGCATTCCTAGTGCCGGGGCACATCATTCCCACTCATTGGTAGCTGGCGGCTTGCCg                                                                                <  1:215481/148‑1 (MQ=255)
                                                                         aTCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAATTCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCt                                                                              <  1:115367/150‑1 (MQ=255)
                                                                         aTCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCt                                                                              <  2:72567/150‑1 (MQ=255)
                                                                            ggACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGAt                                                                           >  1:838221/1‑150 (MQ=255)
                                                                                 gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCAt                                                                      >  2:684093/1‑150 (MQ=255)
                                                                                 gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCAt                                                                      <  1:164486/150‑1 (MQ=255)
                                                                                  cTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATa                                                                     <  2:65164/150‑1 (MQ=255)
                                                                                  cTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGATGGCGGCTTGCCGCTGATGTCATa                                                                     >  2:38397/1‑150 (MQ=255)
                                                                                         aattaatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCAcg                                                              >  1:1273478/1‑150 (MQ=255)
                                                                                             aatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCggg                                                          <  2:839316/150‑1 (MQ=255)
                                                                                             aatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGATGATGTCATACAACACGCggg                                                          >  1:881399/1‑150 (MQ=255)
                                                                                                cGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGaaa                                                       >  2:1022731/1‑150 (MQ=255)
                                                                                                     ttAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCg                                                  >  1:1248083/1‑150 (MQ=255)
                                                                                                     ttAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCg                                                  >  1:1260056/1‑150 (MQ=255)
                                                                                                           aCTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACt                                            >  1:811138/1‑150 (MQ=255)
                                                                                                           aCTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCATGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACt                                            >  1:640988/1‑150 (MQ=255)
                                                                                                            cTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACAACGCGGGAAATACAGTTCACtt                                           >  1:1032572/1‑150 (MQ=255)
                                                                                                              tCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGGCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATTCACCACGCGGGAAATACCGTTCACTTCa                                         <  1:49736/150‑1 (MQ=255)
                                                                                                                aaTCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                       <  1:981950/150‑1 (MQ=255)
                                                                                                                 aTCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTg                                      <  1:574542/150‑1 (MQ=255)
                                                                                                                  tCCTATAATTCTTTAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCAACTCAATGGTAGCTGGAGGATTGCAACTTATGTCATACAACACGCGGGAAATACCGTTCACTTCATTGa                                     >  1:1694518/1‑150 (MQ=255)
                                                                                                                  tCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGa                                     <  1:1174041/150‑1 (MQ=255)
                                                                                                                   ccTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGAt                                    <  1:389502/150‑1 (MQ=255)
                                                                                                                   ccTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGAt                                    >  1:1613206/1‑150 (MQ=255)
                                                                                                                       taATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGTAAATACCGTTCACTTCATTGATAATg                                <  2:1273478/150‑1 (MQ=255)
                                                                                                                         aTTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCg                              >  1:1294205/1‑150 (MQ=255)
                                                                                                                         aTTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCg                              >  1:979477/1‑150 (MQ=255)
                                                                                                                          ttCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCgg                             <  2:1525423/150‑1 (MQ=255)
                                                                                                                            cTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                           >  2:180867/1‑150 (MQ=255)
                                                                                                                            cTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                           >  1:617449/1‑150 (MQ=255)
                                                                                                                                 atgaacaaaaCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATTCCGTTCACTTCATTGATAATGCAGTTGGaaa                      >  2:850545/7‑150 (MQ=255)
                                                                                                                                 aGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGaaa                      <  1:602226/150‑1 (MQ=255)
                                                                                                                                 aGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGaaa                      >  1:1163103/1‑150 (MQ=255)
                                                                                                                                  gATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc                     >  2:204333/1‑150 (MQ=255)
                                                                                                                                    tAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgc                   >  1:1335942/1‑150 (MQ=255)
                                                                                                                                     aaaaaaCCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgcg                  >  2:1316130/1‑150 (MQ=255)
                                                                                                                                     aaaaaaCCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTAATAGTGCAGGGTCAAAACATTCACACTCAATCGTAGCTGGAGGATTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgcg                  >  1:64529/1‑150 (MQ=255)
                                                                                                                                      aaaaaCCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGa                 >  1:1216569/1‑150 (MQ=255)
                                                                                                                                      aaaaaCCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTCGGAAACGCGa                 >  1:1001790/1‑150 (MQ=255)
                                                                                                                                         aaCCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCg              >  1:1479898/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTg                          <  2:691457/137‑1 (MQ=255)
                                                                                                                                          aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTg                          >  1:691457/1‑137 (MQ=255)
                                                                                                                                          aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:957798/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  2:331268/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:411973/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:515587/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACACAATGCTAGCTGGCTTCTTGCCGAATATGTCATACACAACGCGGGAAATAACTTTCACATCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:283760/1‑150 (MQ=255)
                                                                                                                                            ccTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAgg           <  1:987446/150‑1 (MQ=255)
                                                                                                                                             ctctG‑TAGTAACATAGGGTTTTGTTCATTCATAGGGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGa          <  1:1393468/150‑1 (MQ=255)
                                                                                                                                             ctctG‑TAGTAACAGAGGGTTTTGTTGATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGAAGCTGATGTCATACACCACCCGGGAAATTCCGTTCACTTCATTGATAATGCGGTTGGAACCGCGACCGAGGa          >  2:300335/1‑150 (MQ=255)
                                                                                                                                                 g‑tagtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATc      >  1:1162528/1‑150 (MQ=255)
                                                                                                                                                   tagtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCg     <  1:1149281/150‑1 (MQ=255)
                                                                                                                                                     gtaACAGAGGGTTTTGTTCATTCATAGTGAAGGCTCAAATCATTCCCACTCAATGGTAGCTGGCGGCAAGCCGGTTATGTCATACACCACGCGGGAAATATCGTTCACTTCATTGATAATGCGCTAGGAAACGCGACCGAGGAAATCGTa   >  2:536106/1‑150 (MQ=255)
                                                                                                                                                      taACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc                     <  1:1390196/131‑1 (MQ=255)
                                                                                                                                                      taACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc                     >  2:1390196/1‑131 (MQ=255)
                                                                                                                                                      taACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAc  <  2:13579/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAC  >  NC_000913/2630771‑2631069

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑