Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,561 G→C 6.3% intergenic (‑579/+24) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,5610GC 6.3% 221.5 / 5.3 76intergenic (‑579/+24)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base C (2/3);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

CCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTA  >  NC_000913/987413‑987710
                                                                                                                                                     |                                                                                                                                                     
ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAAGAGCCAGCggg                                                                                                                                                       >  2:799861/1‑150 (MQ=255)
ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGgc                                                                                                                                                      >  1:878359/1‑150 (MQ=255)
     aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTgg                                                                                                                                                 >  1:1124859/1‑150 (MQ=255)
     aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTgg                                                                                                                                                 >  1:740862/1‑150 (MQ=255)
          aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGCCCGCTGGCTCtt                                                                                                                                            <  2:849378/150‑1 (MQ=255)
           aTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCTTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCttc                                                                                                                                           >  2:583014/1‑149 (MQ=255)
            tGCTGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACACGATGATCGTTATGGCGGGGGTGCAAATAGATTGCGTGACTTAACAAACAGATAACATGACAGAAAAGAT‑AAAAAAACAGCCAGCGGGCGCTGGCTCTTGa                                                                                                                                          <  1:175323/150‑1 (MQ=255)
            tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                                                                          <  2:639093/150‑1 (MQ=255)
            tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                                                                          <  2:164735/150‑1 (MQ=255)
            tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGCCCGCTGGCTCtttt                                                                                                                                          <  1:587659/150‑3 (MQ=255)
            tGCGGATTAGTTTTTGCGCGTAAATGTCGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTgc                                                                                                                                          >  2:796587/1‑149 (MQ=255)
             gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                                                                                                                         <  1:555983/150‑1 (MQ=255)
             gCGGATTAGTTTTTGCGCGTAAATGTGTCAATGTAATATGCCGTAACAGGATGATCGTTATGGCGGGGTTGCAACTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑TAAAAAAGAGCCAGCGCCCGCGGGCTCttttt                                                                                                                                         <  2:87844/150‑4 (MQ=255)
                      tttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaaa                                                                                                                                >  1:238313/1‑150 (MQ=255)
                        tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAtt                                                                                                                              <  1:824074/150‑1 (MQ=255)
                         ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGATGCAAATAGATTGCTTTACTTAACAAACAGATAACTTGACATAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                                                                             <  1:811465/150‑1 (MQ=255)
                          tGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTat                                                                                                                            >  1:1251183/1‑149 (MQ=255)
                           gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGa                                                                                                                           >  2:1107667/1‑150 (MQ=255)
                            cgcgTAAATGTTGCAATGTAATATTCCGAAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                                                                                          <  1:536130/150‑1 (MQ=255)
                              cgTAAATTTTGCAATGTAATATTCCGTAACACCAAGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAACAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                                                                                                                        <  1:633259/150‑1 (MQ=255)
                                tAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTg                                                                                                                      <  1:313980/150‑1 (MQ=255)
                                   aTGTTGCAATGTAATATTCCGTAACATGATGATCGTTATGGCGGGGTTTCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCg                                                                                                                   <  1:941724/150‑1 (MQ=255)
                                   aTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCg                                                                                                                   <  2:1045631/150‑1 (MQ=255)
                                     gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGAAAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                                                                                                 <  1:536799/150‑1 (MQ=255)
                                      ttGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAATTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGATGCTGGCGTTa                                                                                                                >  2:643694/1‑150 (MQ=255)
                                      ttGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTa                                                                                                                >  1:571111/1‑150 (MQ=255)
                                        gCAATGTAATAATCCGTAACAGGATGATCGTTATGGCGGTGTTGAAAATAGATAGCTAGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcc                                                                                                              >  2:897955/1‑149 (MQ=255)
                                           aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                                                                           <  1:1297968/150‑1 (MQ=255)
                                           aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                                                                           >  2:276773/1‑150 (MQ=255)
                                             ttAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                                                                         <  1:564143/149‑1 (MQ=255)
                                                   ttCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                                                                                                   <  1:860270/150‑1 (MQ=255)
                                                     ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTAGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                                                                                 >  2:243229/1‑150 (MQ=255)
                                                       gTAACAGGAGGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCACTTGAAGACTCTAAATTACAAGCTGGCGATACGCGGAGTACGAGGGac                                                                                               >  1:869462/1‑149 (MQ=255)
                                                        tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAACCAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                                                                              <  1:173442/150‑1 (MQ=255)
                                                             ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                         >  2:313049/1‑150 (MQ=255)
                                                             ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                         >  1:118658/1‑150 (MQ=255)
                                                              gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                                                                        <  1:1100886/150‑1 (MQ=255)
                                                               atgatCGTTATGGCGGGGGTGAAAATAGATTGCTTGACTTAACAAACACATAACTTGACAGAAAAGAT‑AAAAAATGAGCCAGCGGGCGCTTGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                                                                       <  2:470423/150‑1 (MQ=255)
                                                                    cGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATCAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatca                                                                                  <  1:236893/150‑1 (MQ=255)
                                                                     gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatcac                                                                                 >  2:643410/1‑150 (MQ=255)
                                                                        aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgta                                                                              <  2:878359/150‑1 (MQ=255)
                                                                            cGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtt                                                                          <  1:183071/150‑1 (MQ=255)
                                                                               ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                                                                       >  2:1136693/1‑150 (MQ=255)
                                                                                gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                      >  1:310441/1‑150 (MQ=255)
                                                                                gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                      <  1:539118/150‑1 (MQ=255)
                                                                                  tGCAAATAGATTGCTTGACTTAACACACAGATAACTTGACAGAAAAGAT‑AAAAAAACAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                                                                    <  1:168227/150‑1 (MQ=255)
                                                                                  tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGACCGGAATCACATCACGTACGTTTTGCAc                                                                    >  2:1230204/1‑150 (MQ=255)
                                                                                  tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                                                                    <  1:847066/150‑1 (MQ=255)
                                                                                  tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                                                                    <  1:491026/150‑1 (MQ=255)
                                                                                  tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGCCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                                                                    >  2:1235109/1‑150 (MQ=255)
                                                                                   gCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                                                                   >  1:450432/1‑150 (MQ=255)
                                                                                   gCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAATAGAAGCTGGCGATACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                                                                   >  1:7277/1‑150 (MQ=255)
                                                                                    cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                                  >  2:290836/1‑150 (MQ=255)
                                                                                    cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                                  <  1:485119/150‑1 (MQ=255)
                                                                                        tAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGAGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTATGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGt                                                              >  1:802071/1‑150 (MQ=255)
                                                                                                   aCTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCaa                                                   <  1:948949/150‑1 (MQ=255)
                                                                                                   aCTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGACAGTTACGTAAGCaa                                                   >  1:1220168/1‑150 (MQ=255)
                                                                                                    cTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                                  >  1:740680/1‑150 (MQ=255)
                                                                                                           aacaGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGt                                           >  2:1270166/1‑150 (MQ=255)
                                                                                                           aacaGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGt                                           >  2:1306274/1‑150 (MQ=255)
                                                                                                            acaGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCATAGGGCGCTGGCTCTTGAAGACGCAAAAATAGAAGCTGGCGTTACTCGGAGTACGAGGGAACGGAATCACAACACGTACGATTTGCACGCCAGTTACGAAAGCAATCAGACttt                                          >  1:529784/1‑150 (MQ=255)
                                                                                                                  aaCTTGACAGAAAAGAT‑ATAAAAATAGCCAGCGGGCGCTGGCTCTTTAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAAc                                    <  2:122124/150‑1 (MQ=255)
                                                                                                                  aaCTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAAc                                    <  1:610752/150‑1 (MQ=255)
                                                                                                                   aCTTGACAGAAAAGATAAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAAc                                    <  1:799861/150‑1 (MQ=255)
                                                                                                                      tgccAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAACCCatg                                >  2:894621/4‑148 (MQ=255)
                                                                                                                         cAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTTCACGCCAGTTACGTAAGCAATCAGAAGTTCAAAACCAAGAcc                             >  2:908252/1‑150 (MQ=255)
                                                                                                                         cAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGAcc                             >  1:134075/1‑150 (MQ=255)
                                                                                                                          aGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCg                            <  2:238313/150‑1 (MQ=255)
                                                                                                                                  t‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGa                    <  1:290836/150‑1 (MQ=255)
                                                                                                                                  t‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTAAAAACCAAGACCGAAACCTGa                    >  2:123134/1‑150 (MQ=255)
                                                                                                                                    tcttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCCAGACCGAAACCTGaa                   >  2:1031687/6‑150 (MQ=255)
                                                                                                                                    aaaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTGTTGCACGCCAGTTACGTAAGCAATCAGAAGTGCAAAACCAAGACCGAAAccgtaa                   >  1:727820/1‑146 (MQ=255)
                                                                                                                                     aaaacaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTTAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAAt                  <  1:1143474/150‑1 (MQ=255)
                                                                                                                                      ttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTAACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAGCAGACGTTCAAAACCACGAACGAAACCTGAATg                 >  1:337436/4‑150 (MQ=255)
                                                                                                                                      aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATg                 <  2:191445/150‑1 (MQ=255)
                                                                                                                                         aaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCgg              <  2:1121923/150‑1 (MQ=255)
                                                                                                                                              ccAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAg         <  1:888188/150‑1 (MQ=255)
                                                                                                                                                  cGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGCCGTTCAAAACCAAGACCGAAACCGGAATGCGGAACAGTAcc     >  2:329008/1‑150 (MQ=255)
                                                                                                                                                     gcgcTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTa  <  1:706839/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACCGTA  >  NC_000913/987413‑987710

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