Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 140,499 T→C 5.3% S243G (AGC→GGC)  gcd ← quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913140,4990TC 5.3% 291.1 / 4.7 94S243G (AGC→GGC) gcdquinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (41/48);  new base C (2/3);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

GAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCA  >  NC_000913/140352‑140645
                                                                                                                                                   |                                                                                                                                                  
gaaTGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTgg                                                                                                                                                  >  2:1216494/1‑150 (MQ=255)
   tGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTggcgg                                                                                                                                               <  1:70634/150‑1 (MQ=255)
    gATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTggcggc                                                                                                                                              <  2:447823/150‑1 (MQ=255)
    gATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCCagcggc                                                                                                                                              <  1:1309640/150‑1 (MQ=255)
      tACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCAt                                                                                                                                            >  2:634150/1‑150 (MQ=255)
       acgacgCGGGAAATCCGCCATCACTTCCGGCGAAGCGGTTTCTCCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATc                                                                                                                                           >  2:911190/1‑150 (MQ=255)
        cgacgCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCa                                                                                                                                          <  1:655592/150‑1 (MQ=255)
        cgacgCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCa                                                                                                                                          >  1:1030930/1‑150 (MQ=255)
         gacgCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCaa                                                                                                                                         >  1:310359/1‑150 (MQ=255)
         gacgCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCaa                                                                                                                                         >  1:725287/1‑150 (MQ=255)
         gacgCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCCAGCGGCaaaga                                                                                                                                         <  2:350006/150‑5 (MQ=255)
           cgcgGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAgc                                                                                                                                       <  2:810745/150‑1 (MQ=255)
           cgcgGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGGTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGAGCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAgc                                                                                                                                       <  2:58544/150‑1 (MQ=255)
             cgGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAgcgc                                                                                                                                     >  1:1209289/1‑150 (MQ=255)
              gggCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCa                                                                                                                                    >  2:1167007/1‑150 (MQ=255)
               ggCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCaa                                                                                                                                   >  1:1215695/1‑150 (MQ=255)
                 cAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAAc                                                                                                                                 >  2:1286895/1‑150 (MQ=255)
                      cGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGgcg                                                                                                                            >  2:288977/1‑150 (MQ=255)
                              cTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAg                                                                                                                    >  1:38314/1‑150 (MQ=255)
                               ttCCGGCGAAGCGGTTTCTGCTTTGGCGTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGc                                                                                                                   <  2:1460478/150‑1 (MQ=255)
                                tCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTg                                                                                                                      <  1:911978/146‑1 (MQ=255)
                                tCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTg                                                                                                                      >  2:911978/1‑146 (MQ=255)
                                tCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCg                                                                                                                  <  1:1286895/150‑1 (MQ=255)
                                 ccGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCgg                                                                                                                 <  2:856223/150‑1 (MQ=255)
                                 ccGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAACAGACCCCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCgg                                                                                                                 >  2:74454/1‑150 (MQ=255)
                                   ggCGAAGCGGTTGCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTa                                                                                                               <  1:1487980/150‑1 (MQ=255)
                                    gCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTac                                                                                                              <  2:1389813/150‑1 (MQ=255)
                                    gCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTac                                                                                                              <  1:528116/150‑1 (MQ=255)
                                      gAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCATGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTATCTTTGCCGCTGGCGGCATCAAGAGCAAACAGGCTCTGGTGAGCGGTacac                                                                                                            <  2:636087/150‑1 (MQ=255)
                                      gAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTacac                                                                                                            >  2:742788/1‑150 (MQ=255)
                                      gAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCCAGCGGCaa                                                                                                                                            <  1:1114325/118‑2 (MQ=255)
                                      gAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCCAGCGGCaa                                                                                                                                            >  2:1114325/1‑117 (MQ=255)
                                        aGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAg                                                                                                          >  2:626954/1‑150 (MQ=255)
                                        aGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAg                                                                                                          >  2:568770/1‑150 (MQ=255)
                                         gCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAgg                                                                                                         <  2:419679/150‑1 (MQ=255)
                                           ggTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTa                                                                                                       <  1:533252/150‑1 (MQ=255)
                                               tCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAgg                                                                                                   <  2:1209289/150‑1 (MQ=255)
                                                cTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAggg                                                                                                  >  1:452170/1‑150 (MQ=255)
                                                cTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAggg                                                                                                  >  2:388676/1‑150 (MQ=255)
                                                   cTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCATTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGgtgt                                                                                               >  1:1447936/1‑150 (MQ=255)
                                                     ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCg                                                                                             <  1:215496/150‑1 (MQ=255)
                                                        gCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTTAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGccc                                                                                          <  2:144651/150‑1 (MQ=255)
                                                        gCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGccc                                                                                          <  2:223513/150‑1 (MQ=255)
                                                           tCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACt                                                                                       <  1:906241/150‑1 (MQ=255)
                                                             aTGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCACACAGGCGCTGGTGAGCGGTACACAGGTTAAGGGTGTCGccccgttt                                                                                     >  2:1286707/1‑145 (MQ=255)
                                                             aTGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACttt                                                                                     <  1:634150/150‑1 (MQ=255)
                                                             aTGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACttt                                                                                     <  2:420562/150‑1 (MQ=255)
                                                              tGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTa                                                                                    >  1:367876/1‑150 (MQ=255)
                                                               gATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCAGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTaa                                                                                   >  2:1125147/1‑150 (MQ=255)
                                                                 tAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATc                                                                                 >  2:180860/1‑150 (MQ=255)
                                                                  aaGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCg                                                                                >  2:1159330/1‑150 (MQ=255)
                                                                   agagACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCgg                                                                               >  2:415553/1‑150 (MQ=255)
                                                                        caccacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCa                                                                          >  1:897877/1‑150 (MQ=255)
                                                                              ggCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCAt                                                                    >  1:864411/1‑150 (MQ=255)
                                                                              ggCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCAt                                                                    <  2:1215695/150‑1 (MQ=255)
                                                                                 aGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTgg                                                                 <  1:847525/150‑1 (MQ=255)
                                                                                  ggTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGt                                                                >  1:41664/1‑150 (MQ=255)
                                                                                    ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGa                                                              >  1:121790/1‑150 (MQ=255)
                                                                                    ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGa                                                              <  2:888236/150‑1 (MQ=255)
                                                                                    ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGa                                                              <  1:568770/150‑1 (MQ=255)
                                                                                         tGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCa                                                         >  2:257404/1‑150 (MQ=255)
                                                                                          gCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCAc                                                        >  2:1143804/1‑150 (MQ=255)
                                                                                             ggAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCg                                                     >  2:525076/1‑150 (MQ=255)
                                                                                             ggAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCg                                                     <  1:204690/150‑1 (MQ=255)
                                                                                              gAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCgg                                                    >  1:999252/1‑150 (MQ=255)
                                                                                              gAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCgg                                                    <  2:1069554/150‑1 (MQ=255)
                                                                                               aaaGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGa                                                   <  1:388676/150‑1 (MQ=255)
                                                                                                aaGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGAt                                                  <  1:1159330/150‑1 (MQ=255)
                                                                                                      cGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCg                                            <  1:282460/150‑1 (MQ=255)
                                                                                                       gTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCTCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCgg                                           <  1:1216494/150‑1 (MQ=255)
                                                                                                              ttCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAAGAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTc                                    <  1:800126/150‑1 (MQ=255)
                                                                                                              ttCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTc                                    <  1:180860/150‑1 (MQ=255)
                                                                                                                cAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACATGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATc                                                 >  2:780480/1‑135 (MQ=255)
                                                                                                                cAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACATGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATc                                                 <  1:780480/135‑1 (MQ=255)
                                                                                                                  gCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACAt                                >  1:1494269/1‑150 (MQ=255)
                                                                                                                   cTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATc                               <  2:1030930/150‑1 (MQ=255)
                                                                                                                    tCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCg                              <  1:415553/150‑1 (MQ=255)
                                                                                                                      aGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGcc                            <  1:626954/150‑1 (MQ=255)
                                                                                                                        gATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAg                          <  1:1167007/150‑1 (MQ=255)
                                                                                                                          tCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTa                        >  1:1077637/1‑150 (MQ=255)
                                                                                                                            gTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACg                      >  2:1064966/1‑150 (MQ=255)
                                                                                                                             tAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAAGAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCGGACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACgg                     <  1:911190/150‑1 (MQ=255)
                                                                                                                             tAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCTGCTGCTTCACATCGCCAGTACgg                     <  2:795133/150‑1 (MQ=255)
                                                                                                                                tGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTTAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGCac                  >  1:972006/1‑150 (MQ=255)
                                                                                                                                tGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCGCACTTTAATCGGCGTCACTTCATTGGTGATTTAACCCGGATCGTTAGGCTGCTTCACATCGCCAGTACGGAac                  <  2:446044/150‑1 (MQ=255)
                                                                                                                                        cttgatGCCGCCAGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGACCACCCAGGc          >  1:344946/6‑150 (MQ=255)
                                                                                                                                        ctctTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGc          <  2:864411/150‑1 (MQ=255)
                                                                                                                                        ctctTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGc          >  2:1232158/1‑150 (MQ=255)
                                                                                                                                         tctTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGct         <  1:742788/150‑1 (MQ=255)
                                                                                                                                              gCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGACCACCCAGGCTTCttt    >  1:1105216/1‑150 (MQ=255)
                                                                                                                                              gCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCttt    <  1:1272931/150‑1 (MQ=255)
                                                                                                                                              gCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGGGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCttt    <  2:310359/150‑1 (MQ=255)
                                                                                                                                                cGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCa  <  1:593869/150‑1 (MQ=255)
                                                                                                                                                cGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCa  <  2:1235609/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAATGATACGACGCGGGCAATCCGCCATCACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTTGCCGCTGGCGGCATCAAGCGCAAACAGGCGCTGGTGAGCGGTACACAGGTAAAGGGTGTCGCCCACTTTAATCGGCGTCACTTCATTGGTGATTTCACCCGGATCGTTCGGCTGCTTCACATCGCCAGTACGGAACACCCAGGCTTCTTTCA  >  NC_000913/140352‑140645

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