Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,438,536 A→C 6.5% L750R (CTG→CGG)  ydbK ← putative pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,438,5360AC 6.5% 1017.6 / 52.1 340L750R (CTG→CGG) ydbKputative pyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (156/162);  new base C (6/16);  total (162/178)
Fisher's exact test for biased strand distribution p-value = 7.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 349 total aligned reads displayed.

GTTCCAGTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTG  >  NC_000913/1438387‑1438684
                                                                                                                                                     |                                                                                                                                                    
   ccAGTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTAGACATGTTCCAGGCGAGACATCATATTGATGGATTTAAACTCTGGATTCTGACGCTCTTTCGCCGGGCAAACTTCGACGCACAGGt                                                                                                                                                   >  2:2183217/1‑150 (MQ=255)
    cAGTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGtt                                                                                                                                                  >  1:451839/1‑150 (MQ=255)
      gTTTGCTAGGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGAGATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTAc                                                                                                                                                <  1:2469986/150‑1 (MQ=255)
          gCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAAcc                                                                                                                                            <  2:1626732/150‑1 (MQ=255)
          gCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAAcc                                                                                                                                            >  1:1351035/1‑150 (MQ=255)
           cTACGGTCGATTTCTGGCAGGTTGAGTAAGAAATCGTAATTGATTTTCTCTTCTTCGACTTGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGGCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCg                                                                                                                                           <  1:1428116/150‑1 (MQ=255)
           cTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCg                                                                                                                                           >  2:2144360/1‑150 (MQ=255)
             aCGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGt                                                                                                                                         <  1:936422/150‑1 (MQ=255)
              cGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTATGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTg                                                                                                                                        >  1:778778/1‑150 (MQ=255)
               ggTCGATTTCTGGCAGGTTGAGTAAGAAATCGTAATTGATTTTCTCTTCTTAGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGAAGGTCTTTAGCCGGGCAAACTTCGACTCACAGGTTACAACCGCTGc                                                                                                                                       >  2:1402696/1‑150 (MQ=255)
               ggTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGc                                                                                                                                       <  2:696568/150‑1 (MQ=255)
                gTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGGTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcg                                                                                                                                      <  1:1196754/150‑2 (MQ=255)
                gTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCa                                                                                                                                      >  2:40084/1‑150 (MQ=255)
                gTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCa                                                                                                                                      >  1:1154427/1‑150 (MQ=255)
                gTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCa                                                                                                                                      <  2:1095916/150‑1 (MQ=255)
                 gcGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgt                                                                                                                                     <  2:2463388/149‑3 (MQ=255)
                  cGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtc                                                                                                                                    <  1:1114450/150‑4 (MQ=255)
                  cGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAAt                                                                                                                                    <  2:253802/150‑1 (MQ=255)
                   gATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtcg                                                                                                                                   <  1:2482999/150‑5 (MQ=255)
                   gATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATc                                                                                                                                   >  1:1080751/1‑150 (MQ=255)
                     tttCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtcgaa                                                                                                                                 <  1:2082629/150‑7 (MQ=255)
                     tttCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCCTTCGCCGGGCAATCTTCGTCGCACAGGTTACAACCGTTGCAATCtt                                                                                                                                 <  2:472111/150‑1 (MQ=255)
                      ttCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtcgaag                                                                                                                                <  1:1992535/150‑8 (MQ=255)
                      ttCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtcgaag                                                                                                                                <  2:2388840/150‑8 (MQ=255)
                       tCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTcc                                                                                                                               <  2:1892103/150‑1 (MQ=255)
                        cTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCg                                                                                                                              <  2:2318359/150‑1 (MQ=255)
                         tGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCgg                                                                                                                             <  1:54345/150‑1 (MQ=255)
                           gCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTg                                                                                                                           <  2:2265601/150‑1 (MQ=255)
                           gCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTg                                                                                                                           <  2:2072514/150‑1 (MQ=255)
                           gCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTg                                                                                                                           <  1:2160222/150‑1 (MQ=255)
                            cAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGc                                                                                                                          >  1:1685228/1‑150 (MQ=255)
                               gTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCAc                                                                                                                       <  2:279998/150‑1 (MQ=255)
                                  gAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTg                                                                                                                    >  2:2283826/1‑150 (MQ=255)
                                    ggAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCa                                                                                                                  <  1:272445/150‑1 (MQ=255)
                                    ggAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCCGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCa                                                                                                                  >  2:2090940/1‑150 (MQ=255)
                                       agaaATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGa                                                                                                               >  2:967525/1‑150 (MQ=255)
                                          aaTCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATc                                                                                                                                   >  2:613451/1‑127 (MQ=255)
                                               tAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTc                                                                                                       <  1:1792089/150‑1 (MQ=255)
                                               tAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTc                                                                                                       >  1:648076/1‑150 (MQ=255)
                                                 aTTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGa                                                                                                               >  2:204968/1‑140 (MQ=255)
                                                 aTTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTg                                                                                                     <  1:974736/150‑1 (MQ=255)
                                                  ttGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGc                                                                                                    >  2:1758475/1‑150 (MQ=255)
                                                   tGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGcc                                                                                                   >  1:2014937/1‑150 (MQ=255)
                                                   tGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGcc                                                                                                   >  2:774165/1‑150 (MQ=255)
                                                    gATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGccc                                                                                                  <  1:825123/150‑1 (MQ=255)
                                                     aTTTTCTCTTCTTCGACATGTCCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCcg                                                                                                 <  1:1824502/150‑1 (MQ=255)
                                                          ctctTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCat                                                                                            <  1:150740/150‑1 (MQ=255)
                                                           tcttctTCGACATGTTACAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACTGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATAATTCTGCCCGCGCata                                                                                           >  2:100354/1‑150 (MQ=255)
                                                            cttcttCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCatat                                                                                          >  2:332635/1‑150 (MQ=255)
                                                               cttcGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCAc                                                                                       >  1:1471064/1‑150 (MQ=255)
                                                                 tctaCATGTTCCAGTGGAGAGATGATATTGTTTGCTTTAATCTCTGGATTCTGAAGGTCTTTCTGCGGGGAAACTTCGACGCAAAGGTTACAACCGGTGCAATCTTCCGGTGGCACCTGCAAGAGATATTTCTGCCCGCTCATATCAcgc                                                                                     <  2:1848673/147‑1 (MQ=255)
                                                                  cGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCAcgcg                                                                                    <  1:2379642/150‑1 (MQ=255)
                                                                   gACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGAAGCTCTTTCGCCGGGCAAAATTCGACGCACAGGTTACAAACGGTGAAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCAGCGCATATCACGCGa                                                                                   >  2:1130667/1‑150 (MQ=255)
                                                                     cATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGAtt                                                                                 <  1:502213/150‑1 (MQ=255)
                                                                        gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGCCGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                              >  2:2443844/1‑150 (MQ=255)
                                                                        gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                              <  2:820737/150‑1 (MQ=255)
                                                                        gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                              <  2:16920/150‑1 (MQ=255)
                                                                        gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                              >  2:195906/1‑150 (MQ=255)
                                                                          tCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTcaca                                                                            >  2:44488/1‑150 (MQ=255)
                                                                           ccAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACAt                                                                           >  1:2330078/1‑150 (MQ=255)
                                                                              ggCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCa                                                                        >  1:2423519/1‑150 (MQ=255)
                                                                              ggCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCa                                                                        >  2:197848/1‑150 (MQ=255)
                                                                                  agaCATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGa                                                                    <  1:2091816/150‑1 (MQ=255)
                                                                                     catcatATTGATGGCTTTAATCTCTGGAATCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATg                                                                 >  1:1008753/1‑150 (MQ=255)
                                                                                       tcatATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCa                                                               >  2:2457557/1‑150 (MQ=255)
                                                                                          tatTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGc                                                            >  1:1882491/1‑150 (MQ=255)
                                                                                               aTGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCa                                                       <  2:1524459/150‑1 (MQ=255)
                                                                                                   cTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAgggg                                                   >  2:753781/1‑150 (MQ=255)
                                                                                                       aaTCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGtt                                               >  2:1500791/1‑150 (MQ=255)
                                                                                                         tctcTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGtttt                                             <  2:1455775/150‑1 (MQ=255)
                                                                                                           tctGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGtttt                                             <  2:1227024/148‑1 (MQ=255)
                                                                                                           tctGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTcc                                           <  2:243423/150‑1 (MQ=255)
                                                                                                            ctGGATTCTGCCGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCa                                          >  2:971855/1‑150 (MQ=255)
                                                                                                            ctGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCa                                          <  1:194857/150‑1 (MQ=255)
                                                                                                             tGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCAt                                         <  1:342064/150‑1 (MQ=255)
                                                                                                              ggATTCTGACGGTCTTTCGCCTTGCAAACTTAGACGTAGAGGTTACAACCGGTTCAATCTTCCTTTGCCACCTGCAAGAGATATTTCTGGCCGCGAATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATc                                        <  2:1216890/150‑1 (MQ=255)
                                                                                                                  tCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCtt                                    <  2:4585/150‑1 (MQ=255)
                                                                                                                    tGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCa                                  >  2:85951/1‑150 (MQ=255)
                                                                                                                    tGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCa                                  >  2:736157/1‑150 (MQ=255)
                                                                                                                     taCGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCGGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAg                                 <  2:1080751/149‑1 (MQ=255)
                                                                                                                      aCGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAgg                                >  1:2275542/1‑150 (MQ=255)
                                                                                                                      aCGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAgg                                <  1:1402696/150‑1 (MQ=255)
                                                                                                                       cGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAggc                               >  1:667065/1‑150 (MQ=255)
                                                                                                                        ggTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAggcg                              <  2:63910/150‑1 (MQ=255)
                                                                                                                           cTTTCGCCGGGCAAACTTCGACGCACAGGTTACGACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGca                           <  1:1875683/150‑1 (MQ=255)
                                                                                                                               cGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCCGcaccac                       >  1:1228284/1‑150 (MQ=255)
                                                                                                                                gCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACt                      >  1:1749467/1‑150 (MQ=255)
                                                                                                                                    ggCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTg                  >  1:2197356/1‑150 (MQ=255)
                                                                                                                                      cAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTgcg                >  2:490090/1‑150 (MQ=255)
                                                                                                                                       aaaCTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTgcgc               <  2:1333293/150‑1 (MQ=255)
                                                                                                                                         aCTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGa             <  2:1809410/150‑1 (MQ=255)
                                                                                                                                          cTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAATACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGTCTGGCAGGGGCGTTTTCCATCGCTTCAGGCTGCACCACTTTTGCGCGaa            >  1:1223728/1‑150 (MQ=255)
                                                                                                                                          cTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGaa            >  1:1740743/1‑150 (MQ=255)
                                                                                                                                          cTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGaa            <  1:1890391/150‑1 (MQ=255)
                                                                                                                                            tCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATa          <  1:2457557/150‑1 (MQ=255)
                                                                                                                                            agattgCACCGGTTGTAACCTGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATa          >  1:735978/6‑150 (MQ=255)
                                                                                                                                            agattgCACCGGTTGTAACCTGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATa          >  2:1337623/6‑150 (MQ=255)
                                                                                                                                                  cacaGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGGTTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAg    <  1:1604023/150‑1 (MQ=255)
                                                                                                                                                   acaGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAgt   >  2:589035/1‑150 (MQ=255)
                                                                                                                                                    caGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAgtg  >  1:1289576/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTCCAGTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTG  >  NC_000913/1438387‑1438684

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: