Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,364 G→T 6.6% intergenic (+21/+26) slyB → / ← slyA outer membrane lipoprotein/global transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,3640GT 6.6% 1215.4 / 16.3 385intergenic (+21/+26)slyB/slyAouter membrane lipoprotein/global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (221/123);  new base T (2/39);  total (223/162)
Fisher's exact test for biased strand distribution p-value = 3.81e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 402 total aligned reads displayed.

CTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCT  >  NC_000913/1720215‑1720513
                                                                                                                                                         |                                                                                                                                                     
   gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCCGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGttc                                                                                                                                                    >  2:2710112/1‑149 (MQ=255)
   gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTa                                                                                                                                                    <  2:395537/150‑1 (MQ=255)
   gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTa                                                                                                                                                    <  2:4592005/150‑1 (MQ=255)
   gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTa                                                                                                                                                    <  1:6253028/150‑1 (MQ=255)
     aaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACt                                                                                                                                                  >  1:2406753/1‑150 (MQ=255)
       ttCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGa                                                                                                                                                <  1:8693644/150‑1 (MQ=255)
        tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               <  1:7451522/150‑1 (MQ=255)
        tCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                                                                                                                                               <  1:5648386/150‑1 (MQ=255)
          gTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACca                                                                                                                                             <  2:6471918/150‑1 (MQ=255)
          gTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACca                                                                                                                                             <  1:3662798/150‑1 (MQ=255)
           tAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcac                                                                                                                                            >  2:8557882/1‑150 (MQ=255)
           tAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcac                                                                                                                                            <  2:2525663/150‑1 (MQ=255)
            aaaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcaca                                                                                                                                           <  1:203822/150‑1 (MQ=255)
                      gTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTc                                                                                                                                 >  1:10035088/1‑150 (MQ=255)
                        aaTACCTTCATGGTGGTACATAAACAAGGCAACAGTCGTTTCTCGCGGGGCAACGTTGTCTTACTGGCCTGCAATGGCAGTCTGGTGACCGTTTCTCCGCGCTAAATACTTTTGCGTGTGTTCTGGTTACTGACCCCACGCCCCCTTCAt                                                                                                                               <  2:7393875/150‑1 (MQ=255)
                          tACCATCAGGTGTGTACAGAACCACGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCCGTCAGGTGACGGTTTCTCCGCGCTAATTAAGTTTGCGTGTGGTCAGGTTAGTGACCACACGCCCCCTTCAtgt                                                                                                                             <  2:7542032/150‑3 (MQ=255)
                             catcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAc                                                                                                                          >  2:7625776/1‑150 (MQ=255)
                               tcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAccc                                                                                                                        >  2:1260752/1‑150 (MQ=255)
                               tcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAccc                                                                                                                        >  2:1864244/1‑150 (MQ=255)
                               tcatGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCAccc                                                                                                                        >  2:5146267/1‑150 (MQ=255)
                                 ttgttggtACATAAACATGGAAACACTCTTTTCTCTCGGGGCAAAGTTGTGTTACTGGCAAGAATTGGCTGTCTGGTGACGTTTTCTCCGCGCTAAATATTTTTGGGTGTGGTCAGTTTACTGACCACACGCCCCCTTCATTTCACCCtt                                                                                                                      <  2:2314426/146‑1 (MQ=255)
                                  tggtggTACAGAAACAAGGCAACCCTGGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCATTTGCTGTCAGGTGACCGTTTCTCCGCGCTAAATATGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCttt                                                                                                                     <  2:1474100/150‑1 (MQ=255)
                                    gtggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTgg                                                                                                                   >  1:7641273/1‑150 (MQ=255)
                                     tggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGc                                                                                                                  >  1:1642446/1‑150 (MQ=255)
                                     gggtACTGACACCCGGCACCACTCGTTTCTGTCGGGGCCAAGGTGGCGTAGGGGCCTGCCATGGCTGTCAGGTGAGCGTTTCGCCGGGTTACCTATGTGTGGGTGGGGGCAGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGc                                                                                                                  <  2:3358635/149‑1 (MQ=255)
                                        gtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCt                                                                                                                >  2:7837005/1‑150 (MQ=255)
                                        gtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATATCACCCTTTGGCCt                                                                                                                >  1:512448/1‑150 (MQ=255)
                                         tACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTg                                                                                                               >  1:4089904/1‑150 (MQ=255)
                                          aCAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGt                                                                                                              >  2:5502540/1‑150 (MQ=255)
                                           cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTa                                                                                                             >  2:2812972/1‑150 (MQ=255)
                                            aGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTaa                                                                                                            >  1:6432069/1‑150 (MQ=255)
                                            aGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTaa                                                                                                            >  1:3236942/1‑150 (MQ=255)
                                            aGAAACAAGGCAACACTCGGTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACGGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTac                                                                                                            >  1:2647732/1‑149 (MQ=255)
                                              aaaCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACt                                                                                                          >  2:10221030/1‑150 (MQ=255)
                                                aCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCa                                                                                                        >  2:5087014/1‑150 (MQ=255)
                                                aCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCa                                                                                                        >  1:4948071/1‑150 (MQ=255)
                                                 cacggCAACCCGCGTTTCTCTCCGGGCCAAGGTGTGTTAGTGGCCAGCAATGGATGTGTGGTGACCGTTTCACCGCGCTAAATAAGTTGGGGTGTGGTCATTTTACTGACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCaa                                                                                                       <  2:7081277/147‑1 (MQ=255)
                                                  aaGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAAt                                                                                                      >  1:4450711/1‑150 (MQ=255)
                                                   aGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATg                                                                                                     >  1:9801834/1‑150 (MQ=255)
                                                     gCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGat                                                                                                   >  2:7292020/1‑150 (MQ=255)
                                                       aaCACTCGTTTCTCTCCGGGCCAACGTGTCGTCCTGGCCAGCACTGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGatat                                                                                                 >  1:1472239/1‑150 (MQ=255)
                                                          acTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGAtattct                                                                                              >  2:9019495/1‑148 (MQ=255)
                                                          acTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGAAACTCAATGAtattat                                                                                              >  1:6726316/1‑150 (MQ=255)
                                                              gTTTCTCTCGGGGCAACCGTGTGGTTCTGGCCTGCAATGGCTGTCTGGGTACCGTTTCTCCGCGCTAAATAAGTTGGGGTGTTGTCAGTTTACTGACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTc                                                                                          <  2:5966760/150‑1 (MQ=255)
                                                              gTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGTAACCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTc                                                                                          >  1:7340609/1‑150 (MQ=255)
                                                                  ctctcCGGGCCATGGGGTGGTACTGCCCCGAATTGCCTGTCCGTTTACCGTTTCTCCGCGCTAAATTGTTTGGGTTGTGTTCTGTTTACGCCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGt                                                                                      <  2:2977432/150‑1 (MQ=255)
                                                                   tctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGtt                                                                                     >  2:7185126/1‑150 (MQ=255)
                                                                       cGGGCCCACGGTTTCGTAGTGGCCAGCATTGTCCGTGCGGGGCCTTTTTCGCCGCTATAAATATGTTTGCGTGTGGACTGTTTCCTGCACACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTg                                                                                  <  2:3579334/150‑1 (MQ=255)
                                                                          gCCACCGTGTCGTACTGGCCAGCAAGGGCAGTCCGGTGACGGTTTCTCCGCGCTAAATATGTTTGGGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGAt                                                                              <  2:2406753/150‑1 (MQ=255)
                                                                            cAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGa                                                                            >  1:9669413/1‑150 (MQ=255)
                                                                             aaCGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAg                                                                           >  1:5934311/1‑150 (MQ=255)
                                                                                        tACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAg                                                                  >  2:8000982/1‑150 (MQ=255)
                                                                                         aCTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAGGTTTTGCGATGAGCGTAATCAGt                                                                 >  1:10109080/1‑150 (MQ=255)
                                                                                         aCTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGt                                                                 >  2:6832710/1‑150 (MQ=255)
                                                                                           tGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTg                                                               >  1:8239257/1‑150 (MQ=255)
                                                                                             gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCt                                                             >  1:4650868/1‑150 (MQ=255)
                                                                                             gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCt                                                             >  2:5834961/1‑150 (MQ=255)
                                                                                             gCCAGCAATGCCAGTCAGGTCACCGATTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAGGGCACCCTGTGGCCTGTAACTCCATGATCTTATGCTCCAGGGTTGCGATGAGCGTACTCAGTTGCt                                                             >  2:8826113/1‑150 (MQ=255)
                                                                                              ccAGCATGGGCAGTCTGGGGACCGTTACCCCGGGATAAATAGGTTTGCGTGTGGTCTGTTTACGGCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTc                                                            <  2:1058280/150‑1 (MQ=255)
                                                                                              ccAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCGGTGGCCTGTAACTCAATGATATTATGCTCAAGTGTTGCGATGAGCGTAATCAGGTGCTc                                                            >  2:5239457/1‑150 (MQ=255)
                                                                                               cAGCACTGGCTGTCTGGTGACGTTTTCTCCGCGCTAAATATGTTTGGGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTcc                                                           <  2:8315704/150‑1 (MQ=255)
                                                                                               cAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTcc                                                           >  1:518145/1‑150 (MQ=255)
                                                                                                aGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCa                                                          >  2:9549826/1‑150 (MQ=255)
                                                                                                  cAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGt                                                        >  2:5620103/1‑150 (MQ=255)
                                                                                                      ggCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCct                                                    >  2:8619575/1‑150 (MQ=255)
                                                                                                        cAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCctct                                                  >  2:6317374/1‑150 (MQ=255)
                                                                                                         aGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCGTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTg                                                 >  1:7966155/1‑150 (MQ=255)
                                                                                                           tCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCg                                               <  2:9352616/150‑1 (MQ=255)
                                                                                                             tggTGCCCGTTTCTCCGCGATAAATATGTTTGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGga                                             <  1:9963120/149‑1 (MQ=255)
                                                                                                               gtgccgttttCCCCGCTCTAAATATGTTGGCGTGTGGTCTGGTAACTCCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           <  1:1943287/143‑1 (MQ=255)
                                                                                                               gagcccGTTTCGCCGCGAAAATCACGTTGGCGTGGGGTGGGGTTACTCCCCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           <  1:2864151/146‑1 (MQ=255)
                                                                                                               gTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGgaga                                           >  2:773983/1‑150 (MQ=255)
                                                                                                                tgcccGTTTCTCCGCGCTAAATAAGTTGGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGAt                                          <  2:1723398/147‑1 (MQ=255)
                                                                                                                tGACCGTTTCTCCGCGATAAATAAGTTTGGGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGAt                                          <  2:7034792/150‑1 (MQ=255)
                                                                                                                    cGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGGGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGTTGCCa                                      >  1:9066176/1‑150 (MQ=255)
                                                                                                                    cGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCa                                      >  2:4222444/1‑150 (MQ=255)
                                                                                                                     gTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCGTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTTCGATGTGCGTAATCAGTTGCTCCAGTTCCTATGCGGAGATGCCAt                                     >  1:9466140/1‑150 (MQ=255)
                                                                                                                        tctcCGCGCTAAATTTGTTTGGTTGTGTTCTGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTa                                  <  1:7528312/150‑1 (MQ=255)
                                                                                                                            cgcgCTAAATATGTGGTCGTGTGTTAAGTTAACTGACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTACTCTGCGGAGATGCCATGTAatat                              <  1:10284982/150‑1 (MQ=255)
                                                                                                                                  aaaTATGTTTGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgc                        <  2:3237178/150‑1 (MQ=255)
                                                                                                                                  aaaTAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgc                        <  2:1626886/150‑1 (MQ=255)
                                                                                                                                    ttatttttGCGTGTGGTCAGGTTACTTACCACACGCCCCCTTCATTTCACCCTTGGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcgc                      <  2:6726316/145‑1 (MQ=255)
                                                                                                                                    atttttttGGGGGTGGTATGGTTAGTCACCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcgc                      <  1:3088032/145‑1 (MQ=255)
                                                                                                                                    aTAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTACTCTTTCcgcgc                      >  2:8352214/1‑150 (MQ=255)
                                                                                                                                      ttttttGGGTGTGGTCCGGTTAGTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCgg                    <  2:3851284/147‑1 (MQ=255)
                                                                                                                                       aGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCggg                   >  2:6597770/1‑150 (MQ=255)
                                                                                                                                           tgggtgtGGTCTGGTTAGTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGtttt               <  2:3125633/147‑1 (MQ=255)
                                                                                                                                           tGCGTGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGtttt               <  2:492347/150‑1 (MQ=255)
                                                                                                                                             ggtgtgtTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGt             <  2:518145/149‑1 (MQ=255)
                                                                                                                                                  tgGTCAGTTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         <  2:432161/150‑1 (MQ=255)
                                                                                                                                                  gggTCTGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat         <  1:3515300/149‑1 (MQ=255)
                                                                                                                                                    gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTACTATTTCCGCGCGGGTTTTGTtaataa       >  2:8118838/1‑150 (MQ=255)
                                                                                                                                                    gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa       >  2:3502964/1‑150 (MQ=255)
                                                                                                                                                    gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa       >  2:6166590/1‑150 (MQ=255)
                                                                                                                                                    gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa       >  2:9156903/1‑150 (MQ=255)
                                                                                                                                                     tCAGGTTACTGACCACACGCCCCCTTCATTTCAACCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGAGATGAGCGGAATCAGTTGCTCCAGTTCCTCTGCTGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAAc      >  1:70012/1‑150 (MQ=255)
                                                                                                                                                         gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  >  2:1580359/1‑150 (MQ=255)
                                                                                                                                                         gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  <  1:8348995/150‑1 (MQ=255)
                                                                                                                                                         gTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  >  2:4432417/1‑150 (MQ=255)
                                                                                                                                                         gTTACGGACCCCACGCCCCCTGCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCt  <  2:5335116/150‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
CTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCT  >  NC_000913/1720215‑1720513

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: