Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,560 G→C 6.5% intergenic (‑578/+25) 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,5600GC 6.5% 227.5 / 7.1 77intergenic (‑578/+25)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base C (2/3);  total (38/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 = 7.37e-01

ATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACC  >  NC_000913/987411‑987707
                                                                                                                                                      |                                                                                                                                                   
aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCgg                                                                                                                                                     >  2:1081803/1‑150 (MQ=255)
aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTAGCAATGTAATATTCCGTAACAGGTTGCTCGTTAGGGCGGGGTTGCAAATAGATTGCGAGACTTAACAAACAGATAGCTGGACAGAAAAGTT‑AAATAAAGAGCCCGCgg                                                                                                                                                     >  2:858114/1‑150 (MQ=255)
  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                                                                                                        >  2:276773/1‑150 (MQ=255)
                                             aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                                                                        <  1:1297968/150‑1 (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:539118/150‑1 (MQ=255)
                                                                                  gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                   >  1:310441/1‑150 (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:1306274/1‑150 (MQ=255)
                                                                                                             aacaGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGt                                        >  2:1270166/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)
                                                                                                                                                      |                                                                                                                                                   
ATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACC  >  NC_000913/987411‑987707

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: ‑