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

ATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACC  >  NC_000913/987411‑987707
                                                                                                                                                      |                                                                                                                                                   
aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                                                                                                     >  2:1081803/1‑150 (MQ=255)
aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTAGCAATGTAATATTCCGTAACAGGTTGCTCGTTAGGGCGGGGTTGCAAATAGATTGCGAGACTTAACAAACAGATAGCTGGACAGAAAAGTTAAATAAAGAGCCCGCgg                                                                                                                                                     >  2:858114/1‑150 (MQ=255)
  ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAAGAGCCAGCggg                                                                                                                                                    >  2:799861/1‑150 (MQ=255)
  ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                                                                                   >  1:878359/1‑150 (MQ=255)
       aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg                                                                                                                                              >  1:1124859/1‑150 (MQ=255)
       aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg                                                                                                                                              >  1:740862/1‑150 (MQ=255)
            aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtt                                                                                                                                         <  2:849378/150‑1 (MQ=255)
             aTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCTTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCttc                                                                                                                                        >  2:583014/1‑149 (MQ=255)
              tGCTGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACACGATGATCGTTATGGCGGGGGTGCAAATAGATTGCGTGACTTAACAAACAGATAACATGACAGAAAAGATAAAAAAACAGCCAGCGGGCGCTGGCTCTTGa                                                                                                                                       <  1:175323/150‑1 (MQ=255)
              tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                                                                       <  2:639093/150‑1 (MQ=255)
              tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                                                                       <  2:164735/150‑1 (MQ=255)
              tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttt                                                                                                                                       <  1:587659/150‑3 (MQ=255)
              tGCGGATTAGTTTTTGCGCGTAAATGTCGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTgc                                                                                                                                       >  2:796587/1‑149 (MQ=255)
               gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                                                                                                                      <  1:555983/150‑1 (MQ=255)
               gCGGATTAGTTTTTGCGCGTAAATGTGTCAATGTAATATGCCGTAACAGGATGATCGTTATGGCGGGGTTGCAACTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATTAAAAAAGAGCCAGCGCCCGCGGGCTCttttt                                                                                                                                      <  2:87844/150‑4 (MQ=255)
                        tttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaaa                                                                                                                             >  1:238313/1‑150 (MQ=255)
                          tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAtt                                                                                                                           <  1:824074/150‑1 (MQ=255)
                           ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGATGCAAATAGATTGCTTTACTTAACAAACAGATAACTTGACATAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                                                                          <  1:811465/150‑1 (MQ=255)
                            tGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTat                                                                                                                         >  1:1251183/1‑149 (MQ=255)
                             gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGa                                                                                                                        >  2:1107667/1‑150 (MQ=255)
                              cgcgTAAATGTTGCAATGTAATATTCCGAAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                                                                                       <  1:536130/150‑1 (MQ=255)
                                cgTAAATTTTGCAATGTAATATTCCGTAACACCAAGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAACATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                                                                                                                     <  1:633259/150‑1 (MQ=255)
                                  tAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTg                                                                                                                   <  1:313980/150‑1 (MQ=255)
                                     aTGTTGCAATGTAATATTCCGTAACATGATGATCGTTATGGCGGGGTTTCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCg                                                                                                                <  1:941724/150‑1 (MQ=255)
                                     aTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCg                                                                                                                <  2:1045631/150‑1 (MQ=255)
                                       gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGAAAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                                                                                              <  1:536799/150‑1 (MQ=255)
                                        ttGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAATTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGATGCTGGCGTTa                                                                                                             >  2:643694/1‑150 (MQ=255)
                                        ttGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTa                                                                                                             >  1:571111/1‑150 (MQ=255)
                                          gCAATGTAATAATCCGTAACAGGATGATCGTTATGGCGGTGTTGAAAATAGATAGCTAGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcc                                                                                                           >  2:897955/1‑149 (MQ=255)
                                             aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                                                                        >  2:276773/1‑150 (MQ=255)
                                             aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                                                                        <  1:1297968/150‑1 (MQ=255)
                                               ttAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                                                                      <  1:564143/149‑1 (MQ=255)
                                                     ttCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                                                                                                <  1:860270/150‑1 (MQ=255)
                                                       ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTAGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                                                                                              >  2:243229/1‑150 (MQ=255)
                                                         gTAACAGGAGGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCACTTGAAGACTCTAAATTACAAGCTGGCGATACGCGGAGTACGAGGGac                                                                                            >  1:869462/1‑149 (MQ=255)
                                                          tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAACCAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                                                                           <  1:173442/150‑1 (MQ=255)
                                                               ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                      >  2:313049/1‑150 (MQ=255)
                                                               ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                      >  1:118658/1‑150 (MQ=255)
                                                                gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                                                                     <  1:1100886/150‑1 (MQ=255)
                                                                 atgatCGTTATGGCGGGGGTGAAAATAGATTGCTTGACTTAACAAACACATAACTTGACAGAAAAGATAAAAAATGAGCCAGCGGGCGCTTGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                                                                    <  2:470423/150‑1 (MQ=255)
                                                                      cGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATCAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatca                                                                               <  1:236893/150‑1 (MQ=255)
                                                                       gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatcac                                                                              >  2:643410/1‑150 (MQ=255)
                                                                          aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgta                                                                           <  2:878359/150‑1 (MQ=255)
                                                                              cGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtt                                                                       <  1:183071/150‑1 (MQ=255)
                                                                                 ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                                                                    >  2:1136693/1‑150 (MQ=255)
                                                                                  gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                   <  1:539118/150‑1 (MQ=255)
                                                                                  gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                   >  1:310441/1‑150 (MQ=255)
                                                                                    tGCAAATAGATTGCTTGACTTAACACACAGATAACTTGACAGAAAAGATAAAAAAACAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                                                                 <  1:168227/150‑1 (MQ=255)
                                                                                    tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGACCGGAATCACATCACGTACGTTTTGCAc                                                                 >  2:1230204/1‑150 (MQ=255)
                                                                                    tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                                                                 <  1:847066/150‑1 (MQ=255)
                                                                                    tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                                                                 <  1:491026/150‑1 (MQ=255)
                                                                                    tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGCCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                                                                 >  2:1235109/1‑150 (MQ=255)
                                                                                     gCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                                                                >  1:450432/1‑150 (MQ=255)
                                                                                     gCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAATAGAAGCTGGCGATACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                                                                >  1:7277/1‑150 (MQ=255)
                                                                                      cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                               >  2:290836/1‑150 (MQ=255)
                                                                                      cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                               <  1:485119/150‑1 (MQ=255)
                                                                                          tAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGAGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTATGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGt                                                           >  1:802071/1‑150 (MQ=255)
                                                                                                     aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCaa                                                <  1:948949/150‑1 (MQ=255)
                                                                                                     aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGACAGTTACGTAAGCaa                                                >  1:1220168/1‑150 (MQ=255)
                                                                                                      cTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                               >  1:740680/1‑150 (MQ=255)
                                                                                                             aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGt                                        >  2:1306274/1‑150 (MQ=255)
                                                                                                             aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGt                                        >  2:1270166/1‑150 (MQ=255)
                                                                                                              acaGATAACTTGACAGAAAAGATAAAAAAAGAGCCATAGGGCGCTGGCTCTTGAAGACGCAAAAATAGAAGCTGGCGTTACTCGGAGTACGAGGGAACGGAATCACAACACGTACGATTTGCACGCCAGTTACGAAAGCAATCAGACttt                                       >  1:529784/1‑150 (MQ=255)
                                                                                                                    aaCTTGACAGAAAAGATATAAAAATAGCCAGCGGGCGCTGGCTCTTTAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAAc                                 <  2:122124/150‑1 (MQ=255)
                                                                                                                    aaCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAAc                                 <  1:610752/150‑1 (MQ=255)
                                                                                                                     aCTTGACAGAAAAGATAAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAAc                                 <  1:799861/150‑1 (MQ=255)
                                                                                                                        tgccAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAACCCatg                             >  2:894621/4‑148 (MQ=255)
                                                                                                                           cAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTTCACGCCAGTTACGTAAGCAATCAGAAGTTCAAAACCAAGAcc                          >  2:908252/1‑150 (MQ=255)
                                                                                                                           cAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGAcc                          >  1:134075/1‑150 (MQ=255)
                                                                                                                            aGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCg                         <  2:238313/150‑1 (MQ=255)
                                                                                                                                    tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGa                 <  1:290836/150‑1 (MQ=255)
                                                                                                                                    tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTAAAAACCAAGACCGAAACCTGa                 >  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)
                                                                                                                                                      |                                                                                                                                                   
ATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACC  >  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: