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.

CTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAC  >  NC_000913/2630771‑2631069
                                                                                                                                                      |                                                                                                                                                     
cTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTgtagt                                                                                                                                                       <  1:61495/150‑1 (MQ=255)
cTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATAAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTgtagt                                                                                                                                                       >  1:1685913/1‑150 (MQ=255)
  gACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTaa                                                                                                                                                     >  2:1149281/1‑150 (MQ=255)
    cGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACa                                                                                                                                                   <  2:1513811/150‑1 (MQ=255)
    cGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGAAGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACa                                                                                                                                                   <  2:907431/150‑1 (MQ=255)
      ccAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCATATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACaga                                                                                                                                                 >  1:1142464/1‑150 (MQ=255)
        aaCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAgg                                                                                                                                               >  2:1174041/1‑150 (MQ=255)
         aCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAggg                                                                                                                                              <  1:480417/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAATCAGAAGAGAACAAATCGGACAGCTTGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAGGGt                                                                                                                                             <  1:369892/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAGGGt                                                                                                                                             <  1:222482/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAGGGt                                                                                                                                             <  2:588583/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAGGGt                                                                                                                                             <  2:65775/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGt                                                                                                                                             <  2:1632635/150‑1 (MQ=255)
          cAGAAATATGACTGAAGGCGGGCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAGGGt                                                                                                                                             <  1:203040/150‑1 (MQ=255)
               atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAGGGtttttt                                                                                                                                        <  2:693626/150‑3 (MQ=255)
               atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGt                                                                                                                                        <  2:941929/150‑1 (MQ=255)
               atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGt                                                                                                                                        <  2:948060/150‑1 (MQ=255)
                tatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCAATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGtt                                                                                                                                       <  2:91066/150‑1 (MQ=255)
                tatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAACAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAGGGtttttta                                                                                                                                       <  1:197818/150‑4 (MQ=255)
                 atGTCTGAATGCGTTCTCTATCTGATGAAAAGAGTGGGCAAATCAGAAGAGATCAAATGGGAGAGGTGTTTTAATTAATCGATGTTAGTAACTTGAATCCTTTAATTCTTTAAGATAAAATACCCTCTGTTAGTACAGAGGGttttttat                                                                                                                                      <  2:1551642/150‑5 (MQ=255)
                 atGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAGGGttttttat                                                                                                                                      <  2:436221/150‑5 (MQ=255)
                  tGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttca                                                                                                                                     >  2:678090/1‑150 (MQ=255)
                   gACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcat                                                                                                                                    >  1:890571/1‑150 (MQ=255)
                   gACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcat                                                                                                                                    <  2:602840/150‑1 (MQ=255)
                      tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcattca                                                                                                                                 >  1:839316/1‑150 (MQ=255)
                       gAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTtcattcat                                                                                                                                >  2:1270497/1‑150 (MQ=255)
                       gAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTtcattcat                                                                                                                                >  1:13579/1‑150 (MQ=255)
                             ggTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCa                                                                                                                          >  2:115367/1‑150 (MQ=255)
                              gTCGCGATCTGATGAAAAGCGTGGGCAAAGTAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAg                                                                                                                         >  2:1692660/1‑150 (MQ=255)
                              gTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAg                                                                                                                         <  1:854173/150‑1 (MQ=255)
                                cgcgATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAggg                                                                                                                       <  2:458825/150‑1 (MQ=255)
                                      cTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                                 >  1:1525423/1‑150 (MQ=255)
                                       tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATc                                                                                                                <  1:1270497/150‑1 (MQ=255)
                                       tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACACTCTGTTGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATc                                                                                                                >  1:846972/1‑150 (MQ=255)
                                        gatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCa                                                                                                               <  1:450416/150‑1 (MQ=255)
                                           gaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTc                                                                                                            >  1:20863/1‑150 (MQ=255)
                                             aaaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAATATAAAAAACCCTCTGTAGTAAAAGAGGGTTTTGTTCATTCATAGTTCAGGGTCAAATCATTccc                                                                                                          <  1:738238/150‑1 (MQ=255)
                                                       cAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGt                                                                                                >  2:113301/1‑150 (MQ=255)
                                                        aaaGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTa                                                                                               <  1:471005/150‑1 (MQ=255)
                                                          aGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAg                                                                                              <  1:573670/149‑1 (MQ=255)
                                                            cAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACACTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTg                                                                                           >  1:1214951/1‑150 (MQ=255)
                                                             agaagaGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTgg                                                                                          >  2:187708/1‑150 (MQ=255)
                                                                agagAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAAAAGAGGGTTTTTTTCATTCATATTGAAGGGAAAAATAATTCACAATAAATGGTAGATtgcgg                                                                                       >  1:18529/1‑150 (MQ=255)
                                                                  agaACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCt                                                                                     >  2:43229/1‑150 (MQ=255)
                                                                       acaTCGGACAGCTGGTTTACTTACTCGATGGTAGTCACTTCACTCCTCGACTGCTTGCAGATACACAACCCTCTGTCGTCCCAGAGGGTTTTGTGCATTCCTAGTGCCGGGGCACATCATTCCCACTCATTGGTAGCTGGCGGCTTGCCg                                                                                <  1:215481/148‑1 (MQ=255)
                                                                         aTCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAATTCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCt                                                                              <  1:115367/150‑1 (MQ=255)
                                                                         aTCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCt                                                                              <  2:72567/150‑1 (MQ=255)
                                                                            ggACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGAt                                                                           >  1:838221/1‑150 (MQ=255)
                                                                                 gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCAt                                                                      >  2:684093/1‑150 (MQ=255)
                                                                                 gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCAt                                                                      <  1:164486/150‑1 (MQ=255)
                                                                                  cTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATa                                                                     <  2:65164/150‑1 (MQ=255)
                                                                                  cTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGATGGCGGCTTGCCGCTGATGTCATa                                                                     >  2:38397/1‑150 (MQ=255)
                                                                                         aattaatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCAcg                                                              >  1:1273478/1‑150 (MQ=255)
                                                                                             aatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCggg                                                          <  2:839316/150‑1 (MQ=255)
                                                                                             aatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGATGATGTCATACAACACGCggg                                                          >  1:881399/1‑150 (MQ=255)
                                                                                                cGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGaaa                                                       >  2:1022731/1‑150 (MQ=255)
                                                                                                     ttAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCg                                                  >  1:1248083/1‑150 (MQ=255)
                                                                                                     ttAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCg                                                  >  1:1260056/1‑150 (MQ=255)
                                                                                                           aCTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACt                                            >  1:811138/1‑150 (MQ=255)
                                                                                                           aCTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCATGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACt                                            >  1:640988/1‑150 (MQ=255)
                                                                                                            cTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACAACGCGGGAAATACAGTTCACtt                                           >  1:1032572/1‑150 (MQ=255)
                                                                                                              tCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGGCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATTCACCACGCGGGAAATACCGTTCACTTCa                                         <  1:49736/150‑1 (MQ=255)
                                                                                                                aaTCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                       <  1:981950/150‑1 (MQ=255)
                                                                                                                 aTCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTg                                      <  1:574542/150‑1 (MQ=255)
                                                                                                                  tCCTATAATTCTTTAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCAACTCAATGGTAGCTGGAGGATTGCAACTTATGTCATACAACACGCGGGAAATACCGTTCACTTCATTGa                                     >  1:1694518/1‑150 (MQ=255)
                                                                                                                  tCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGa                                     <  1:1174041/150‑1 (MQ=255)
                                                                                                                   ccTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGAt                                    <  1:389502/150‑1 (MQ=255)
                                                                                                                   ccTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGAt                                    >  1:1613206/1‑150 (MQ=255)
                                                                                                                       taATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGTAAATACCGTTCACTTCATTGATAATg                                <  2:1273478/150‑1 (MQ=255)
                                                                                                                         aTTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCg                              >  1:1294205/1‑150 (MQ=255)
                                                                                                                         aTTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCg                              >  1:979477/1‑150 (MQ=255)
                                                                                                                          ttCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCgg                             <  2:1525423/150‑1 (MQ=255)
                                                                                                                            cTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                           >  2:180867/1‑150 (MQ=255)
                                                                                                                            cTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                           >  1:617449/1‑150 (MQ=255)
                                                                                                                                 atgaacaaaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATTCCGTTCACTTCATTGATAATGCAGTTGGaaa                      >  2:850545/7‑150 (MQ=255)
                                                                                                                                 aGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGaaa                      <  1:602226/150‑1 (MQ=255)
                                                                                                                                 aGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGaaa                      >  1:1163103/1‑150 (MQ=255)
                                                                                                                                  gATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc                     >  2:204333/1‑150 (MQ=255)
                                                                                                                                    tAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgc                   >  1:1335942/1‑150 (MQ=255)
                                                                                                                                     aaaaaaCCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgcg                  >  2:1316130/1‑150 (MQ=255)
                                                                                                                                     aaaaaaCCCTCTGTAGTAACAGAGGGTTTTGTTCATTAATAGTGCAGGGTCAAAACATTCACACTCAATCGTAGCTGGAGGATTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgcg                  >  1:64529/1‑150 (MQ=255)
                                                                                                                                      aaaaaCCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGa                 >  1:1216569/1‑150 (MQ=255)
                                                                                                                                      aaaaaCCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTCGGAAACGCGa                 >  1:1001790/1‑150 (MQ=255)
                                                                                                                                         aaCCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCg              >  1:1479898/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTg                          <  2:691457/137‑1 (MQ=255)
                                                                                                                                          aCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTg                          >  1:691457/1‑137 (MQ=255)
                                                                                                                                          aCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:957798/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  2:331268/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:411973/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:515587/1‑150 (MQ=255)
                                                                                                                                          aCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACACAATGCTAGCTGGCTTCTTGCCGAATATGTCATACACAACGCGGGAAATAACTTTCACATCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:283760/1‑150 (MQ=255)
                                                                                                                                            ccTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAgg           <  1:987446/150‑1 (MQ=255)
                                                                                                                                             ctctGTAGTAACATAGGGTTTTGTTCATTCATAGGGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGa          <  1:1393468/150‑1 (MQ=255)
                                                                                                                                             ctctGTAGTAACAGAGGGTTTTGTTGATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGAAGCTGATGTCATACACCACCCGGGAAATTCCGTTCACTTCATTGATAATGCGGTTGGAACCGCGACCGAGGa          >  2:300335/1‑150 (MQ=255)
                                                                                                                                                 gtagtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATc      >  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)
                                                                                                                                                      |                                                                                                                                                     
CTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAC  >  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: